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ABSTRACT 


The physical mechanisms of the various processes which may occur as a 
result of encounters between low velocity atoms or ions are discussed, 
particular attention being paid to the factors which determine whether 
or not a collision is likely to be near-adiabatic. The following are the 
processes considered : 

excitation, A+B+A+B’, At++t B+At+ B’, At+B+At'+B; 

charge transfer, At+t+B>A+Bt, A-+B+A+B-, A+B>At+B-, 
Att++ B+At+ Bt ; 

tonzation, A+B>A+Bt+e, At+B>At+ Btte, At++BsAtt++ Be; 

detachment, A-+B+>A+B+e, A-+B+AB-+e. 


§ 1. InTRODUCTION 


EXCITATION and ionization of atoms or ions by the impact of other atoms 
or ions may be discussed in a similar way to excitation and ionization by 
the impact of electrons provided the ratio of the relative velocity of the 
colliding particles to the relevant orbital velocity is high. When this 
ratio is low the situation is entirely different for the transitions are then 
only energetically possible if the encounter is between two massive systems. 
Such inelastic encounters are of interest in many connections. They are, 
for example, responsible for the luminosity and other effects resulting from 
the passage of a meteor through the earth’s atmosphere (cf. Herlofson 
1948). Again it is important to know whether a significant contribution 
to ion production in a discharge can arise from slow collisions between the 
positive ions and the gas atoms (cf. Llewellyn Jones 1953) ; and should 
negative ions be present to know whether they can suffer appreciable 
destruction through detachment taking place in slow collisions (cf. Massey 
1950). Many technical applications require the production of intense 
ion beams. One of the limitations on the intensity which may be achieved 
is the dispersing effect of space-charge. This effect may be neutralized to 
some effect at least by the electrons liberated in collisions with atoms of the 


* Communicated by the Authors. 
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residual gas. Information on the dependence of the efficiency of such 
collisions on the relative velocity, and on the nature of the ions and atoms 
involved, is necessary in order to examine the conditions under which 
self-neutralization is likely to be adequate. 

In a first attempt at classification of the possibilities it has been custo- 
mary to refer to slow collisions as near-adiabatic. This implies that they 
occur so gradually that the chance of any ultimate transfer of energy is 
slight. There is experimental evidence that in some cases the associated 
cross sections are indeed small, and fall off rapidly with decrease of the 
relative velocity (cf. Massey and Burhop 1952). However, this does not 
seem to be truein general. Although the results obtained in the laboratory 
by different investigators are sometimes in conflict they leave little doubt 
that many slow collisions have not these characteristics (cf. Massey and 
Burhop 1952). The same conclusion is indicated by the remarkable 
effectiveness with which ion beams in a good vacuum (10-° mm Hg) can 
neutralize their own space charge. Further the light and ionization 
produced by meteors appears to be much more intense and much less 
dependent on the velocity than would be expected if all slow collisions 
were near-adiabatic (Opik 1953, Kaiser 1953). 

Progress in the study of slow collisions clearly requires elucidation of the 
reasons why they are not always near-adiabatic and recognition of the 
factors which determine whether or not they can be so classified. The 
present paper represents an exploratory effort in this direction. After 
noting the predictions of Born’s approximation (according to which all 
slow collisions are near-adiabatic) we consider the effects arising from the 
crossing of the various potential energy curves (including radiationless 
transitions into the continuum). Charge transfer, on which reliable 
experimental data is now becoming available (Hasted 1951, 1952) is 
amongst the collision processes discussed. No attempt is made at 
carrying through precise quantal calculations for any specific case. 


§ 2. Born’s APPROXIMATION 


Results based on Born’s approximation cannot of course be expected to 
be accurate at low velocities but nevertheless they provide a useful basis 
for discussion. Consider a system A of effective nuclear charge Z,e 
incident on a system B which, for simplicity, will be taken to be comprised 
of a nucleus of charge Ze with an electron in the 1s state. By expanding 


them suitably it may be shown (Bates and Griffing 1953) that in the low — 


energy region the Born approximations to the cross sections for the 
excitation of B to the 2s and 2p states reduce to 


((1s—2s)=={1:3 x 10-°Z .2H4/ M4 AE(1s—2s)®}rra,2, :. i Ge 
Q(1s—2p)={7-4 x 10°°Z ,2H5/M5 AE(1s—2p)7}na,2 ; 2) Eee 


and similarly that the Born approximation to the cross section for ioni- 
zation reduces to 


and 


Q(1s—c)={1-9 x 10-"Z 24) M4 AE (1s—c)®\na,2, ar ae 
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where a is the radius of the first Bohr orbit (so that Ty" iS 8-8 X 10-1” cm?) ; 
# is the impact energy and AEF is the energy change, both inev ; M is the 
mass of the incident system relative to the mass of the proton. It is 

-apparent from the formulae and from table 1, which gives numerical 
results for some representative cases, that the cross sections are rapidly 
increasing functions of H ; that they are rapidly decreasing functions of 

AF and of M ; and that they are greater for excitation than for ionization 
but are extremely small for both. These features are typical of near- 
adiabatic collisions, that is, collisions whose effective duration, 7’, is long 
compared with the time period, 7, associated with the internal motion 
(cf. Massey 1949). It is usual to take 


= la/v: ct) AE, ee ee) 


where / is a number of the order of unity measuring the range of the 
interaction, and v is the velocity of approach. On this basis the near- 
adiabatic condition Tir >I, he eer 


becomes — 6LABM12/H12 51, ry te ce (0) 


“energy and mass being in the same units as before. For all cases treated 
in table 1 7/7 is 21 or greater. 


Table 1. Data on Representative Near-adiabatic Collisions 


Incident particle :* mass, M=30 
square of effective nuclear charge, 
(Zuejt—10e" 


Excitation Energy H (ev) 
eas or ionization. 
poeEe ae 100 200 400 800 
(ev) 
itati Cross section (7r”) 
Excitation 0 
ls—2s 2 Ole 4 x 102 |. 63,2010 Px 10=* 
3 Deel erles x 105 | 0-5 L059, x LO? 
4 Alm 6 coe l0=e el clas iki 10 A 
= =3 
ls—2 2 Deel ae nite Omelet SLO Pea LO 
oe 3 1, 10-8 AOC a Mala X10? | 4s x 10s 
4 eeet0- | oe x 10g etl 10 oH) Cy x10 
‘Tonization 7»: Fe 
ls—c Bale So Om gel) 27, «10 


6 

8 9-, 10-20] 1-,x 10-8 | 2-3x10-7 | 3-,x 10-8 
10 2+, 10-29] 3+,x 10-8 | 69x 10-8 | 9410-7 
12 7-,X10-12| 1-3X10-9 | 2%)X10-8 | 3-2 10-" 


* Here, and in the later tables, results for other cases may easily be obtained 
by changing M and ZH by the same factor. 
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As mentioned in § 1, though many experimentally determined cross 
section versus energy curves have the characteristics described above 
others do not even well within what might at first be thought to be the 
near-adiabatic region. There is no real anomaly. As can be seen most 
clearly from the perturbed stationary state treatment the criterion derived 
from the definitions of 7’ and 7 given in (4) is much too crude (Bates, 
Massey and Stewart 1953). Briefly, it is not permissible to regard A and 
B as separate atoms or ions ; they must instead be regarded as forming a 
quasi-molecule and account must be taken of the dependence of the various 
parameters on R, the distance between the nuclei. Because of this it is 
not possible to express the near-adiabatic condition simply and compre- 
hensively. However, in many cases it is sufficient to modify (6) by 
replacing 4U( co) (for which the designation 4H was used) by 4U(Rj,) the 
minimum separation between U,(R) and U,(R&), the initial and final 
potential energy curves ; and by replacing / by fl where f is the fraction 
of the range of interaction over which the separation is near the minimum 
value. The effect may of course be slight. It will often happen that 
U(R) and U,(R) come closer together as the systems approach but that 
AU(Ry,) and fl remain large enough for the near-adiabatic condition to be 
satisfied even at quite high energies. In this event a detailed perturbed 
stationary state treatment would yield much the same type of result as 
would the Born approximation, the main difference being that the cross 
section would become appreciable at rather lower energies. Much more 
interest is attached to the possibility of collision processes for which 
AU(fR),) and fl are so small that the near-adiabatic region is virtually 
non-existent ; and in considering this possibility special attention must 
obviously be paid to the pseudo-intersection of potential energy curves. 


§ 3. TRANSITIONS AT AN INTERSECTION OF PoTENTIAL ENERGY CURVES 
Let U,°(R) and U,°(R) be the zero-order approximations to two potential 
energy curves, U ,(#) and U,(R), of the same electronic species. As is well 
known U,° and U,® may cross at some internuclear separation, R,dp, 
but U; and U, are prevented from so doing by the mutual interaction 
(von Neumann and Wigner 1929). However, if the velocity of relative 
motion, v, is finite there is a certain probability, P, that in passing through 
f,, the colliding systems transfer from one potential energy curve to the 
other ; and taking account of the fact that passage through R,, must be 
made twice, there is therefore a probability 
PaeP l= hy eee ee see eee (7) 
that systems which originally approach along U, finally recede along U,, 
one or both of them being excited. Clearly the cross section for the 
process is approximately 
Q= (GFR ad sine et aerate or ke em ees 
where g is the ratio of the statistical weight associated with U ; to the sum 


of the statistical weights associated with all possible initial potential energy 
curves. 
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Landau (1932), Zener (1932) and Stueckelberg (1932) (whose treatment 
of the whole problem is the most complete) have shown that 


Po — 

ae OX Die 0) aumee euiaeene Sa eta s. r.4¥ (9) 
1a 
5={ 7] 402/20 dn a) eee een LO) 
«z 
which may be expressed more conveniently as 

d 

= {oun| ave [ie ap (Ue—U,;") is aie oF (LL) 


the mass of the incident particle being measured relative to that of the 
proton, the energies in ev and the length in atomic units as usual. 


3.1. Hacitation 
In simple excitation 

eee er eee eee, oo oo (19) 
At++ B>A++B’, |) 
AMR SAY EB 
crossing obviously can only occur at small or moderate internuclear 
separations where the potential energy curves are distorted from their 
asymptotic form. Consequently the interaction is in general likely to be 
strong, making 4U(&,) comparable with the distortion. It is important 
to note that this is not always the case. For instance there may be 
chance cancellation within the relevant matrix integral ; and again vio- 
lation of the selection rules may cause 4U(R,) to be small but yet not 
vanishingly so if the quantum numbers are poorly determined. As 
regards (d/dk)(U ;°—U,°)p_p,, little can be said beyond that in many 
instances it must be several ev per atomic unit oflength. To allow for the 
uncertainties the probability, 7, was computed from formulae (7), (9) and 
(11) for various choices of the parameters. The values obtained are 
presented in table 2. It will be observed that over quite a wide range, P 
(and hence the cross section), is large and is a slowly varying function of 
the energy of the incident particle. This suggests that very marked 
departure from the characteristic near-adiabatic behaviour is not unduly 
rare when crossing occurs (though it must frequently be prevented by the 
interaction being too strong). Unfortunately the information available 
on the potential energy curves of even the better known systems is meagre, 
and it is not easy to predict in any particular case whether there is, or 
there is not, crossing. In general, however, the phenomena would seem to 
be most likely if the colliding systems have a number of low lying levels 
and least likely if they have only high levels : that is, it would be expected 

to arise in processes such as 


O+Fe+0-+ Fe’ , shah ieee eee ees Cla) 


rather than in processes such as 
H--He>H+He’.-". ee. ee) 


(13) 
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3.2. Charge Transfer 


Collisions of the types 
At++B+>A+B*, rs ep 


A*-BSA+B, 2 


Table 2. Data on Representative Curve-crossing Collisions Yielding 
Excitation 


Incident particle : mass, 1J=30 


a U4 T,9 energy, E (ev) 
| 4u(R,)| ||an Ce bt) Per, 
ev ev per atomic unit 200 400 

of Jength 


Probability of excitation 
3°9X10-1 | 3-2«10-1 | 2:5x10-1 | 1910-1 
2°510-1 | 1910-1} 1-4x1071 | 1-0 10-1 
1-4x10-1 | 10107! | 7-4x10-2 | 5-3 10-2 
74x 10-7 | 5-3 10-* | 3-6X10— | 2-710] 


rare 


4-1X 10-1 | 4:9x10-1 | 5-0x10-! | 46x10-19 
5010-1 | 4-6 1077 | 3-9 10-1 3-210" 
39x 10-4 |. 3:2x10-* | 2510-4 | 1-9 se Gee 
2-5 10-1 | 1:910-1 | 1-410" | LOsc1G= 


1 
2 
4 
8 


1310-2 | 5-610-2 | 1:510-1 | 2-8 10-1 9 
1-5 10-1 | 2-8x10-1 | 4-1x 107-1 | 4:9 x 10-2 
4-1x10-1 | 4:9x10-1 | 5-0x10-1 | 4-6x 10-1 
5:0x10-1 | 4-6x10-1 | 3-9 10-1] 3-2 10-2 


Oreo 


4-1X10-® | 1-4x10-® | 9-1x10-5 | 1-7 10-3 
9:1%10-5 | 1-710 | 1-3 10-2 oe 
1:3 10-? | 5-6 10-2 | 1510-1 | 2-8 10-1 
1510-1 | 2-8x 10-1 | 4:1x10-1 | 4:9 10-1 


OPW rR 


<10-10 < 10-10 <10-10 < 10-10 

<1 Q>520 <10-10 4:1x10-° | 1-4 10-8 
4-1x10-® | 1-4x10-6 | 9-1 «10-5 | 1-7 10-3 
9-1 x 10-5 | 1:7 10-8 | 1-3 10-2 | 5-6 10-2 


ore 


are essentially similar to excitation collisions and hence do not require 
separate discussion. This is also true of many of the collisions 

Aq BRAT BR os ee Sic dS 
but in others the Coulomb attraction between the positive and negative 


ions formed causes crossing to occur so far out that special considerations 
enter: thus for the endothermic processes* 


* The electron affinities and ionization potentials used throughout are taken 
from Massey (1950) and from Moore (1951). 
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Ke Br Ke Bre (0-8 ey). 2 4 os 719) 

Na+Cl+Nat+Cl (=1-4ev), . . . . . (20) 

Al+0+Al++0-  (—3-8ev),_. (21) 

f, is about 35, 20 and 7 respectively. Much the same remarks apply to 
the collisions At++-_B>At-+_ Bt (22) 


as here again distant crossing may occur : thus for the exothermic processes 
Sit+++N->Sit+N+  (+1-8 ev), (23) 
At++He+At+Het (+3-0 ev), (24) 

f,, is about 15 and 9 respectively. 
In cases such as those cited above (and even when R, is rather less) it is a 


sufficient approximation to represent the potential energy curves by their 
asymptotic forms so that (11) may be written 


8=2-5 x 102M V2 AU(R,)?/B¥2AU( «0? (25) 
which, when substituted in (7) and (9), gives the probability F in terms of 
a convenient set of parameters. Some computations were carried out to 


Table 3. Data on Representative Curve-crossing Collisions Yielding 
Charge Transfer and Involving a Coulombic Potential Energy Curve 


Incident particle : mass, M/=30 


energy, H (ev) 


AU(R,) 
AU(c) 100 200 | 400 | 800 
Probability of charge transfer P 

0-05 A104 3-4 10-4 2-6 x 107? 2-0 x 10-2 
0-10 ator < 1h Oa COE as DO On* 4-7 x 10-1 
0-15 9036104 ax 10s4 3-4 10-1 4-5x 10-1 
0-20 862910 4-210? 12 10=* 2510-4 
0-25 AT x10 49x10 228x102 TO x 104 
0-30 1-0 105% ool RO: 4-4 10-3 2-6 x 10-2 
0-35 Te25cl0-4 eo LO @ 4-8 x 10-4 5-5 x 10-3 
0-40 TS yc 1 OR 200 eds e107? 3-8 x 10-° 2:0 <10=* 


enable the general position to be assessed easily. From the results of 
these, presented in table 3, it can be seen that F is a slowly varying 
function of the impact energy provided 
AU(B,,)/4U(«0)<~0°15. 
It may be observed that if 4U( co) is 3 ev crossing occurs at an internuclear 
separation of 9-Ia) and for condition (26) to be satisfied 4U(R,) must be 


less than about 0-45 ev; if JU(c) is 4 ev the corresponding values are 
about 6-8a, and 0-60 ev ; and if it is 5 ev they are 5-4a, and 0-75 ev. 


(26)* 


* Ag usual MV is taken as 30 for the sake of definiteness. 
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The figures just quoted suggest that the condition is quite severe as far 
as process (18) is concerned, for it appears likely that in many instances 
AU(R,,) must exceed the relevant limit.* However, it must be borne in 
mind that a number of potential energy curves are usually associated with 
a given pair of systems in given states ; and that it is sufficient for the 
interaction to be less than the appropriate limit at only one crossing point. 
Some of the possible effects have already been mentioned in connection with 
excitation collisions. Special attention may be drawn to collisions, such 
as (21), in which a p (or d) electron is transferred from one system to the 
other. If 4U(c) is, say, 4ev 4U(R,) may well exceed the limit at the 
crossings associated with the o orbitals but be less than the limit at the 
crossings associated with the 7 orbitals ; for the overlap integral (which 
gives some measure of the interaction}) is in general very much smaller 
for the latter than it is for the former. 

In the case of processes (22) condition (26) is unlikely to exert so 
restrictive an influence since the greater compactness of the wave functions 
tends to make the interaction weaker. 

If § is small compared with unity formulae (7), (8), (9), (11) and (25) may 
be combined to give 


M24U(R,)? 
sf ( zr} naa? (27) 


ca {30x 10°" PEZU(co)* 


the notation and units being as before. Some numerical values are given 
in table 4. 

Magee (1940) has made estimates of the interaction energy between each 
of the various alkali metal-halogen pairs at the (A, B)—(A*, B-) crossing 
points. He finds that when £,, is 20 4U(R,) is usually at least 2x 10-2 ev 
and that when R, is 30 it is usually at least 1x10-%ev. From these 
results and the entries in table 4 it would appear a type (18) collision which 
is endothermic by between perhaps 1 and 3 ev is quite likely to have a 
large and slowly varying cross section.t Because of the increased prob- 
ability of the interaction being too strong the chance of this is less for 
one which is endothermic by a much greater amount. 

Calculations on the interaction energy at (A++, B)—(A+, Bt) crossing 
points have not been performed but in view of the compactness of the 


_ * It is of course certainly satisfied if JU(0o) is sufficiently small because R 
is oe large and hence, as the interaction falls off exponentially, AU(R,) is 
minute. 

_ } Indeed Magee (1940) proposes that if the transition is allowed, and R, 
is not too extreme, 4U(Rz) may usually be taken to be of order (10S8/R,) ev 
where S is the overlap integral between the initial and final wave functions of 
the active electron. Tables of overlap integrals have been published by 
Mulliken, Rieke, Orloff and Orloff (1949). 

{ The two limits given here should not be treated as other than illustrative. 
It is of course impossible to formulate any precise general rule since 4U (R,) 
depends not only on &, but also on other factors such as the orbitals involved 
and the energies associated with each (which control the amplitudes of the 
wave functions in the region of most interest) 
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relevant wave functions there is no doubt that at the corresponding 
internuclear separations 4U(R,,) is much lower than the values obtained 
by Magee. Consequently a type (22) collision process will almost certainly 
have a small cross section if it is only slightly exothermic. One which is 
exothermic by rather more than perhaps about 2 ev* is quite likely to have 
a large and slowly varying cross section ; but it is less probable that one 
which is very exothermic has such a cross section. 


Table 4. Further Data on Representative Curve-crossing Collisions 
Yielding Charge Transfer and Involving a Coulombic Potential 
Energy Curve 

(6< 1) 


Incident particle : mass M=30 
energy H=100 ev 


Internuclear Cross section per 

Energy change, 4U(co) | separation at unit statistical 
(ev) crossing, Ry weighting ratio, Q/g 

(4) (7ra19”) 

1-0 20-2, DO ORAL (ie) 
1-5 18-1 30610) PAU Ra) 
2-0 13-6 Tied See eal) (1b,.) 
2:5 10-9 50x10? AU(R,)* 
3-0 Veil 2-410 AU(R,)? 
3-5 7:8 ox dO aU Ch, )e 
4-0 6:8 to O2ge een lone 


(with 4U(R,) in ev) 


It-is perhaps worth drawing specific attention to the fact that a few 
collision processes of type (18) are exothermic, and many of type (22) are 
endothermic. Most of these should exhibit near-adiabatic characteristics 
since curve-crossing must be very rare. The former group does not include 
any important cases but the latter group includes 


Cat++O-+Cat+O+r (—1-7 ev), ee ate (5) 
and several others which are of geophysical or astrophysical interest. 


§ 4, TRANSITIONS INTO THE CONTINUUM 


4.1. Ionization 
Tonizing collisions Mey BAR Bice ale aries ere 29) 


can take place as a result of curve-crossing but they must be treated rather 
differently from excitation or charge transfer owing to the fact that the 


final state is in the continuum. 
Suppose that at an internuclear separation of FR, there is an inter- 
section of one of the potential energy curves along which A and B may come 


* Once again a definite limit cannot be given. 
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together and one of the potential energy curves of the complex comprised 
of A, Bt and an ejected electron of zero energy. At smaller internuclear 
separations radiationless transitions may occur, that is, the quasi-molecule, 
(A+B), may suffer auto-ionization. If tis the lifetime towards this 
process (measured in sec) the cross section associated with (29) is 
approximately 
Q=({gR,,2[1—exp (—R,a,/vt)]}wag2= . - - - - + - 680) 
= (gR,>[1—exp (—3-8 X 10-5 M1?R,/B¥t)}}aa?_. (31) 
(provided of course EF is sufficiently high to permit the necessary close 
approach). In many instances the exponential coefficient in (31) is small 
so that the formula simplifies to 
Q={3-8x 10-Mg MAAR AEA ra) > 
It is clear from the above that cross sections for ionization (or simul- 
taneous ionization and excitation) may be slowly varying functions of the 
impact energy ; and moreover, since in favourable cases t may be as short 
as 10-15 sec (Shenstone 1931, 1948, Wu 1944, Bransden and Dalgarno 
1953) they may be quite large. As in the case of simple excitation reliable 
predictions on the prevalence of curve-crossing are difficult to make. 
It may, however, be noted that because of the character and position of the 
final potential energy curves the phenomena must be rare in collisions of 


the types At+5BSAt EB ve, a tae aca 

At++ B+>A+++B-+e, » oo, el oo) ee 
so that a departure from near-adiabatic behaviour would not be expected 
in either. In certain cases the sequence of process (16) followed by (29) 


might of course make a slowly moving positive ion an indirect but quite 
effective liberator of electrons. 


4.2. Detachment 
Finally we will briefly consider detachment 
A-~+-B+A+ B46... se Le 


This may of course occur in the manner just described. In many instances 
the crossing point may be reached even in collisions of thermal energies.* 
For example the theoretical work of Eyring, Hirschfelder and Taylor 
(1936) on the structure of H,~ shows that H atoms and H- ions approach 
along an attractive potential energy curve which intersects the H, 
ground state potential energy curve at an internuclear separation of about 
34a. Again, though proper calculations have not been carried out it 
seems likely (Bates and Massey 1943) that several of potential energy 
curves along with O atoms and O- ions approach are algo attractive and 


intersect the lower O, potential energy curves at comparable internuclear 
no teint tee LS Te te 

* At such energies the reaction path is A~+B—> AB-+e. 

+ This is naturally also possible in de-excitation collisions and in exothermic 
charge transfer collisions. 


Slow Inelastic Collisions between Atomic Systems 121 


Separations. Such thermal processes (which are referred to as associative 
detachment) may be important in the upper atmosphere and in discharges. 
Because of lack of knowledge of ¢ it is not possible to give a reliable estimate 
of the cross section in any particular case; but it is apparent from 
formula (30) (modified to allow for the fact that the variation of the 
_ kinetic energy of relative motion with the internuclear separation cannot 
here be neglected) that cross-sections of order 10~4 7a)" are possible even if 
this time is as long as 10-!° sec. 

Another detachment mechanism remains to be described. The system 
formed on bringing two nuclei together may have no bound state for an 
excess electron. If the affinity becomes zero at an internuclear separation 
of Ra, the excess electron is free to leave the system at closer separations, 
and it is apparent therefore that the detachment cross section is of order 
(gh,")za)*._ The phenomenon is essentially the same as the ionization 
of atomic hydrogen by negative mesons which has recently been 
discussed by Wightman (1950). It should be noted, however, that the 
united limits which may be reached adiabatically may all be doubly 
excited ; and doubly excited states of a negative ion may be bounded even 
though the normal state is unbounded. Because of this, detachment 
through lack of affinity will not occur in many cases where it might at 
first be expected. However in such cases, there must of necessity be 
curve-crossing so that detachment may proceed through radiationless 
transitions as described above. 
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ABSTRACT 


Gallium has a marked anisotropy in its electrical resistivity ; these 
experiments, which cover the temperature range 0-2°K to 1°x, show 
that any dependence of superconductive transition point on crystallo- 
graphic direction is less than c. 0-002°K. 


§ 1 

GALLIUM is a superconductor with a transition point just above 1°x 
(De Haas and Voogd 1929, Shoenberg 1940, Goodman and Mendoza 
1951). Since it possesses the largest known anisotropy in electrical 
resistance, which has been shown to persist down to the superconducting 
transition point (Powell 1951, Olsen-Baer and Powell 1951), the question 
arises whether the superconductive properties are affected by this aniso- 
tropy. (Itis true that de Haas and Voogd (1931) could find no anisotropy 
in the superconductive properties of single crystals of tin, but the aniso- 
tropy in the resistivity is very much less marked than in gallium.) 

Two possible effects must be considered : 

(a) An anisotropic single crystal might at a certain temperature be 
superconducting in one direction and not in another. 

(6) The temperature at which the superconductive phase-change takes 
place in an external magnetic field might depend on the angle between 
the field and the crystal axes. 

Although it seems rather improbable that the first effect occurs, this 
is not so in the case of the second. It is true that both effects might be 
very small but in the absence of a complete theory of superconductivity 
we decided to look for them. 

Two different cryostats were used in both of which the specimen is 
immersed in a bath of liquid helium and is therefore in as strain-free a 
state as possible. One apparatus consisted simply of a bath of liquid 
helium which could be pumped down to just below 1°K ; the other was a 
demagnetization apparatus in which a bath of liquid helium could be 
cooled to 0:2°x (Croft et al. 1953). 


* Communicated by Professor F. E. Simon, C.B.E., F.R.S. 
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§ 2. MEASUREMENTS ABOVE 1°K 

Three gallium single-crystals grown along the directions of greatest, 
least and intermediate resistivity were prepared as already described 
(Olsen-Baer and Powell 1951) and were mounted in the liquid helium 
bath. Currents of up to 100 ma were used and the potentials measured 
with a galvanometer amplifier circuit as described in the paper just 
mentioned. The temperature was determined by measuring the vapour 
pressure and allowing for the Knudsen effect. 

It was found that in zero field and in a field of 10 gauss the transitions 
occurred at temperatures which were the same to the limit of accuracy 
of the experiment: that is, a few thousandths of a degree. 


§ 3. MEASUREMENTS BELOW 1°K 

Some different specimens of similarly high purity ‘were used in these 
experiments. Preliminary experiments showed that the transition 
points of the different crystals in various fields at temperatures down to 
0-2°K were very nearly the same and that they agreed fairly well with the 
data published by Goodman and Mendoza (1951). 

Two crystals with their long axes in the high and low resistivity 
directions were then mounted close together in the apparatus so that 
an accurate comparison could be made in a magnetic field which was at 
different angles to the crystal axes. Each specimen had an independent 
current supply and the potentials were joined in series opposition; in 
this way, if as the bath warmed up from 0-2°K one specimen became 
normally conducting before the other, a deflection in the galvanometer 
would be seen first in one direction and then in the other. By timing the 
interval between the swings, one could deduce the temperature interval 
from the rate of warming up. After demagnetization, a field of 55 gauss 
was applied by means of a solenoid and the galvanometer watched while 
the temperature passed through the transition point in this field, i.e. 
c. 0:24°K,. (The transitions were very sharp; this is a guarantee of the 
homogeneity of the field.) The field was then reduced and the same 
procedure repeated several times. Temperature measurements were 
taken magnetically at intervals. The time taken to warm up from 
0-2°K to 1°K was about 20 minutes. 

At all temperatures between 0-2°K and 1°K in several experiments 
consistent small differences were found of about 0-:005°K. To see whether 
these were due to inhomogeneity of field, temperature or the specimen, 
one of the crystals was cut in half, the two halves mounted exactly as 
before and the experiment repeated. We then found differences of 
about 0:003°K. Any real difference in transition temperature must 
therefore have been less than, say, 0-002°K. 


§ 4. CONCLUSION 
From these two sets of experiments, we conclude that any difference 


between the transition points of gallium single crystals measured along 
the different crystal axes is less than a few thousandths of a degree. 
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This work, which was carried out at the Clarendon Laboratory, also 
formed part of the general research programme of the National Physical 
Laboratory and is published with the approval of the Director. 
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ABSTRACT 


A theory is given of the statistical fluctuation of the number of 
neutrons in a pile. An experimental verification of this theory is 
described, and the results of the theory are used to evaluate the accuracy 
of the measurement of neutron absorption cross sections with a pile 
oscillator. 


§1. [yTRODUCTION 


For some years, a ‘ pile oscillator ’ has been used on the Graphite Low 
Energy Experimental Pile (G.L.E.E.P.) at Harwell for measuring thermal 
neutron absorption cross sections. A sample of the element to be 
measured is moved periodically in and out of the pile and produces a 
modulation in the number of neutrons in the pile. The depth of modulation 
produced by the sample is compared with that produced by a boron 
sample of known cross sectional area, and hence the total cross sectional 
area of the unknown sample may be determined. In this way the thermal 
neutron absorption cross sections of about forty of the elements‘ have 
been measured (Colmer and Littler 1950). 

The ultimate accuracy of the method depends on the fluctuations in 
the number of neutrons in the pile. This paper gives a theory of these 
fluctuations, describes an experimental verification of the theory and 


shows how the fluctuations affect the accuracy of the absorption cross 
section measurements. 


§2, DESCRIPTION OF THE G.L.E.E.P. Oscrnuaror 

A small graphite boat holds the sample to be measured and is 
oscillated between the outside and the centre of the pile. Transfer 
between these two positions is as rapid as possible (the twenty feet 
between the two positions are covered in one second), so as to minimize 
the effects on pile multiplication of neutrons scattered by the sample ; 
this effect vanishes when the sample is at the centre of the pile. At each 
of the two extreme positions the boat is stationary for an adjustable 
time, usually much longer than the transit time. Under these conditions 


* Now at University of Cambridge. 
+ Communicated by the Authors. 


> aes 


On Pile Modulation and Statistical Fluctuations in Piles 7 


the motion of the sample produces an almost perfect square wave 
modulation in the reproduction constant of the pile. The boat, which 
can carry loads up to about three kilogramms, is driven by a conventional 
type of servo mechanism, the only novel feature being that the position 
of the boat at the centre of the pile is always reproduced to within 
+ § inches. The servo mechanism is triggered twice during each period 
of oscillation by a master control unit which accurately governs the total 
period of the oscillator. 


§3. THEORY oF PLE MopULATION 


Since there are delayed neutrons in the pile, the square wave modulation 
in the pile reproduction constant produces a complex periodic modulation 
of the number of neutrons in the pile, superimposed upon the steady 
level NV, obtained when the reproduction constant has its mean value K 9. 
Suppose that the change in the reproduction constant produced by a 
sample going from the outside to the centre of the pile is 4. Let N, be 
the amplitude of the first harmonic component in the modulation of the 
number of neutrons in the pile, then, 


Dae On\ Ala). kl Vee 6 2-2 (3-1) 


A(w) is a function which is usually called the pile amplification 
factor. If A is small, it may be shown that A(w) is given by (Gilbert and 
Fergusson 1948) 


1/A?(w)=[wr—towp1,Aj/(Aj2+ wo?) 2+ [ew*Zy4,/(AP+o2)+1—KyP . (3.2) 


where w/2z is the frequency of oscillation of the sample, 

7 is the mean life time of a neutron in the pile (time being between being 
born as a fission neutron and being captured as a thermal neutron), 

K, is the reproduction constant of the pile when there is no modulation, 

cK, is the total number of delayed neutrons produced per thermal 
neutron captured in the pile, 

pu; is the fraction of ¢ corresponding to the production of delayed 
neutrons with a decay constant A;, and the summations extend over all 
the delayed neutron periods. 

K, is not quite equal to unity since some neutrons are produced in the 
pile by other means than neutron induced fission (e.g. by spontaneous 
fissions, alpha-neutron reactions, cosmic rays). However, in practice, 
provided that the pile is run at more than a few watts power we may 
put K,=1 in (3.2). 

Experimentally, Gilbert and Fergusson have determined the values 
of A(w) for periods ranging between 1 and 500 seconds by making 
measurements on the Zero Energy Experimental Pile (Z.E.E.P.) at Chalk 
River, Canada. Their experiment consisted in rotating a cadmium 
shutter at a fixed position in the pile in such a way as to produce a pure 
sinusoidal change in the pile reproduction constant, and then measuring 
the modulation produced in the number of neutrons in the pile. If 


K2 


128 O. R. Frisch and D. J. Littler on 


allowance is made for the different values of 7 in the two piles (G.L.E.E.P. 
is a graphite-moderated pile and Z.E.E.P. a heavy water moderated pile), 
their values may be used to give the values of A(w) for G.L.E.E.P. ; for 
modulation periods longer than 20 seconds, the values of A(w) for the 
two types of pile differ by less than 3%. 

Although the expression (3.2) is obtained from the pile kinetic equations 
by assuming that 4 is small, the Harwell Computing Group have shown 
by numerical integration that (3.2) gives the value of A(w) to within 3% 
for values of 4 up to 3 10-, and for all values of modulation period up 
to 100 seconds. 

For values of 4 up to 3x10-% therefore and for a fixed period of 
oscillation, we see from (3.1) and (3.2) that V,/N, is proportional to 4 and 
therefore to the cross sectional area of the sample being measured. 


§4. MEASURING ARRANGEMENT 


Experimentally, therefore, we measure the ratio N,/N, and take this 
to be proportional to the cross sectional area of the sample, the constant 
of proportionality being determined by a calibration with boron as 
mentioned in §1. The determination of the ratio N,/N,) may be 
understood by referring to fig. 1. The neutrons in the pile are detected 
by a boron trifluoride ionization chamber in series with a 108 ohm resistor. 


Fig. 1 
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Schematic diagram of electrical equipment. 
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MULTIPLIER 


To the earthy end of this resistor a potentiometer is connected, which 
‘ backs off’ the voltage produced by the current component corresponding 
to No. A pre-amplifier followed by a main amplifier then amplifies only 
the difference in voltage between the balancing potentiometer and that 
produced across the 10° ohm resistor. The voltage produced at the 
amplifier output is then of the same form as the modulation of the 
number of neutrons in the pile. A high speed recorder is connected to 
the amplifier output so that the pile modulation may be observed. 
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To obtain the desired Fourier component there is connected across the 
output of the amplifier a potentiometer with a sliding contact which 
moves with sinusoidal motion. This sliding contact is driven by the 
master control unit and has the same period as the sample oscillator. 
Because of the delayed neutrons in the pile, the first harmonic component 
in the pile modulation lags in phase behind the first harmonic component 
in the change in reproduction constant produced by the motion of the 
sample. By suitably advancing the phase of the motion of the sample 
with respect to that of the sliding contact on the potentiometer, it can 
be arranged that the first harmonic component in the voltage produced 
at the amplifier output is multiplied by a sine term in phase with it. 

After the sine multiplication the output voltage of the amplifier is 
fed into a velodyne integrator* and integrated over one modulation cycle. 
The sine multiplication followed by the integration thus performs an 
electro-mechanical Fourier analysis on the pile modulation. The answer 
recorded by the integrator is therefore proportional to N, multiplied by 
half the modulation period. If we denote the recorded integrator result 
by I(w), we will have from (3.1) 

I(w)/V =(2«a4/w)A(w), eee oo eee V4.1) 
where V is the value of the balancing voltage and « is a constant which 
converts the amplitude of sinusoidal voltages applied to the pre-amplifier 
input into answers recorded on the integrator. « may be measured by 
applying a known sinusoidal voltage to the amplifier input ; in fact this 
is frequently done during the course of modulation experiments, in order 
to ensure that the electronic equipment maintains a constant calibration. 

Thus for a constant modulation period the values of J(w)/V are 
proportional to the total cross sectional area of the sample being measured, 
provided that 4 does not exceed 3x 10-3. 

Statistical fluctuations in the number of neutrons in the pile will 
produce fluctuations in the recorded values of [(w)/V. Since we perform 
a Fourier analysis on the pile modulation in order to obtain the amplitude 
of its first harmonic component, it follows that the fluctuations in /(w)/V 
will only be due to the Fourier component of the pile fluctuations which 
has the same periodicity as the motion of the sample being measured. 
In other words, the method we use automatically eliminates a considerable 
proportion of the pile fluctuations. 

There still remains, however, the question of whether there is an 
optimum modulation period which can be used. That is to say, is there 
a modulation period such that the signal given by (4.1) divided by the 
‘noise’ due to the pile fluctuations has a maximum value? In order to 
answer this question we must know the Fourier spectrum of the pile 
fluctuations, and this is deduced in the next section. 


* Basically this is a system in which an electric motor rotates at an angular 
speed which is strictly proportional to the output voltage from the amplifier ; 
the angle through which the motor turns is therefore proportional to the 
integral of the output voltage. The system has been extensively described in 
engineering literature (e.g. J. Instn. Elect. Engrs., 1947, 94, Part IT, 112). 
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§5. Fourter Spectrum or STATISTICAL FLUCTUATIONS IN 4 PILE 


Since a pile is a chain reacting system, we must use that branch of 
probability theory which deals with stochastic processes. A compre- 
hensive account of the earlier work in this field is given by Bartlett, 
Kendall and Moyal (1949). ‘ 

For the sake of completeness we will state several theorems concerning 
stochastic processes ; the proofs of some of these theorems will be found 
in the papers by Bartlett et al. 

We may define a probability generating function (conveniently 
abbreviated to p.g.f.) by 

A=a,)+a,¢-+a.u74+ .... a,” ve oot eee ee 
where the coefficient of x” is the probability of obtaining an answer 7 in 
an experiment, and thus 2a,—1. If we denote (d4/dx),, by A, and 
(624A /ex),_, by A,,, then the mean value M and the variance v in the 
experiment are given by 

M=n=A,, \ 
v=(n—n=A,,—A,2+4,: 

If we consider two experiments A and B, and the p.g.f. for the sum — 
of their outcomes is C, then 


(5.2) 


C2-AB. 
Ce pe «ery, o> ee ee 
U,=V, +03: 


(It is assumed that the outcome of one experiment has no influence on 
the other.) 

As an extension of (5.3) we see that if we perform the same experi- 
ment K times (p.g.f. is A each time), then the p.g.f. is A™. 

We shall now consider what we may call a ‘ chain experiment’ in 
which the outcome of a first experiment decides how many times a 
second experiment is to be made. This type of combination will occur 
whenever we divide a chain reaction process into two stages. The first 
stage ends, for instance, with the production of m neutrons, and the 
probability distribution of the values of n is described by a p.g.f. 
P= p,x", In the second stage, every one of these neutrons will produce 
an independent reaction of its own, resulting in the occurrence of 
i fissions, @=2¢,a' being the appropriate p.g.f. To obtain the p.g.f. 
for the complete process, we remember that the p.g.f. for the added 
results of n identical experiments Q is Q”. Now if the first experiment 
always had the result », there would always be n neutrons at the beginning 
of the second stage, and the answer to our problem would be Q”. 
Actually, the first experiment allows a number of alternatives, each with 
its probability p,. Now the probability of any event is equal to the 
sum of the probabilities of all the alternatives which lead to it. Hence 
our required p.g.f. is 


R=pQ+p10!+p.Q?+ 5 oom & 
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This however, can be written more briefly as R=P(Q), i.e. letting P 
operate on Q. 


We may then show that 


eG) | 
Vp=P Qt Q,20p- 

In dealing with fluctuations in a reactor, the difficulty arises that what 
actually happens at any time depends on the number of neutrons present 
at that time, while what we are interested in is the number of fissions 
which have happened up to that time, or what we might call the fission 
score. As long as we merely want to know average values we can assume 
that the fission score grows at a rate proportional to the neutron number, 
but not when we are interested in the fluctuations. Instead we shall 
introduce a new type of p.g.f., a function of two variables « and y: 
P(x, y, t)=Z2f,,,(t)x"y!, where f,,(t) is the probability of finding, at the 
time t, a number n of neutrons in the reactor, and a fission score of . 

It is easily verified that the mean values and the fluctuations can again 
be found by forming the appropriate derivatives 

Fe =(0F 0x), y-1= 1, F,=(0F /0y).-y1= 1, 
FE, S(@ Fea), = n(n—1), Fy =(02F/0y?) . ya =Ul— 1), 
F j=(0?F /dxcy) 

The quantity nl is of no interest in itself but will appear as a link in 
the calculation. 

Consider a small time interval dé after a time ¢t, and what happens to 
one neutron in that interval.* 

If + is the mean neutron lifetime, the probability of the neutron 
number being zero at the end of dt is 


(dt/r)(1—K/v) 


where K is the reproduction constant of the pile, y is the mean number of 
neutrons per fissions, and the term in brackets is the probability of 
unproductive neutron capture. 

The chance that the neutron score is still one is 1—dt/r and the chance 
of fission taking place is (dt/7).K/v. Thus the p.g.f. for the chain 
experiment with one neutron in the time dt is (¢ if the p.g.f. for the 
number of neutrons per fission) 


(dt/r)(1—K/v) + (1—dt/r)a+ (dt/7)(K/v)o(y) . «+ (5.5) 


since if fission does take place we have a fission score of unity. In order 
to study the fluctuations of the Fourier component we consider a 
detector which has a chance y to detect a fission, and will contribute an 


(5.4) 


* We will neglect the presence of delayed neutrons in order to simplify the 
treatment of the problem. Their presence can be allowed for (unpublished 
calculations by O. R. Frisch) and does not alter the significance of the present 
results. 
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amount cos wt to the result for every fission detected. The p.g.f. for 
this detector will then be 1—y+yy°s and we may rewrite (5.5) as 
x-+(dt/7)[1—K|v—a+ (Ky) p(l—y ty )J=atg dt. . . (5.6) 

Now actually we start at time ¢ with n neutrons ; that is we can consider 
the pile plus measuring equipment as a chain experiment, where the 
exponent of x (the neutron number) decides how often the experiment 
described by the p.g.f. (5.6) is to be done in dt. Hence, we obtain the 
p.g.f. for the system at t+dt, that is P(t+dt, x, y), by letting F(t, x, y) 
operate on x+gdt; that is, we replace the x in F(t,2, y) by x+g dt 
(cf. the paragraph preceding eqns. (5.4)). 

Hence, F(t-+dt, x, y)=F(t, x+g dt, y). 

Also we must have an independent source emitting neutrons at a 
rate S (if dt is small enough for only one of these neutrons to be emitted 
during this time, the p.g.f. for this independent source will be 
1—S dt-+-S dt x), then, by (5.3), 

F(t+dt, x, y)\=F (t, x+g dt, y)(1—S dt+S dt, x). 
Using Taylor’s theorem, and discarding higher powers of dt, we get 
OF /ot=g OF /ox+S(x—1)F. ce 3 ee ee chee 
Differentiating this equation partially with respect to x gives 
(02F'/dx dt)=(0g/dx)(0F/0x)+-9(0?F /dx?) + S(x—1)(0F/0x)+SF 
and since ¢=1, g=0, F=1 when s=y=1 


Also from (5.6), 
Y= (K/v)6,— 1=K—1 


hence F,=F,(K—1)/t+8. 
If the pile is running steadily F 2=9, and 
F,=S7/(1—K). Te oy 


By differentiating (5.7) partially with respect to y, twice with respect 
to y, twice with respect to x, once each with respect to x and y, we may 
obtain 


Py=9,F 7, ee etry es keh) 
Fiyy=Gayl at 2G yF xy, SI ex ety Re Te Oey 
Pe= Qc tS) Fog ah ee ee ene Ta 
Bef igh pO el ny POEh ae SEO =) na) 
In which F’, is a constant given by (5.8) and 
Joa—=K$ p9]V75 Jay=(Ky/7) cos wt, 


\ (5.13) 


Jy=(Ky/vt) cos wt, — Jyy=(Ky/v7)(cos? wt—cos wt). 
Integrating (5.9) we have, since F,=0 when t=0, 
Py=(Ky/wv7)(sin wt)F,. >. . 2... (5.14) 


If the pile is running steadily, Pe as well as (oe must be zero. 
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Therefore from (5.11) we have 


Bp — Py Gan 28) 29... pe es r(5 510) 

Since at t=0, F, and v, are both zero, it follows that Fy, and ie are 
both zero at t=0. Also for t=0, Jyy=9;, and therefore from (5.10), 
i — Oat i=—0. 

Thus solving (5.12), using this boundary condition and eqns. (5.14) 
and (5.15), we have 
F =F ,(w?+d?) [A sin (wt—a)—Bsin (wt-+-B)— py/27?(w?— A?) 1/2] (5.16) 
where A=[y/2(1—K)][(p/7)+(2SK/»)], B=SKy/avt, 

tana=g,/w, tan B=w/q,, 
A=(1—K)/r, p=(K/v)?6,21+ 2K (1—K). 

We can now solve eqn. (5.10) for F,,,, using the boundary condition 
that F,,,=0 at t=0. In the sclution we insert the conditions that 
t=27m/q (i.e. an integral multiple of the modulation period) and that 
t>1/A. We then obtain finally 
oe EF Kyt 


py 
F 2Qvr {1 To 


yy 


since w/7mA <1. 
Again since F',=0 for t=27m/w, v,=F,,, so that, 
_ ytNy py 
v= mf 1+ ces} i oe a) 
where V,=F',K/vr, the mean fission rate in the pile. 

Referring to eqn. (3.2), we see that for the case we are considering, 
namely a model in which there are no delayed neutrons, C=0, and 
A®(w)=1/(1—K)*(1+*/A?). 

Thus the variance in the fission score is given by 

etry Nig Uni yA (CU) in to oat a) ee (5-18) 
We have seen that the value of p is given by p=¢,,,(K/v)?+2K(1—K), 
and as we have already mentioned, K is very nearly equal to unity for all 
pile powers above a few watts, so that writing v, for the variance in the 
number of neutrons per fission, 
p=(v,+v?—1)/v?. eee ge eee wa cL 9) 
Equations (5.18) and (5.19) therefore give the form of the Fourier 
spectrum of fluctuations in the pile. We have taken a model with no 
delayed neutrons, but will assume that (5.18) still expresses the form 
of the pile fluctuations correctly, provided that A(w) is now of the form 
given by (3.2). 

If we are only concerned with one modulation cycle of period 7 and 
we write v,=0,7(w) (o,(w) is the standard deviation in the fission score), 
then 

o,2(w)=TyN pl + pyA?(w) }/2. 


oc ert 
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Now if we are recording these fluctuations on the velodyne integrator, 
the integrator answers will be proportional to o,(#) with the same 
proportionality factor that the balancing voltage will bear to yN,—since 
this latter quantity is the mean fission rate recorded by the ionization 
chamber. Thus, calling o(w) the standard deviation recorded on the 
integrator, we have 

o(w)/V=a{T (1+ pyA2(w))/2yN pp}? . . . . (5.20) 
where « is a constant of the measuring equipment as explained in § 4. — 

When the pile fluctuations came to be measured, it was found that the ~ 
results were complicated by current noise fluctuations from the 10° ohm 
resistor. Current noise is observed only when current is flowing through 
the resistor and is said to be due to variation in contact resistance between 
the carbon composition granules which form the resistor (Bisby, Brown — 
and Brownrigg 1952). 

We conclude this section by showing how eqn. (5.18) is modified by — 
current noise fluctuations. | 

Suppose a pulse distribution with p.g.f. A be presented to a chamber 
which has a chance y to detect a pulse. The p.g.f. for the chamber will — 
be 1—y+ya, and from eqn. (5.4) we see that the mean value Z,, and the | 

| 
| 
: 


variance vz—o,” recorded by the chamber are given by (in the first 
instance we neglect the current noise) 


Z,=yA,, 
vz=Z,{1+y(v4—Az)/Az}s 
o7/Z4,={1+y(v4—A,)/A,}V7/Z,,2. Ss) eh Otanies 


Let us denote the p.g.f. for the current noise by J, the mean value of 
the current per pulse by M,, and the variance in the current by v;. 
Then if C denotes the p.g.f. for the experiment when current noise is 
considered, we have from (5.4) 

C,=MZ,=yA,My, 
0 ,=yA,v,+U/7Z,{1+y(v,—A,)/Ax}; 
o/C,={1+0,)/MP+y(v4—A,)/A,}VZ.4%. 2. . 2... (5.22) 


We see from (5.21) and (5.22) that the effect of current noise is to insert 
the quantity (1+-v,/M;*) instead of unity in to the expression for the 
‘noise to signal ratio’ for the experiment. 


: ns ; : ' : 
We may therefore take current noise fluctuations into account in 
eqn. (5.20) by writing 


o(w)/V=afT(1+Blw)+ pyA%w))/2yNp}¥2 .. (5.23) 
where the term B(w) represents the effect of the current noise fluctuations, 
and is a function of w since the current noise has a Fourier spectrum. 


For a given modulation period and a constant current through the 


resistor, B(w) will be constant ; it would, of course, be zero if the resistor 
were free from current noise. 


: 
: 
: 
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§6. VERIFICATION OF THE FoRMULA FOR THE PLE FLUCTUATIONS 


Before making measurements on the pile, it was thought necessary 

to test the electronic equipment on a fluctuation measurement that had 
been previously studied. To do this, the boron trifluoride ionization 
chamber was replaced by a high-pressure argon chamber close to' a 
200 millicurie radium source. 
Measurements were made on the integrator using a three-second 
Integration period, and the integrator answers were found to follow a 
Gaussian distribution. 

In this case, it can be shown that 


o(w)/V=a{T(1+0/m?)/2yr}1/? 


where y is the chance the argon chamber has to detect a gamma ray, 
m and v are the mean and variance of the pulses produced in the chamber, 
and r is the rate of emission of gammas by the source. 

Due to the difficulty in estimating y and v/m? exactly, it is not easy to 
check the width of the Gaussian against the foregoing formula, but the 
measured value of o(w) was certainly of the correct order of magnitude. 
We were therefore confident that the electronic equipment measured 
fluctuations correctly. Also, with the boron chamber re-connected but 
still out of the pile (no current from the chamber) it was found that 
the integrator recorded zero results. This showed that the pile measure- 
ments would not be complicated by any form of electronic noise other 
than current noise. 

Returning now to the measurements on the pile, it is, of course, not 
necessary to actually modulate the pile in order to measure o(w); a 
measurement with the pile running steadily at any given power level 
will yield the same values as a modulation experiment. Since it is much 
easier to make all the fluctuations measurements with the pile running 
at steady power, this was what was done. 

Our measurements were made with a boron neutron detector and not 
with a fission detector as assumed in the last section. To calculate y 
(eqn. (5.23)), therefore, we observe that it will be the ratio of the rate of 
detection of neutrons by the chamber, to the fission rate in the pile. 
Also the rate of detection of neutrons by the chamber is proportional 
to the chamber current and therefore to the balancing voltage, so that 
we may write y=CV/N >». 

It is easy to see that C is the reciprocal of the chamber resistor 
(108 ohms) multiplied by the mean charge produced per boron 
disintegration in the chamber. This latter quantity is only known with 
an accuracy of +30% for our type of chamber, so that y can only be 
calculated approximately. We have, therefore, avoided inserting 
numerical values for y. except to obtain an order of magnitude. 

We may now rewrite (5.23) as 


o?(w)/ VT =e {1+ B(w)4 pyA?(w)}/2C. ee ee OMe 
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Since we are trying to make an estimate of the standard deviation in 
the pile fluctuations, it is useful to know how many observations it is 
necessary to take in order to obtain an accurate measure of o(w). If the 
individual answers on a set of measurements taken on the pile follow 
a Gaussian distribution, then the standard deviation in o(w) is given 
by o(w)//(2n) (Whittaker and Robinson 1937), where n is the number 
of observations. 

To make sure that the individual answers recorded on the integrator 
do give a Gaussian distribution, one set of measurements was taken with 
a fixed chamber position and a three-second integration period. About 
2000 readings of o?(w)/V7' were obtained, and these are plotted in fig. 2 ; 
it will be seen that the distribution is Gaussian. 


No. OF OBSERVATIONS 


-36 -30 —24. -18 =(2 -6 ie) 6 12 18 24 30 36 
INTEGRATOR READING 


Gaussian distribution from pile fluctuations, 


Thereafter, only enough observations were taken for an estimate of 
a(w) to be made with about 10°% accuracy, the accuracy being computed 
from the formula-standard deviation is o(w)/+/(2n). 

It was also found that in many sets of measurements, the mean 
integrator answer obtained was not zero; this would be expected for a 
small number of observations. However, there was no significant 
difference whether o(w) was computed from the measured mean value 
or from a mean value of zero. 

The first experiments were made with the boron trifluoride chamber 
situated near the centre of the pile (y~10-4). It will be seen from 
eqn. (6.1) that for a constant modulation period (A(w#) and B(w) constant) 
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and a fixed chamber position (constant y), the values of o*(w)/VT will 
be independent of pile power. This was checked by a series of runs made 
at various pile powers with a constant integration period of 3 seconds, and 
found correct within the accuracy of the measurements. 

Next, the integration period was varied between 1 second and 
100 seconds. This range of periods is the range we are interested in 
using with the pile oscillator ; periods of less than 1 second cannot be 
used because of the time taken to move the sample in and out of the pile, 
and periods of greater than 100 seconds are unsuitable because it is 
found in practice that atmospheric pressure variations upset the 
measurements by changing the amount of nitrogen in the pores of the 
graphite and hence changing the reproduction constant of the pile. 


Fig. 3 


MODULATION PERIOD (SECONDS) 


AX) 
Pile fluctuations as a function of A?(w). 


In order to eliminate the term B(w) in eqn. (6.1), the measurements 
were repeated with a different chamber position such that y~10~°, the 
pile power being raised so that the chamber current, and therefore B(w), 
was the same as before. The term pyA?(w) in eqn. (6.1) is now negligible, 
and so we have 

{o2(w)/VT }o=a? {1+ B(a) }/2C. eerie eens ve sat Ose) 
Subtracting the individual sets of measurements in the two experiments 


cee PefoXlalfVI}4— {ow VT —atpVAI0) Ny. (6.3) 

The values of the fluctuation difference / are plotted in fig. 3 against 
the values of 42(w) determined by Gilbert and Fergusson. It will be 
seen that F varies linearly with A?(w) as predicted by (6.3). 
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From the slope of the line in fig. 3, and the known values of « and. 
N,/V, we may determine p. It was found that 


p=1-21-0-15 


the error in p being assessed somewhat arbitrarily since the errors in the 
experimentally determined values of A(w) are not accurately known. 
This fits the assumption that the distribution in the number of neutrons 
per 2U fission is Poisson, when v,=v, and therefore from (5.19) p=1. 
It therefore appears that the formula (5.23) correctly predicts the pile 
fluctuation spectrum and that p=1. 


§ 7. SENSITIVITY OF THE G.L.E.E.P. OscILLAToR 


We are now in a position to answer the question posed at the end of 
§4 as to what is the optimum modulation period to use with the oscillator. 
First, we notice that we should use a chamber resistor which does not 
exhibit any current noise, so that the term B(w) in (5.23) will be zero. 
This has now been done on the G.L.E.E.P. oscillator; the carbon 
composition resistor which was originally used has now been replaced 
by a 108 ohm wire-wound resistor, from which there is no observable 
current noise. 

Putting B(w)=0 and p=1, we see from (4.1) and (5.23) that the 
signal/noise ratio for a single modulation experiment of period 7’ is 


I(w) A (2yNypTA?(w)) V2 
= 1+yA*(w) } 

If we perform m modulation experiments each of period 7’, the 
signal/noise ratio for the mean of the m experiments, S(w), will be 
improved by a factor 4/m. So that writing 7)=m7' for the total time 
taken for these m experiments, we have for the signal/noise ratio of the 
mean, 


a(w) 7 


_ . Af 2yN pT A%(ew)) 172 
Srey rerun se 


For a fixed pile power (N, constant) a fixed modulation period 
(A(w) constant) and a fixed time 7’), S(w) will increase as y is increased 
and will tend to the limiting value 


S=A(2N pT 9)¥2/z. iy Se ee 


It will be seen that this expression is independent of modulation period 
and depends only on the total time taken for the set of experiments. 
We have the result, therefore, that provided that y is large enough, 
the signal/noise ratio is independent of the modulation period chosen. 

On the G.L.E.E.P. oscillator itself, we have made y as large as possible 
by placing the detecting chamber as near to the centre of the pile as 
possible. We have also replaced the boron trifluoride chamber by a 
boron coated chamber, since using this latter chamber we are able to 
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use high pile powers whilst still keeping the chamber saturated. The 
highest pile power we can use with the boron coated chamber in its 
position near the centre of the pile is 3kw, at which power 
NV p~10"* fissions/sec. 

Tn its position near the centre of the pile the value of y for the boron 
coated chamber ~610-® With this value of y we have the values 
(table 1) for the factor by which the theoretical maximum value for the 
signal/noise ratio (7.2) must be multiplied in order to obtain the actual 
signal/noise ratio for a given modulation period. 


Table 1 
Modulation period (seconds) 1 10 20 40 60 80 100 © 
yA) 12 ; \ 
4 Ey AX) 0:24 0:33 0-40 048 0:54 0:60 062 1 


We have chosen to use a 40-second period with the G.L.E.E.P. 
oscillator since this gives us about 0-5 of the theoretical maximum 
signal/noise ratio, and at this period our measurements are not disturbed 
very much by atmospheric pressure fluctuations, which at 100 seconds 
modulation period would have been much more troublesome. 

A normal run on the G.L.E.E.P. oscillator consists of twenty measure- 
ments each of 40 seconds period (7’)=800 seconds) done at a pile power 
of 3kw. We would therefore expect the signal/noise ratio for the mean 
of the twenty measurements to be 


S=0:48A(2 x 10" x 800)¥2/x 
=6x 107A. Rt eg AP e578) 


It is of interest to know what is the standard error in measuring an 
absorption cross section in such arun. From the boron calibration of the 
oscillator it is known that 1 cm? of cross sectional area gives a value of 
A~1-4x10-*. (Since, as mentioned in §3, we can use values of 4 up to 
3x 10-3, this means that we can, in principle, measure cross sectional 
areas up to 200cm?.) From (7.3) we have, therefore, S=8-5 x 1072, 
where » is the cross sectional area being measured (expressed in cm?). 
The standard error in determining a cross section is the value of 2’ given 
by this expression when S=1,; the value found is 0-12 mm’. 

It is found in practice that the standard error of twenty measurements 
is of the order of +0-16mm?. This is remarkably close to the value of 
+0-12 mm? which would be expected if the pile fluctuations were the 
only cause of variations in the integrator answers. 
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SUMMARY 


The conductivity and Hall constant of grey tin alloys containing anti- 
mony and galiium were measured over a temperature range of 77°K 
to 286°K. From the results the number of charge carriers and their 
mobility were calculated as functions of temperature. Large variations 
of conductivity and Hall constant in a magnetic field were noted. 
A comparison is drawn with the similar semiconductors silicon and 
germanium. 

In the intrinsic range the conductivity is given by 

o=1-6X 104 exp [—0-064/2K7). 
The energy gap between the filled and unfilled bands, 0-064 ev, does not 


vary significantly with temperature. The number of charge carriers is 
therefore given by 


NN ,= 2:4 108! T® exp [—0-064/KT] 
combined with the expression for electrical neutrality 
N,—N},=Ng 
where n, is the number of charge carriers due to the impurity atoms. 
In the p-type alloys the Hall constant reverses sign, and from the 


temperature of this reversal the ratio of the lattice scattering mobilities 
for electrons and holes can be calculated as 


C=p,[Ma= 1-3. 
In the lattice scattering range the mobilities are given by 
[p= 11-6 x 108 T-3?, 
fin= 8:9 108 T-3?, 


§1. INTRODUCTION 


ELECTRONIC conduction in grey tin has been studied by Busch, Wieland 
and Zoller (1951), Blum and Goryunova (1950), and Kendall (1950). 
It has been shown to be a semiconductor similar to silicon and germanium 
with a separation between the filled and unfilled bands of approximately 


* Now at The Plessey Company Ltd., Caswell, Towcester. 
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0-lev. In the present investigation measurements of electrical conduc- 
tivity and Hall constant have been made over a temperature range from 
77°K to 286°K on grey tin specimens containing known amounts of antimony 
and gallium. One alloy containing lead was also investigated. The data 
are analysed to determine the variation of the concentration of current 
carriers and their mobility with impurity concentration and temperature. 
The results are interpreted in terms of existing theories of conduction in 
semiconductors. 

The results presented here differ considerably from the above earlier 
results, particularly in the lower value obtained for the intrinsic activation 
energy (0-064 ev) and in the much more coherent picture obtained of the 
variation of the number of charge carriers and their mobility with 
temperature. 


§2. PREPARATION OF SPECIMENS 


‘Spectroscopically pure’ tin was obtained from Johnson-Matthey 
Ltd. Chemical analysis by the suppliers gave the following results :— 


Lead 0-0012% Antimony 0-001% 
Tron 0-00027% Copper 0-0002°% 
Arsenic 0-0002% Bismuth 0-0001% 
Sulphur 0-0003% Tin 99-997, 


The alloys shown in table 1 were made by melting the constituents in 
vacuo, and were then converted to the grey modification by seeding with a 
minute quantity of pure grey tin and storing at —50°c in vacuo until 
conversion was complete, as determined by an x-ray powder photograph. 
During the conversion the metal crumbles to a fine powder. 


Table 1. Compositions of Tin Alloys 
1 2 3 4 5 6 7 8 
Ata Carri ¢| Sign) Carriers 
Sample|Impurity| 9%, "©! Atoms/ce ea ei “| of | due to Carriers/ 
added | 2dded temp. car-| added atoms 
: riers| impurity 
(1) — 0 0 3°25 x 1017] — | 0 — 
(2) Sb 0:0143 | 4-18 x 1018 | 2-52 x 1018| — | 2-20x10!8| 0-53 
(3) Sb 0-0279 | 8-17 4-08 — | 3-76 0-46 
(4) Sb 0-0740 | 2-16 x 1019 | 8-47 — | 8-15 0-38 
(5) Sb 0-310 | 9-07 2-16 10!9| — | 2:13x10!9| 0-24 
(6) Ga 0-0060 | 1-76 x 1018| 7-0 1017) + |1:02x1018| 0-58 
(7) Ga 0-0119 | 3-48 1-24 1018} + | 1-56 0-45 
(8) Ga 0-0806 | 2-36 x 1019] 8-15 + | 8-48 0-36 
(9) Ga 0-322 | 9:42 2-70 X 1078} + | 2-7310!9| 0-29 
(10) Ga 0-063 | 1-84 101° 
Sb 0-056 | 1-63 x 101° 
Difference 2-1 x1018| 4-40 x 1017} + | 7-65x1017| 0-36 
(11) Pb 0-025 | 7-32 x 1018] 8-00 x 10!7| — | 4-75 1017] 0-065 
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Unsuccessful attempts were made to prepare crystals by electrolytic 
deposition at a temperature below the transition point (13-2°c), evapora- 
tion in vacuo on a cold surface of grey tin, crystallization from 
concentrated hydrochloric acid, and crystallization from mercury. 

It was found possible to press grey tin powder into cylindrical specimens 
(0-25” diameter and slightly over 0-5” long) by carrying it out at a tempera- 
ture of approximately —125°c. A pressure of 10 tons/sq. in. was used 
and the mould was cooled with liquid nitrogen. The specific gravity of the 
specimens was 5-75, and no lines of tetragonal white tin were observed in an 
X-ray powder photograph. Usually some reversion to white tin occurred 
on the side-surfaces of the specimens. This was easily removed either by 
dissolving the surface in dilute hydrochloric acid or by storing the specimen 
at —50°c overnight. 

If grey tin powder at ordinary temperature is subjected to pressure 
it reverts back to white tin. A study of the dependence of the transition 
temperature on pressure has been made by Komar (1936) and he finds that 
it changes according to 07'/0P=—5 x 10-2 °c/Atmosphere. A pressure 
of 10 tons/sq. in. will thus lower the transition temperature to —62°c. 
Using a mould cooled with liquid nitrogen a thermocouple measurement 
of the temperature of the inside of the mould filled with grey tin powder 
was —125°c under working conditions. 


§3. EXPERIMENTAL PROCEDURE 
3.1. Conductivity 


The electrical conductivity was measured by a potential-probe method. 
The voltage drop along the length of the specimen was measured indepen- 
dently of the current carrying circuit by the introduction of probe 
electrodes. The specimen holder is shown in fig. 1, and it is seen that the 
polystyrene tube enclosing the specimen is perforated by two holes 
through which thin platinum-iridium wires (0-06 mm diameter) are pushed 
into contact with the grey tin. These holes were placed 0-50" apart and 
arranged to be near the end surfaces of the grey tin cylinder. 

The specimen holder was immersed in a vacuum flask filled with methyl 
cyclohexane for temperatures down to about —80°c and with liquid 
propane for lower temperatures. This liquid medium was kept stirred 
with a small motor-stirrer and was cooled by additions of liquid nitrogen. 
For the very lowest temperature the specimen holder was immersed directly 
in liquid nitrogen. 

The temperature was measured by means of an Fe/NiCr thermocouple 
recessed into the polystyrene cylinder as shown in fig. 1. The thermo- 
couple was calibrated at several fixed temperatures. After addition of 
liquid nitrogen to the bath an interval of about 5 minutes was required for 
the attainment of temperature equilibrium. In practice about 10 minutes 
was allowed. It was found convenient to carry out measurements at a 
series of fixed temperatures corresponding to 103/7’'=13, 12, 11, 10, 9, 8, 7, 
6-5, 6, 5-5, 5, 4:5, 4, 3-5. 

L2 
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The electrical measuring circuit was conventional, a d.c. potentiometer 
being used for measuring all voltages. In practice the use of the probe 
electrodes did not appear to be really necessary as the voltage drop across 
the current-carrying contacts was always very small compared with the 
voltage drop along the length of the specimen. 


3.2. Hall Effect 


Platinum-iridium probe electrodes were placed on opposite sides of 


the grey tin specimen as shown in fig. 1. The specimen holder and its — 
container were placed between the poles of an electromagnet, the Hall — 


probes being transverse to the direction of the field. The current through 
the sample and the direction of the magnetic field were both reversed 


Fig. 1 
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independently and the four Hall voltages thus obtained generally differed 
slightly (presumably due to thermo-electric effects). The mean value 
was taken. The magnetic field strength was determined by measuring 
the current through the windings after calibrating by means of a small 
search coil and ballistic galvanometer. 

The Hall constant R, in em?/Coulomb, was calculated from the relation 


R=10VUHI eee 


where V is the measured Hall voltage in volts, ¢ is the thickness of the 
sample in cm, H the magnetic field strength in gauss, and J the current in 
amperes. 

The use of a cylindrical specimen for Hall measurements is unusual, 
but it is not thought necessary to introduce any geometric factor into the 
above formula. Isenberg (1948) has shown that for a rectangular 
specimen with a ratio of length to width of 2 (as used here) the Hall 
constant is reduced to 0-93 of its normal value. However, for a 
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cylindrical specimen this reduction will be smaller, and we have therefore 
taken the geometric factor to be unity in the present case. 


ee Variation of R and o with Magnetic Field Strength 
considerable variation of & and o with H was noticed with the purer 
a especially at low temperatures. This effect was used to estimate 
mol ilities of charge carriers. The current through the specimen was 
ee constant, and values of the potential difference along the 
ength of the specimen and of the Hall e.m.f. were measured as functions 
of field strength. 
§ 4. EXPERIMENTAL RESULTS 
4.1. Conductivity (at zero magnetic field) 
The results are given in fig. 2. For the straight-line portion of curve (1) 
(pure grey tin) the conductivity is given by 
Ue-GGOxXD|——e/2K 1 | oe. ts ca cs 4 (2) 
where o)7=1:6X< 104 ohm-! cm-! and «0-064 ev. 
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This value of ¢ is also obtained from the other curves, except the 
alloy containing approximately equal amounts of gallium and antimony. 
Slightly different values of o, are obtained for the different alloys, the 
smallest value being for the purest material. 


4.2, Hall Constant (at zero magnetic field) 
The results are given in fig. 3. It will be seen that as the temperature 
is decreased R increases exponentially until a saturation value is reached 


corresponding to the number of impurity charge-carriers. 
In the case of alloys (6)-(10) positive values of were obtained, 
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indicating hole-conduction, except at high temperatures in the case of 
alloys (6), (7) and (10). In these curves a reversal of sign took place 
and from the curves a temperature can be found for which R=0. 


4.3. Variation of R and o with H 


The results for alloys (1) and (6) are given in figs. 4 and 5. The effect 
was also measured for the other alloys. In the case of alloy (1) the 
variation of R is plotted against H? for low fields in figs. 6 and 7. This 
procedure is necessary in order to extrapolate to H=0, thus finding the © 
Hall constant at zero field. The same procedure was followed for all the 
alloys. 

Fig. 3 
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Variation of Hall constant with temperature. 


§ 5. INTERPRETATION OF CoNDUCTIVITY AND Hat CoNsTANT 
MEASUREMENTS 
(a) It will be noted that « does not vary with temperature (within the 
limits of experimental error). 
We may therefore use the relation derived from statistical mechanics 
(Fowler) : 
NN,= 2-4 X 1081 x T3 xX exp [—e/KT]. . . (3) 


This relationship will hold whatever the purity of the sample. 
In the case of Si a multiplying factor of 32-5 is needed on the right-hand 
side of this equation, and in the case of Ge a multiplying factor of 3:5. 
If this is taken to be due to a linear variation of ¢ with 7, 


i.e, €p=€9—B pr. e - e e e e e e (4) 
Then we have 


exp [—«7/KT]=exp [8/K]xexp[—e«,/KT]. . . . (5) 


The multiplying factor is thus exp [8/K’]. 
Since in the case of grey tin 8 must be very small we may take the 
multiplying factor to be unity. 
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(b) From fig. 3 (Hall constant for tin-gallium alloys) the maximum 


values of FR are: 
Alloy (6)~10-6, Alloy (7) 6-0. 


Assuming these alloys are non-degenerate at low temperatures and that 
n,<n, then we may calculate the number of charge carriers due to 
impurities as 

ne= —7-4X 1018/R=(6) 7-0 1027, (7) 1:24 1028. 


Also from fig. 3, R=0O when 103/7'=(6) 5-65, (7) 4:90. 


Fig. 6 
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Variation of Hall constant with magnetic field, alloy (1). 
Using eqn. (3) we may calculate n,n, at the temperature at which R=0, 
.n,=(6) 2:08 x 1088, (7) 5-55 x 1086, 


Also, at all temperatures 
Ne=N),—N,. ey ee a Be (i). 


From eqn. (3) for n,n, and (6) for n,—n, we may now calculate 
n, and n,, separately at all temperatures. 


—_—_ 
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For temperatures at which R=0 we obtain : 
(6) n,=1-05 1018, »,—1-75< 1028, 
(7) 2,=1-82 1018, m,=3-05 x 108, | 
Now at R=0, the Hall constant is given by 
R=—T-4xX 1018(n,C?—n,)/(n,C-+n,2?=0 . . . . (7) 


where C=p,/u,. 


Therefore, at these temperatures, C= \/(n,/n,)=(6) 1-29, (7) 1:30. We will 
take C=1-3. 


R (em Coulomb.) 
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Variation of Hall constant with magnetic field, alloy (1). 


(c) The critical temperature (7’,) for the transition from non-degenerate 
to degenerate conditions is that for which the energy of an electron or 
hole corresponding to the energy at the surface of the Fermi distribution 
is equal to KT,. 


We thus have 
T= (h?/8 km)(3/7)??9n?89=4-2 x 10-11y2/8, het es) 
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(d) Using eqn. (6) and the relationship 
R=—T7-4 x 1018(n,C?—n,)/(m,C-+m,)?> «. . « + (9) 
we may now calculate n, and n,, separately at all temperatures for the 
non-degenerate samples. 

For the degenerate samples we may substitute the numerical factor 
6-25 in place of 7-4 in eqn. (9). For samples which are on the border-line 
between degenerate and non-degenerate we may use a weighted mean 
of the two numerical factors. 

In the lattice scattering range we take C=1-3. In the impurity | 
scattering range we take C=1. In the transitional range a weighted 
mean (obtained graphically) must be used. In practice this variation — 
of C makes only a slight difference in the cases of some of the alloys — 
over a range where n, = 7, corresponds with the transitional region. 


Fig. 8 
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Variation of number charge-carriers with temperature. 


The values obtained for n, and n, are not shown separately, but 
values of n,+m, are given in fig. 8. The boundary between the 
degenerate and non-degenerate samples is shown by plotting 7’... It is 
seen that the number of charge carriers increases. exponentially in the 
intrinsic range and is constant in the impurity range. 

(e) Using the relationships 


OPEN Le HEN gion 5) ie 0 | ee LO 
and Be=1:3Xp, . . . . . . . . (11) 


we may now calculate values for 4, and j«,, in the lattice-scattering range. 
In the impurity-scattering range we must take 1,—.,. For intermediate 
ranges we must take a weighted mean of the numerical ratios 1-3 and 
1-0 (found graphically). In practice this makes only a slight difference 
and then only when n,n, in the intermediate range. The results 
are given in fig. 9. 
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It is seen that in the lattice-scattering range , and ., are given by 
the equations 


Pgh RO Er 1k ena Ee): 
Wepre sO C1 00 nS ee epee et at (19) 


Impurity scattering gives rise to lower mobilities and the greater the 
number of impurity atoms the smaller the mobility. 


§6. Discussion 
6.1. Ratio of Charge Carriers to Impurity Atoms 
In column 5 of table 1 the number of charge carriers per c.c. has been 
calculated from the value of the Hall constant at low temperatures. 
A positive value indicates that the carriers are holes, a negative value 
that they are electrons. 
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Variation of carrier mobility with temperature. 


The pure sample, to which no impurity has been added, contained 
3:251017 free electrons at low temperature. This represents the 
carriers due to impurities originally present in the ‘ spectroscopically 
pure’ sample. The value of 3-25 x 1017 is therefore subtracted from the 
n-type samples in column 5, and added to the p-type samples, to give 
the values in column 7 which are the carriers due to the added impurities. 

Tt will be noted that the values in columns 7 and 4 are not the same, 
and their ratio is given in column 8. If every impurity atom added 
gave rise to one carrier this ratio would be unity. Tt is less than unity 
because a proportion of the added impurity atoms is ineffective in giving 
rise to carriers. Presumably this is due to a proportion of the added 
impurity atoms being precipitated at crystal boundaries or elsewhere 
and thus not entering the grey tin lattice. It will be noted that the 
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“= 


more impure the sample the smaller the proportion of effective impurity 
atoms. In this, antimony and gallium appear to behave similarly. 

In the case of lead, either only a very small proportion enters the 
lattice substitutionally or some other factor is operating to cut down 
the effectiveness. It is difficult to understand why lead should give rise 
to a carrier at all, as it has the same valency as tin. It does not seem 
likely that the carriers are due to an impurity in the lead as ‘ spectro- 
scopically pure ’ material was used. 


6.2. Concentrations of Electrons and Holes 


The values of n, and n, as functions of temperature may also be 
calculated without knowing the Hall constant by making use of the 
theoretical equation (3) for the product n,n,. The values so calculated 
come out to be within the limits of experimental error in the measurement 
of the Hall constant throughout the entire range of temperature. 

In the case of silicon, Pearson and Bardeen (1949) did not obtain such 
close agreement and therefore cast doubt on the adequacy of the theory 
of the Hall effect in the intrinsic range. In the case of grey tin the theory 
seems to be quite adequate. 


6.3. Mobilities of Electrons and Holes 


For the non-degenerate samples the mobilities of electrons and holes 
approach the same values at high temperatures for all impurity concen- 
trations. This is due to lattice scattering and it is found that w varies 
asl 32: 

At lower temperatures impurity scattering becomes important, 
although the temperature range investigated did not extend sufficiently 
low for an accurate evaluation to be made of this effect. Making a rough 
comparison with the impurity scattering observed in silicon alloys by 
Pearson and Bardeen (1949) it is remarkable that in grey tin the scattering 
by impurity atoms seems to be less effective by a factor of about 100. 
The reason for this is not at all apparent. 

For the degenerate samples the mobility should be relatively insensitive 
to temperature, and this is seen to be generally true, although alloy (3) 
which should be degenerate at low temperatures shows a significant 
decrease in mobility at low temperatures, i.e. it seems to behave as though 
it were non-degenerate. 


6.4. Tin—Gallium—Antimony Alloy 


This alloy has a comparatively large impurity content, but only a small 
number of charge carriers owing to the mutual neutralization of n- and 
p-type impurity centres. The p-type impurities are in slight excess. 
The conductivity-temperature relation is unusual in that the slope of 
log o/7'" is not constant in the intrinsic range (fig. 2). 

. The Hall constant has the value one would expect from the difference 
in the Sb and Ga contents. The mobility approximates to the value one 
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would expect from the sum of the Sb and Ga contents. It would thus 
appear that neutral scattering atoms are just as effective as ionized ones. 


6.5. Tin—Lead Alloy 

Although the impurity concentration is comparatively large the 
number of free charge carriers produced is small and the Hall constant is 
therefore large. The mobility corresponds to approximately what one 
would expect from the added impurity concentration so that a much 
larger proportion of the lead must go into substitutional positions than is 
effective as a carrier producer. 

The simplest assumption to make would be that lead itself contained 
an impurity (e.g. antimony) which was the effective electron-producing 
agent, while the lead remained only a neutral scattering agent. However 

about 10% of impurity in the lead would be required, and in fact spectro- 
- scopically pure lead was used. 


6.6. Variation of Conductivity and Hall Constant in a Magnetic Field 
(a) Conductivity 
At low fields we have 


(092) /oy= 3 x (47) x2? i aah pects.» Le) 


(Pearson and Suhl 1951, Estermann and Foner 1950). 

Values of (cy—c)/o, for alloys (1) and (6) are plotted against H in 
figs. 4 and 5. It is seen that at low fields the H? law is obeyed but at 
higher fields the relationship falls to a lower power of H. From the 
values in the low-field region an estimate of the mobility of the charge 
carriers may be obtained. Introducing numerical values into eqn. (14) 
and changing from e.s.u. to practical units, we obtain 

(o9—0)/op=0-38u7H? x 10-18. ee eee (ko) 

Values of mobility calculated from this equation are given in table 2. 
It is seen that much higher values of mobility are obtained by this 
method, especially at low temperatures ; also that the mobility of holes 
appears to be higher than that of electrons. These effects have also 
been observed by Pearson and Suhl (1951) in germanium. 


(6) Hall Constant 
Electronic theory gives a similar formula for the variation of the Hall 
constant with transverse magnetic field 


9772 
(Ry—R)/Ro= Gre, eh we es fa) 
or in practical units 
(B,—R)/Ro=1:38u2H?x 10-8, 2. (17) 


Values of (Ry>—R)/R, for alloys (1) and (6) are plotted against H in 
figs. 4 and 5. It is seen that at low fields an H? law is obeyed, but at 
higher fields the relationship falls to a lower power of # (even lower 
than in the case of the magneto-resistance effect). From the values in 
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the low-field region an estimate of the mobility of the charge-carriers 
may be obtained, and these are given in table 2. It is seen that the 
values in general agree with those estimated from the magneto-resistance 
effect. 


Table 2 
Mobilities calculated from :— 
(R,—R)/Ry)H?=1-38uR? x 10-18 
. (¢9—2)/o,H? =0-38u,? x 10-16 
. fe (from § 5 (e)) 
| py (from §5 (e)) 
: Hm 


Ou Go bo 


(mean mobility) =(%,u, +7, p)/(M +p) 


Alloy (1) Alloy (6) 


10°/T | wR bo Me Pn BR bo Ln Lm 

13 35000 | 35700 | 16000 | 16000 | 37000 | 41800 | 13000 | 13000 
12 33400 | 33800 | 14800 | 14800 | 35000 | 39600 | 11700 | 11700 
11 30600 | 32300 | 13300 | 13200 | 31000 | 36200 | 10400 | 10400 
10 28100 | 28400 | 11900 | 11600 | 27500 | 32000 | 9000 | 9100 


9 25400 | 25000 | 10000 | 9550 | 24000 | 27600 | 7500 | 7650 
8 18000 | 20400 | 8400 | 7800 | 18500 | 21700 | 6400 | 6700 
ul 16300 | 16200 | 6800 | 6200 | 12000 | 15900 | 5200 | 5600 
6-5 10400 | 13100 | 6300 | 5600 12600 | 4600 | 5100 
6 6300 | 11200 | 5400 | 4700 10600 | 3800 | 4400 
5:5 5000 | 7200 | 4750 | 4200 200 | 3600 | 4000 
5 4100 | 5600 | 4000 | 3500 5600 | 2950 | 3400 
4-5 3200 | 4800 | 3600 | 38000 | 2600 | 4800 | 2750 | 3000 
4 2100 | 3600 | 2900 | 2600 3400 | 2250 | 2550 
35 2800 | 2400 | 2150 2800 | 1900 | 2100 


6.7. Activation Energy of Impurity Electrons 


The activation energies for the ionization of donators and acceptors 
can be estimated as approximately 0:005ev. This follows from the 
Bohr formula 


EL = 27? me*/K*h?. 
Setting « (dielectric constant) equal to 50 (Moss 1952) gives 
E,=13-5/502=0-005 ev. 
This corresponds to a temperature of 60° abs. so that we may assume 


that all the impurity centres are ionized throughout the temperature 
range investigated. 


6.8. Comparison with other Homopolar Semiconductors 


Some of the fundamental electronic properties of the homopolar semi- 
conductors are given in table 38. 

The variation of the energy gap ¢« with the crystal lattice constant a 
can be empirically correlated by means of the relationship 


e=6-2 exp [—a>*/6400]. 
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This relationship is shown plotted in fig. 10 and it is seen that all the 
points lie exactly on a straight line. It is interesting that the substance 
silicon carbide, which is nearly homopolar in its properties, also lies very 
near the straight line. Moss (1952) gives a similar empirical relationship, 
but the one given here is more accurate. 


Fig. 10 
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Relationship between activation energy and lattice constant, 


The mobility of holes in grey tin is the same as in germanium in the 
lattice scattering range. The mobility of electrons is somewhat smaller. 
In the impurity scattering range the mobilities in grey tin are very much 
larger than in germanium or silicon. 

In table 4 the intrinsic conductivities of the homopolar elements are 
compared at various temperatures. There is seen to be an enormous 
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variation from diamond to grey tin. In table 5 the number of free 
intrinsic charge carriers is similarly compared. 


Table 4. Intrinsic Conductivity 
(Calculated from o=o, exp [—«/2K7']) 


ohm—! em-1 


T(x) C Si . Ge Sn 
1000 1-2x10-1 14 480 11000 
500 1:3 x 10-24 23x 10 7 7500 
290 2-4 10-48 2-1 10-6 1-4 x 10-2 4500 
200 9-1 «107-11 2-0x 10-5 2500 
100 9-1 x 10-25 1-3 10-14 400 
Table 5. Number of Charge Carriers (calculated from n=2(n,n,)!/?2= 


2(2°3 X 103173 exp [—e¢/KT'] exp [8/K])!”, 
for C, 32-5 for Si, 3-5 for Ge, 1-0 for Sn) 


where exp [8/K]=7:4 


T (°K) C Si Ge Sn 

2000 7-6x 10" 

1000 8-710? 2-8 x 1018 8-0x1018 | 2-2 102 
500 3-3 x 10-8 15x 1015 4-0x1018 | 5-1x1029 
290 3 x10-25 | 6-0x101 4-010 | 1-3x 1029 
200 1-6 x 10° 1-410 | 4-3x 1018 
100 5-5 x 10-9 5-6 1-2 x 1017 

20 8-6 x 10° 
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XIX. A Method for the Measurement of Multiple Scattering 
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Fysisk Institutt, Norges tekniske hogskole, Trondheim, Norway* 


[Received November 5, 1953] 


SUMMARY 


A new method for the determination of the multiple scattering of 
short, highly scattered discontinuous tracks is described. 


§ 1. INTRODUCTION 


For the scattering measurement in nuclear emulsions there are many 
methods in vogue. In most cases the Fowler coordinate method with its 
modifications is the best one. However, for the measuring of short, highly 
scattered, discontinuous tracks, e.g. from electrons near the minimum 
ionization, this method is not suitable. A method will be demonstrated 
here which, in this particular case, may be of some use. 


§ 2. THEORY 


Imagine an infinitely long track the projection of which in the 
x-y-plane is represented as a Fourier integral 


hee 
ae} | o g(w) exp (twx) dw. 


Then g(w) contains the statistical properties of the track. To this track 
we lay a chord of length 2A (fig. 1). t(s,A,é) is the distance from the 
chord to the track, measured at a distance s from the middle of the chord, 
which lies at ~=€. Since the angle between the chord and the x-axis is 
small, we can write 


A+s A—s 


t(s, A, €)=y(€+s)— aa YlE+A)— DA y(€—A) 
+ 00 vie 
= ib ; g(w) exp (tw€) [exp (iws)— se exp (twA)— A exp (—iwd) | dw. 


By squaring and integrating with respect to €, using 


ato) Peo OM? 


ee ee ee eee 
* Communicated by Dr. H. Wergeland. 
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(Olsen, Wergeland and Overas 1953), we obtain the mean square distance 


<t7(s, A) )= se =| ott Ce aee (A—s) w 


i igeco® 12 P2212 A 
Ss a cos(h-+5)0-+ ieee cos 20} 
or 
(P(s, A) )= aH (3) ) = «#0, d)) (1 (3) ) ean) 


<P(s, A)) is thus the value to be expected from measurements of the 
projection of many tracks or track pieces produced by electrons with 
equal energy. 

($7(2X) )/2A is the characteristic quantity of the scattering process, viz. 
the mean square of the projected scattering angle per unit length. It is 
related to the mean square distance at the middle of the chord <t2(0, A)) by 


2(2d 
— 2 — © (0, d)). Pao og to (2) 
Fig. 1 


The relations (1) and (2) can also be found from the diffusion equation or 
by a direct method developed by H. Wergeland (unpublished). 

To find <#2(0, A)) experimentally one may proceed in many ways, for 
instance straight forward by measuring 4(0, A) for the tracks available. 
In order to get more information out of the tracks, we should rather 
measure all the grains in each track and use (1) to find <#(0, A)). The 

game information can, however, more easily be obtained in the following 
way : . | 
We imagine all the projected tracks under consideration placed in such a 
way that their chords of length 2A coincide. Then the grains belonging 
to all these tracks will have a certain density p(s, t, A) at a point (s, t), and 


we may write 
J + co 
{ p(s, t, A) dt 
OSA) = oon eee ae 
| p(s, t, A) dt 


M2 


(3) 
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For the reason of continuity of the paths the integral 
+0 
{ p(s, t, A) dt=po 


should be independent of s. | 
By integrating the two expressions for ¢#?(s, A) ), (1) and (3), with respect 


to s, we get | 


[ "ee. A))ds= ¢#2(0, ))20 [1-5 +3 =| ==|" #2 P(t, A, o) dt. 


: (4) 
+o 
P(t, A, 0) dia at | o(s, t, A) ds 
is the number of grains between —o and +o (o<A) in a strip dt parallel : 
to the z-axis. We do not include the regions near the ends of the chord, © 


because of relatively large measuring errors there. 2cp,) is the total 
number of grains between —o and +o. 


Be bra Bie 
ib #2P(t, A, o) dt, : 
20 po 
which we shall denote by <7?(A,c)), is the mean square of ¢, taken — 


together for all the grains between s=—o and s=-o, irrespective of their 
s-coordinates. 
Combining (2) and (4 


(4) 
KP*(2A))_ {TA @)) 8 
2 


we then obtain 


Mi 2erabes] a 
etait 


The relation (5) tells us what we should expect to find after measuring 
an infinite number of tracks, provided there were no measuring errors. © 
In practice the large grain size introduces some uncertainty, so that — 
instead of ¢7(A, c)) we find a value (7'(A, o))* from the measurements. — 
It may be shown with reasonable assumptions that : 


(D(A, o) = (T(A, 2) )*— (A?%(a/A)), 
where ¢4*(c/A)) is of the order of magnitude (D?)/8, slightly varying 
with o/A. D is the grain diameter. The corrected formula is 
(P2(2A)) (TA, a) )>*— (A?(a/X)) 6 ' 
Dy ae ; SP wre ye Eats Las 8 . . . . . (6) 
| os al 


From (6) we see that the correction term in <6?(2A))/2A varies as A-3, 
and we therefore get the best determination with fairly long chords. 

(A*(c/A)) may be estimated using two different chord lengths rx 
and Ag, with 03/Ay=0,/A,=o/A 


(4%(o/)) = Ee ee 


= (une? eo ee a | 
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§ 3. EXPERIMENTAL VERIFICATION 

The relation (6) was examined experimentally by means of 40 tracks 
produced by 0-59 mev electrons. 

($7(2A))/2A was determined from (6) for values of o/A ranging from 
1/15 to 12/15 (table 1). Apart from statistical fluctuations the values of 
($?(2A))/2A should be independent of o/A. The largest value of o/A gives 
statistically the best determination, but how the standard deviation of 
($7(2A))/2A decreases with increasing o/\ is not easily predicted, since 
measurements on the same track here are statistically dependent. If we 
measure NV grains in each of M tracks, it can, however, be safely stated 
that the standard deviation lies between what would be expected for 
and N x M independent observations. 


Table 1 Table 2 


(a/A) x 15 


1 
2 
3 
4 
5 
6 
ti 
8 
9 
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i 
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DAAADRAAABDBDINDD 
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The constancy of pp) was also checked (table 2). In these measurements 
the correction ¢4?(c/A)) was of the order of 5% of <7(A, o))*. 


§ 4, THz Meruop IN PRACTICE 


For practical measurement a micrometer eyepiece as indicated in fig. 2 
would be convenient. The line a is fixed and has two marks at a distance 
2\. The line b has marks at +c and —o and can be displaced in the 
direction of the arrows by a micrometer screw. 

Before measuring, the stage of the microscope is moved until the track 
crosses the line aat +Aand —A. Then the micrometer screw is turned and 
reading made whenever the line b passes the centre of a grain between —-o 
and —oc. Such readings are taken for all the tracks to be measured, 
and the mean square of all these readings gives (7?(A, o) )*. 

The method is limited to chord lengths less than the diameter of the 
microscope field, but the stage needs not to be moved, except for putting a 
new track into position. We therefore get no stage noise. 


162 On a Method for the Measurement of Multiple Scattering 


As an example the method was used to determine the multiple scattering — 
of 0-59 Mev electrons in nuclear emulsions, and the result obtained by 


Hisdal (1952) was confirmed. 


Fig. 2 
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SUMMARY 


Using the electrical resistance as an indicator for ordering it is found 
that Cu,Au specimens prepared in a state of long-range order and then 
water quenched from 380°c (8°c below the critical temperature) show an 
ordering effect at about 100°c. After this has decayed a similar effect is 
again observed after an electron bombardment (energies of 0-5 and 1 Mev) 
or after a small deformation at room temperature. An activation energy 
H~0-9 ev for all three cases is determined. Attributing the ordering 
effect to the migration of vacancies or interstitials the experiments 
provide evidence that one or both of these defects are obtained in 
thermodynamic excess by quenching or deforming. 


§1. INTRODUCTION 


In a previous letter (Adam, Green and Dugdale 1952) the effect of a 
1 mev electron bombardment on a highly-ordered Cu,;Au specimen was 
reported. Using the electrical resistance at 0°c as an indicator it was 
shown that a further ordering at 100°c took place after the bombardment. 
This was attributed to the migration, at 100°c, of vacancies or interstitials, 
produced by the bombardment, which allowed some of the remaining 
wrong atoms to move to right positions. The experiment also indicated 
that the ‘ active’ defects were gradually used up. 

Apart from the technique of bombardment one may hope to introduce 
a thermodynamic excess of vacancies into a Cu;Au specimen by two 
other methods. The first of these is by quenching from a suitably high 
temperature. If U is the energy required to form a vacancy the 
equilibrium concentration at a temporature is approximately 
exp (—U/kT) (see Mott and Gurney 1948). Thus, if the vacancy 
concentration is allowed to approach its equilibrium value at 7’, a rapid 
cool to a lower temperature should produce the desired eXcess. | The 
second is by a deformation at a suitably low temperature. As indicated 
by Mott (1951) and Seitz (1952) the movement of dislocations occurring 
during deformation may be expected to generate vacancies. Thus, a 


* Communicated by the Authors. 
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specimen deformed at a temperature low enough to prevent a rapid 
approach to equilibrium should contain the desired excess.* 

The aim of the present experiments was to find out whether, in fact, 
a quench or deformation would produce an effect similar to that of 
electron bombardment, i.e. a low-temperature ordering at c. 100°C. 
In order that the conditions for comparing the three effects should be as 
nearly as possible similar all specimens were first quenched and then 
used for electron bombardment or deformation studies. 


§2. APPARATUS 


The alloy used was made by Johnson, Matthey & Co., Ltd., in the form 
of wire of diameter 0-25 mm from specimens of Cu and Au of high purity. 
An analysis made by the Chemistry Division at A.E.R.E. gave 
25-8-+0-2 atomic °%% Au with Ag, Fe and Si as major impurities with a 
concentration of 0-01 to 0-1%. 

All the specimens were initially constructed as shown in fig. 1 (a). 


Fig. 1 
=A~2cm— mica 
pare ales a 
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spot welds Al. clips 
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They were heated in an argon atmosphere at about 900°c for two hours 
after which coring could not be detected by x-ray diffraction methods. 
Specimens intended for electron bombardment were then bent into an 
M shape and beaten flat to a thickness of 0-025—-0-:05 mm as shown in 
fig. 1 (b). This shape could be circumscribed by a circle of 8 mm diameter. 
All specimens were annealed in vacuo at 600°c for half an hour to remove 
the effects of cold work (which, in the case of the straight specimens, 
occurred to some extent on their removal from the homogenizing furnace). 
The straight specimens (denoted by ‘A’) were mounted in a mica 
sandwich as shown in fig. 1 (c). The electron bombardment specimens 
(denoted by ‘ R’) were attached to bakelite mounts after heat treatment. 
Great care was taken to avoid bending these specimens. 

The critical temperature 7’,, upon cooling at 1-5°c per hour, was 
measured on one straight specimen against a Pt—10°% Pt Rh thermo- 
couple. It was found to be 388°c according to a standard calibration 
table for the couple. 

To obtain a high degree of order specimens were maintained at 380°c 
for three days and then cooled to room temperature at 20°C per hour. 


* It is possible that interstitials may also be produced by these methods. 
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For quenching, a vertical furnace with a temperature regulation of 
--l°c was used. The temperature was measured by means of a 
Pt—13% Pt Rh couple permanently attached. Specimens to be 
_ quenched were placed in a small Al box and this, together with a steel 
plunger and some Mg ribbon to act as a getter, were sealed inside a 
glass tube evacuated to 10-2mm Hg. The tube was attached to a wire 
and drawn up into the furnace from which it could be dropped into a 
water bath by simply cutting the wire. The presence of the steel plunger 
ensured an immediate breaking of the glass tube on impact with the water. 
The Al box, the interior of which was easily accessible to the 
water, prevented the specimens from being damaged. It was necessary 
to maintain ordered specimens for two hours at 391-392°c in this furnace 
before quenching to remove their super lattice lines. This is in rough 
conformity with the measured value of 7', for slow cooling, mentioned 
above, and served as a rough check on the furnace temperature 
measurement. 

The resistance at 0°c of the specimens was measured with a precision 
better than +1 part in 20 000. They were not suitable for a determination 
of their resistivity p. However, a suitable specimen was made and after 
quenching from 450°c gave, for the disordered resistivity (0°c), 
p=11-6u2cm-1. The values of p for the other specimens, in a particular 
state of order, could then be determined from their disordered resistance. 

Isothermal curves of resistance & (0°c) v. time ¢ at a temperature 7’ 
were obtained by immersing the specimens in thermostated oil baths for 
suitable periods of time. The oil bath temperatures were measured to 
0°1°c. The temperature regulation was +0-1°c up to 100°C and 
+0-2°c above 100°c. 

A Van der Graaf machine was used for electron bombardments. The 
specimen was placed (in air) in a uniform beam defined by an 8 mm 
diameter hole in a brass plate. The temperature of the specimen (which 
was cooled only by convection in the air) was continuously monitored by 
measuring its resistance. It was not allowed to rise above 40°c. 
A measurement of the electron flux at the specimen was made on several 
occasions by replacing it by a suitably designed constant flow water 


calorimeter. wie 
The displacement of atoms by electron bombardment is discussed in 


the Appendix. 


§ 3. EXPERIMENTAL 
(i) The Effect of a Quench Followed by an Electron Bombardment 
An initially highly ordered specimen Rl (p=5-4y2cm™') was 
maintained at 380°c for two hours then quenched (p=9-5 pu cm~?), 
On heating at 130°c its resistance decreased as shown in fig. 2 (a). At 
the point B it was bombarded with 10” electrons/cm? of 1 Mev energy. 
During the bombardment the resistance decreased. about 0-5% as 
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indicated by D. Subsequent heating at 100°c and 130°c produced the 
fall in resistance shown. 

Specimen R25 was maintained at 380°C for four Cae then quenched 
(p=7:542cm-1). Subsequent heating at 100°c and 130 c produced the 
fall in resistance shown in fig. 2 (b). At the point B it was bombarded 
with 5106 electrons/cm? of energy 0-5 Mev. Subsequent heating at 
100°c caused the resistance to fall as shown. 


Fig. 2 
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Specimen R2 was maintained at 450°c for three hours then quenched 
(p=11-642cm~'). The small decreases in resistance shown in fig. 2 (c) 
occurred on subsequent heating at 100°c and 130°c. A bombardment 
with 101 electrons/cm? of energy 1 Mey at the point B caused the further 
small decreases on heating at 100°C and 130°c shown. 


(ii) The Effect of a Quench Followed by a Plastic Deformation 
An initially highly ordered specimen A15 (p=5-42cm-!) was 
maintained at 381°c for three hours then quenched (p=7:9 1.2 cm), 
Heating at 100°c afterwards caused R to fall as shown in fig. 3 (a). After 
500 min at 100°c the specimen was removed from its mica sandwich 
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and given a small deformation at room temperature. The deformation 
consisted of bending the central portion of the specimen into a single 
loop (~8 mm diameter) and then straightening it again. The whole 
operation, including remounting and returning to the ice bath, took 
about 5 min. When R was next measured (about 30 min later) an 
increase of nearly 0-5°% was found. Subsequent heating at 100°c for 
a further 500 min produced the fall in R shown. 


Fig. 3 
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An initially highly ordered specimen A3 (p=5-4 2 cm™) was ye 
tained at 450°c for three hours then water quenched (p=11-6 p2 cm). 
Heating at 100°c for 500 min produced hardly any change in R as shown 
in fig. 3(b). The specimen was then deformed in the same manner 
as Al5 with the resulting increase in R of about 0-2%. Further heating 
at 100°c had very little effect as shown. 


(iii) The Determination of an Activation Energy E 
The preceding experiments indicate that relatively large effects at 
c. 100°c may be induced in specimens in a state of long-range order by 
a quench, a bombardment or a deformation. Specimens without long- 
range order (those quenched from 450°c) show relatively small effects. 
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In this section, more detailed experiments on specimens in a state of 
long-range order are described. The object was to determine activation 
energies for the three types of effect on the assumption that 


dR 

TTS) exp (HET). 6 5) eee 
where R is the resistance at 0°c. To obtain # it is necessary to measure 
dR/dt at two temperatures 7’, and 7’, for the same value of R, when, 


from (1) Ae ve 
_ kT ,T, ak dk ] ‘ 
Po 77,8 a) a) 
rhe hie Ie 


For this purpose, sets of isotherms of R (0°c) v. t, at several tempera- 
tures in the range 70-130°c, were obtained on two specimens, A2 and R25, 
after quenching from approximately 380°c. A set of activation energies 
was determined for each using the values of the gradients at the ends 
of the curves (measured with a ruler) in (2). Similarly, further sets of # 
were obtained from A2 after a small deformation at room temperature 
(similar to that, described above, for Al5 and A3) and from R25 after 
a 1 Mev electron bombardment.* 

The isotherms obtained (except those from R25 after quenching) are 
shown in fig. 4.+ There was a considerable scatter in the values of # as 
exemplified by the detailed values for A2 shown in table 1 (a). 
Table 1 (b) shows the mean values of # for the two specimens. ‘These 
all agree within the R.M.S. deviation and a figure of approximately 
0-9ev may be taken as the value of the activation energy for the 
low-temperature ordering process induced by the three treatments. 

Some additional results on quenched specimens are now mentioned. 
Isotherms were obtained from three specimens, A6, All and A16, in the 
quenched ordered condition, under identical conditions and at the same 
time as A2. A6 had received an initial heat treatment with A2 (four days 
at 379°c then quenched) whereas All and A16, after being initially 
highly ordered, were maintained together at 379°c for two hours then 
quenched. The values of H, initial p and total percentage decrease in R 


*The detailed treatment of these two specimens was as_ follows. 
A2: Maintained at 379°c for four days then quenched (p=7-3 p.2 cm-}). 
After obtaining the quench isotherms, was heated at 130°o (curve 9, 
fig. 4 (a)) until dR/dt was very small, then deformed (D, fig. 4 (a)). R25: 
Originally maintained at 380°c for four days, quenched and used in various 
preliminary bombardment experiments. Later, maintained at 381°c for 
two hours, quenched (p=7-7 12 cm-!) and heated at 130°c until dR/dt was 
small. Bombarded with 81016 electrons/em? at 1 Mev (during which 
k decreased 0:3%) then used to obtain the bombardment isotherms. 
Maintained at 381°c for three hours, quenched (p=8-3 «4 cm-1) and used 
for the quench isotherms. 


Tit is interesting to note that the rate of ordering at 0°c was sufficiently 
high for an isotherm to be obtained in a reasonable time (~6000 min) after 
deformation and after bombardment. This was not the case after quenching. 
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during the isothermal treatments obtained from these specimens are 
compared with those of A2 in table 2. The values of EZ agree, within the 
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The following key to the curves gives the temperature in °c and total time in 
minutes at temperature. 


(a) (0) 
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scatter, with those already mentioned and show that # is not very 
sensitive to the type of heat treatment before quenching. 


Table 1 (a). Values of H from A2 
(i) After Quenching 


°o) | 68-8 89-4 68-8 89-4 107-9 128-3 110-8 
°o) 89-4 68-8 89-4 107-9 128-3 110-8 128-4 
| 0-77 0-91 0-79 0-95 0-92 0-88 0-94 


(ii) After Deforming 


°a) 0 68:8 O 68:8 89-4 68:8 110-8 1284 110-8 


( : 
(2a) 68:8 OO 68:8 89-4 688 110-8 1284 110-8 1285 
(ev) 0:86 0-91 0-90 0-87 0-96 1216 eS) 1-14 1-04 


Table 1 (b). Mean Values of # from A2 and R25 


Specimen Mean £ (ev) | R.M.S. deviation | No. of values 
A2 after quenching 0-88 +0-07 a 
2 after deforming 0-97 +0-12 9 
R25 after quenching 0-88 +0-14 13 
R25 after bombarding 0-84 +0-08 18 


Specimen | Initial p (u2cm-!) | %AR E (ev) R.M.S. deviation 


A2 7:3 1-81 O:SSaamai| +0-07 
A6 7:2 1-60 0-88 +0-04 
All 8-0 2°32 0-94. +0-06 
Al6 8-1 2-22 0-92 +0-10 


These results also show a variation in °/,4R between specimens given 
an identical preparation which may be due to insufficient precision in the 
. . > 
quenching technique. 


§4, Discussion 


The experiments demonstrate that an appreciable resistance decay at 
c. 100°C, somewhat similar to that following an electron bombardment, 
does occur after quenching or deforming specimens in a state of long-range 
order. A plausible explanation is, as in the bombardment case, that 
these decays are ordering effects. The agreement between the activation 
energies is sufficiently good to suggest that the ordering mechanism is the 
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same in all three cases and therefore that a thermodynamic excess of 
vacancies (and/or interstitials) is indeed produced by a quench or a 
deformation. If the estimate made in the appendix is not too far out it 
would appear that these ordering effects may be brought about by a 
vacancy concentration of less than 1 in 10° or 106. The magnitude in 
the fall in resistance suggests therefore that each vacancy must, during 
its migration, put many wrong atoms right before it is captured. 

If this explanation is correct the physical significance of H(~0-9 ev) 
possibly lies in some rate controlling process in which the vacancies are 
involved. The scatter in the values of EH is attributed mainly to errors in 
f and dR/dt but these may obscure a possible fine dependence on, for 
instance, the type of treatment causing the effect and the degree of 
order of the specimen. 

The small effect in specimens without long-range order (quenched 
from 450°C) appears to be consistent with the above explanation since 
one would expect their resistance to be relatively insensitive to any 
ordering effects which may have occurred. Indeed, Sykes and Jones 
(1936) showed, by observations on the specific heat, that ordering does 
occur at c. 100°C in specimens quenched from above the critical 
temperature without affecting the resistance appreciably. 

The increase in resistance upon deforming is possibly due to three 
things : (i) a decrease in order (Dahl (1936) showed that an ordered Cu,;Au 
specimen could be completely disordered by sufficient deformation) ; 
(ii) an increase in the number of dislocations or similar large faults ; 
(iii) a change in the physical dimensions of the specimen. The two 
latter would seem to be the cause of an increase in Rk when a disordered 
specimen (e.g. A3) is deformed. 
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APPENDIX 
The Displacement of Atoms by Electron Bombardment 


The maximum recoil energy W,, which an atom can receive when its 
nucleus elastically scatters an electron of energy His, very approximately, 


2H (H+ 2m,c?) (3) 


Wi Mc 
where m,c2 and Mc? are the rest energies of electron and atom, 
respectively. 
The number of atoms displaced in a given bombardment can be 
estimated if the minimum displacement energy W, and displacement 
cross section op» are known. W, can, in principle, be measured, op, the 
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cross section for the recoil energy to lie between W, and W,,, may be 
estimated by suitably integrating the differential scattering formula of 
Mott (1932). Because of mathematical difficulties the integration of an 
approximate formula is necessary. Since, from (3), W,, for Au is about 
one-third that for Cu it is probable that, at low electron energies, displaced 
Cu atoms will predominate in Cu,Au. Thus, only op for Cu need be 
determined. For this purpose the following expression is used. 


1—f? Wm 1 (B2-4-na) log = 
0 


B* LW 
+2nap(/ GP ) [pars ere eh iS 


where Z is the nuclear charge, .=Z/137 and B=v/c. This formula is 
derived from the « approximation of McKinley and Feshback (1948). 

W, has not yet been measured. However Seitz (1949), has suggested 
that it will be of the order of 25ev. That it is certainly less than 26 ev 
is demonstrated by the fact that ordering can be induced by a 
0-5 Mev bombardment for which W,, is 26ev for Cu. The fraction 
of atoms f displaced during a 1 Mev bombardment of integrated flux 
¢=101" electrons/em? may be estimated if we assume that W)—25 ev. 
From (3) W,,=68 ev for Cu which gives in (4) ¢p=30 barns. Thus 

J=60,=3 10% 

Dixon, Meechan and Brinkman (1953) bombarded a well-ordered Cu,Au 
specimen at a temperature below —185°c with 1 Mev electrons of 
integrated flux 3-4 10! electrons/em?. According to the above calcula- 
tion this should have displaced a fraction of approximately 10~° of the 
Cu atoms. An increase in resistance expected to be about 5° (on the 
basis of calculations made by Dexter (1952)) due to these displaced atoms 
was looked for but not found. These authors concluded therefore that 
many of the interstitial vacancy pairs formed must have been unstable 
at the temperature of the bombardment and so recombined. Thus, it 
is possible that the number of interstitial vacancies pairs formed during a 
room temperature bombardment which are capable of dissociating to 
form free vacancies and interstitials may be very much less than the 
above calculation indicates. 


op=0:2522 
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ABSTRACT 


In the production of radioactive isotopes by bombardments with 4N 
ions accelerated in the Birmingham sixty-inch cyclotron, total initial 
activities of some hundreds of thousands of counts per minute have been 
produced in aluminium targets. The following radioactive nuclei, in 
order of yield, have been identified with certainty among the products : 


ied Ee 28 A], 501i BAS sop 3S; iS, 24N a. 


In addition to these isotopes, 27Mg, #3N and “C have been detected with 
certainty ; the evidence for 8F, ®8Cl and 2?Na is inconclusive. 

Bombardment of aluminium with 1°O ions was shown to give 38K, #4Cl, 
iP and Na. . 

Most of the assignments of activity have been confirmed chemically. 
The experimental observations are not easily explained if the colliding 
nuclei simply fuse into a compound nucleus which then undergoes 
evaporation or fission; it is suggested that usually only part of one 
nucleus sticks to the other. 


§ 1. INTRODUCTION 


ALTHOUGH in the last two decades an enormous amount of work has been 
done on the nuclear reactions produced by artificially accelerated particles, 
almost all of this has been concerned with bombarding particles no heavier 
than «-particles. This emphasis is natural in view of the complexity of 
the reactions which might be expected in bombardments by heavier 

“nuclei, but it is also a reflection of the difficulty encountered in producing 
useful beams of the heavier nuclei. Considerable progress has now been 
made, however, in the acceleration of ions of elements such as carbon, 
nitrogen and oxygen, so that studies of nuclear reactions produced by such 
ions incident on target nuclei ranging from deuterium to uranium have 
become possible (Miller et al. 1951, Ghiorso et al. 1951, Wyly and Zucker 
1953, Reynolds et al. 1953, Chackett et al. 1953). In a theoretical paper, 

Breit ef al. (1952) have also pointed out that reactions between nuclei in 
the region of nitrogen might be expected to give useful information on 


nuclear structure. 
* Communicated by the Authors. 
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Some of the radioactive products from the bombardment of aluminium 
with 14N ions have already been noted (Chackett ef a/. 1953). The present 
paper is concerned mainly with a detailed examination of these products 
and possible explanations of the relative yields observed. 


§ 2. ACCELERATION OF MULTIPLY-CHARGED [ons 


The chief problem in accelerating heavy ions in a cyclotron is that of 
obtaining ions in a high state of ionization. The method we have used for 
accelerating ions carrying 6 electronic charges (such as *N** and 1®O®*) in 
the Birmingham 60 in. fixed-frequency cyclotron involves a preliminary 
‘third harmonic ’ acceleration of the ions in the 2-charged state (Walker 
and Fremlin 1953). Through collisions with gas atoms in the cyclotron 
tank the 2-charged ions become stripped to 6-charged ions. A number of 
experiments, described elsewhere, have been carried out on this mechanism 
(Walker et al. 1954). We give here only the data essential to an under- 
standing of the main topic of the present paper. 

All bombardments of aluminium have been made inside the cyclotron at 
a radius of 25in. Here the circulating beam contains up to 101° 4N&+ 
ions per sec having energies over 50 Mev ; the number of ions per unit 
energy interval is roughly halved for each rise of 10 Mev above 50 Mev, 
while the maximum possible energy is about 125 Mev. Similarly for 
1606+ jons, beam intensities of 3 x 108 ions per sec with energies over 40 Mev 
have been obtained, the maximum energy being about 140 Mev. The 
strong decrease of intensity with increasing energy means that the ions 
with energies not very far above the nuclear potential barrier will produce 
the greatest number of reactions. This will hold good at all depths in the 
target. 


§3. BoMBARDMENT OF ALUMINIUM witH “4N§&+ Ions 


The choice of aluminium as a target was due mainly to simple practical 
considerations. Aluminium is available in thin strong foils of high 
purity ; it consists of but a single isotope and the chemical difficulty of 
separating the possible products is not great. The target usually consisted 
of a number of foils each 12-5 « thick, one above the other and all covered 
by a protective aluminium foil about 2 thick. This was needed because, 
besides the desired beam of six-charged ions, relatively large beams of 
slower doubly-charged particles exist, moving with energies approaching 
the potential barriers of elements such as carbon, which are found in the 
residual gas and are liable to be deposited on the target. The cover foil 
would reduce to negligible proportions in the main part of the target any 
products arising from the low-energy part of the beam. In fact, but for 
the single product Na, control counts of the protective foils show no 
significant qualitative differences from those of the main foils, which shows 
also that no serious activity arises due to slow non-resonant contaminating 
ions. 

The first few foils of the stack usually formed the ‘ window ’, about 
2cemxlcm, in a shielded probe (see fig. 1), any further foils being 
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mounted on the inner current-collector. The outer tube of the probe 
carried two windows, on opposite sides, so that two successive bombard- 
ments could be carried out without withdrawal, by turning the probe 
through 180° ; this was often useful when short well-defined bombardments 
were needed, as the preliminary adjustment of the ion beam could be done 
using one window and the main bombardment using the other. No 
significant radiation hazard exists, of course, with such small beams, so 
that the probe could be rotated by hand while the beam was present. The 
probe could be removed through an air lock and the foils detached and 
placed beneath a counter within 14 minutes so that reaction products with 
half-lives down to a little less than 1 minute might be observed. 

The radioactive products of the bombardments were identified by 
analysis of gross decay curves obtained after bombardments lasting for 
periods ranging from two minutes to many hours, and also by radio- 
chemical analysis of the foils. In some cases, e.g. ?*Mg and !4C, the 
evidence was entirely from chemically separated specimens. 
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Target used in the bombardment of aluminium. 


3.1. Identification of Products by Analysis of Gross Decay Curves 

In fig. 2 is shown a typical curve for a foil bombarded for 2 min., observed 
with a Geiger counter with an aluminium end-window 7 mg/cm? thick. It 
will be seen that several decay periods are involved. Curve 2 has been 
analysed in figs. 3 (a) and 3 (6) and it will be seen that, besides longer 
periods, it breaks down perfectly, within the statistical errors, into a set of 
curves of half periods 150 min, 33 min, 7} min and 2-3 min. 

Taken in order of decreasing half-life, the products were provisionally 
identified as 32P, 24Na, 348i or !8F, *8Cl or *4Cl, 8K, and 78Al or 2°P. While 
the lack of alternatives made several of these identifications seem fairly 
reliable, three pairs of alternative possibilities existed which could not 
easily be separated simply by analysis of decay curves. These were 
319i and 18F, 38Cl and °4Cl, 9°P and Al. The first pair could easily be 
distinguished chemically, the second pair could not be so distinguished 
and the third pair involved some difficulty owing to their short half-life. 


N2 
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An elementary positron-electron discriminator was therefore made up, — 
using a permanent magnet with a flux density of about 1800 gauss in the | 
14 in. gap between its poles. This gave a discrimination of about 7 to lL 
between f-particles of opposite sign. For the sign of charge which it was 
desired to accept, the counting rate was reduced to about 1/5 of that 
obtained with the source as close as possible to the counter window in the 
absence of the magnet. 

Positive and negative decay curves could be obtained simultaneously 
from this device. As an example, analysis of curves obtained from a 
source bombarded for 2 minutes gave the following information. Firstly, 
both positive and negative curves showed a weak activity of 2-3 hours 
half-life, suggesting that both 18F and *1Si were present. Secondly, not 


Fig. 2 
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Decay of activity in aluminium bombarded for 2 minutes by ™“N ions. 
only the 7; minute activity but also the 33 minute activity was 100° 
positive (probable error under 5°) showing that the main, if not the Rie 
chlorine isotope was *4Cl. Thirdly, the 24 minute activity was 90-45% 
negative, thus confirming that the main 2} minute activity was 28A] hue 
suggesting the possible presence of a trace of a short-lived positron emitter 
presumably either 9°P or 4O. ; 
This charge-discriminator was also used to check the sign of the 
particles emitted by each of the chemically separated fractions mentioned 

below except for the short-lived aluminium and phosphorus fractions. 
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The possible occurrence of ?2Na calls for special comment. The first 
curve analysed showed an extraordinarily large quantity, the yield in 
number of atoms formed being 10 times that of any other isotope. This 
_ curve, however, was derived from the first long bombardment, in which no 

precautions were taken to protect the foil from surface contamination. 
The next bombardment, performed with a cover foil as already described, 
showed a yield less than 1/500th as great. By counting also the cover 
foil itself, it was found that almost all of the 22Na appears in the top 2 u 
of target. This may be due to any of several causes irrelevant to the main 
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Analysis of the later part of fig. 2. 
Curve A shows points obtained after subtraction of the activities due to **P 
and #4Na from the gross decay curve. The straight part represents 


318i with a trace of F’. fs 
Curve B shows the points obtained after subtraction of the activity due to 
319] 4.18F from curve A. The straight part (about 6 half-lives) represents 


220]. 

investigation. The two most likely are (1) that ?2Na is sputtered from the 
dees, which are almost certainly contaminated, as some 40 mC of this 
isotope have been made in the cyclotron in the last two years, or (2) that it 
is produced by the interaction of the very large current of 14 Mev 14N?* ions 
with carbon from the pump oils. A serious discoloration due to the latter 
is always observed on the cover-foil after long runs. 

Absorption measurements have been made of the radiations from a 
four-month-old target foil which had been properly protected from such 
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sources of contamination. This, on analysis, showed some 10 counts/min. 
from #P, about 20 counts/min from *S and less than 2 counts/min of **Na. 
If indeed this last isotope is produced at all, therefore, by the bombardment 
of aluminium by the high energy ™N¢+ ions, the yield must be less than 
0-003 of the #2P yield. 


3.2. Chemical Identifications 


Standard radiochemical methods were used in the identification of the 


radioactive products, the only unusual feature of the investigation being 


Fig. 3 (5) 
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Further analysis of fig. 2. 

Curve C shows the result of subtracting activities due to **P, *4Na, 21Si-+-18F, 
and “Cl from the original gross decay curve. The straight part (about 6 
half-lives) represents °K. 

Curve D shows the result of subtracting the activity due to °8K from Curve C. 
The straight part (about 8 half-lives) represents 8Al with a trace of 2°P, 


the devising of successful methods for the separation of the various 
activities from such a wide variety of radioactive products in this part of 
the Periodic System. Fortunately aluminium is dissolved very rapidly 
by both hot caustic soda solution and by hydrochloric acid. Radio- 
chemical analysis of the solutions so obtained, after adding suitable carriers, 
as well as of the gases evolved during solution, clearly indicated the presence 
of radioactive isotopes of the following elements (in order of increasing Z). 
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(1) Carbon (11C): counted as K,"CO,. This activity is present in the 
gases evolved when a bombarded foil is dissolved in hot sodium hydroxide 
solution. The 1C is presumably carried as a hydrocarbon which, after 
combustion with excess of oxygen, is converted to 14CO, and may now be 
trapped on powdered potassium hydroxide. 

(2) Nitrogen (183N): counted as (18NH,) ,SO,. This activity also 
appears in the gas phase on dissolving bombarded foils in boiling sodium 
hydroxide solution. It may be trapped and counted on a filter pad 
moistened with dilute sulphuric acid. 


Counts per minute 


minutes 
Decay of activity in Zr (HPO,), (after subtraction of *?P) showing *°P. 


(3) Fluorine (}8F) : precipitated and counted as Cal, from neutral 
solutions after removal of the aluminium as hydroxide. 

(4) Sodium (74Na) : counted as NaCl. Here no specific test was applied 
~ but a solution of a bombarded foil was treated to remove all cations 
other than the alkali metals and then evaporated to dryness. The 
residue showed a clear 14-5 hour decay characteristic of 24Na. 

(5) Magnesium (27Mg): precipitated and counted as Mg(NH POR 
care being taken to remove nitrogen and phosphorus activities. . 

(6) Silicon (748i) : a aaa and ae as SiO, from acid solutions, 

‘s Sik, after distillation of Sif’,. 

ee eee (0p and 22P): precipitated and counted as Zr (HPO a)2 : 
22P also as (NH,)3P0,. 12 MoOs, from solutions of bombarded foils in 
aqua regia. (For decay of a zirconium phosphate sample see fig. 4.) 
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(8) Sulphur (8) : precipitated and counted as BaSO,. 


(9) Chlorine (24Cl) : precipitated and counted as AgCl. The decay of. 


one sample was followed over twelve half-lives but the slope remained 
very constant, indicating again that the proportion of 38C] is quite small— 
probably below 5% (fig. 5). 

(10) Potassium (38K): precipitated and counted as KCIO,. 

Tn addition to the above, activities in calcium were looked for but not 
found; activities in argon, neon, oxygen and beryllium have not been 


looked for as yet. 


Fig. 5 
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Decay of activity of AgCl, showing *4Cl. 


3.3. Relative Yields of the Radioactive Products 


Several stacked-foil bombardments have been carried out for different 
bombardment times, and the relative yields of the main radioactive 
products have been calculated from the resolved gross decay curves of the 
individual foils. The relative yields are the same within a factor two for 
all the foils examined (down to 50 below the target surface), and are 
shown in table 1. They are expressed as numbers of atoms formed per 
100 of Cl. This isotope is chosen as an arbitrary standard as, owing to 


its suitable half-life and high yield, it was usually determined with greatest 
accuracy. 


EOE 
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Both “Mg and 11C must be formed in relatively low yield as their 
activities are not detected in the gross decay curves. However, if their 
yields (in numbers of atoms formed) were 10°% of the ®8K or 34Cl yields, 
phey could easily be missed as their half-lives lie between these two 
( : KK 7-5 min, “Mg 9-5 min, '4C 20 min, °4C133 min). Similarly, the yield of 

B must be considerably less than that of ®4Si (18F 112 min, 348i 157 min). 
It is in fact difficult to be sure whether 18F is a true reaction product at all, 
since some trial experiments in which various oxygen compounds were 
bombarded with N showed very strong yields of 8F, and it is possible 
that the *F observed in bombardments of aluminium foil arose from a 
surface layer of oxide. 


Table 1. Yield of Active Nuclei Observed and to be Expected on Various 
Assumptions, from the Bombardment of Aluminium with 4N 


Evaporation Buckshot 
Observable Evaporated @ Whether Group 1 Actual Yield 
nucleus particles (Mev) observed (?4Cl=100) 
Se) 300 
oP Pp 2a — 2-1 Yes 28 Al 140 
858 2p a — 26 Yes *4Na 12 
88 3p — 64 No Group 2 
eK 2n p — 90 Yes Baik 75 
ets 2p 2a —10:‘7 Yes eeCl 100 
ecAI p 3a —12:0 Yes ae 25 
2201 2n p a —16-4 Yes Group 3 
ce n 3p a —20-4 No aS es 
"Mg 2p 3a —21:6 Yes StS. 15 
evP. 2n p 2a —22:5 Yes "Mg small (<10) 
*4Na p 4a —22:9 Yes Group 4 
22Al n 3p 2a —31 No 22Na aay 
22Na 2n p 4a —42:3. No 2N Undetermined* 
2°Na n 3p 3a —44 No UG small 
nO 2n 6a —47-8 q The following 
aot 2n p 5a —50-8 2 should be rare 
a1¢ 2n Ta —58:1 Yes 0) 7 No other 
aa) 3n p 5a —59:9 No 1H Naw nuclei 
zo. 3n p 6« —65:7 Yes 28C lia 2°41) cobserved but 
eae unconfirmed 
possibility 
of ep 


+ See discussion in text of the special case of "Na. 

* The total yield of this isotope may be considerable, but the 8N nuclei appear 
to have a large velocity after formation, the amount found in any foil not being 
simply related to the amount produced therein. This point is being further 
investigated. 

§ 4. BoMBARDMENT OF Macnestum By “N Tons 


The lack of active calcium isotopes from the bombardment of aluminium 
was believed to be due to the lack of B-emitters of suitable half-life. It 
was realized, however, that it might be difficult to avoid the loss of one or 
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more charged particles from a compound nucleus with as small a 
potential barrier as 41Ca, so it was thought desirable to make a specific 
search for 38K in magnesium bombarded by nitrogen ions. The possible 
products of complete fusion of the 14N with the isotopic target nuclei are 
38K 3°K and “°K ; to give 8K one of the latter two must be able to get rid 
of its excitation energy by neutron loss alone. 

Four ten-minute bombardments were carried out and in each case the 
potassium, chemically extracted as before, showed an activity of about the 
expected half-life. In one case the charge discriminator was used to 
confirm that the particles emitted were positrons. Complete confirmation 
was thus obtained that all seven protons could be made to ‘ stick ’, though 
the yield was very small, less than 5% of the *4Cl yield. 


§5. OrHER EXPERIMENTAL EVIDENCE CONCERNING THE KINDS OF 
ReEactTION PRoDUCED BY NITROGEN IONS 


A preliminary examination has been made of some 250 nuclear inter- 
actions observed in an I}ford C2 emulsion bombarded by N ions of 
60-140 Mev. 

As was found by Miller (1951) using !*C ions, the commonest events 
apart from simple scattering are three-pronged stars. Most of the prongs 
are «-particles but a large number of our stars show one proton track. — 
With two possible exceptions, no cases of fission are observed in the 
emulsion either of the heavier compound nuclei formed from silver and 
bromine or of the lighter ones formed from carbon, nitrogen and oxygen. © 
(The fairly common elastic collisions which occur at low energy can readily 
be distinguished by the fact that the centre-of-mass energies are well 
below the Coulomb barrier concerned.) 

This suggests strongly that fission processes are rarely, if ever, involved : 
in the production of the active bodies identified in the work described - 
in previous sections. : 


§6. Discussion OF THE REACTIONS PRODUCED IN ALUMINIUM BY 
NirroGEn BOMBARDMENT 


In reactions produced by light particles of moderate energies, only a few 
reaction products are energetically possible. In the reactions described 
here, however, there is enough energy to produce almost every one of the 
known nuclides lighter than 4!Ca. 

In fact, by no means all of the energetically-possible nuclei appear in 
appreciable quantity and a particular feature stands out very strongly. 
This is the prominence of B~ emitters among the products (see table 1), 
particularly *°P, although the earlier Berkeley work on the bombardment 
of heavier elements with carbon ions showed a high probability of neutron 
evaporation (Miller et al. 1950). 

Let us consider some of the possible ways in which reactions may take 
place. The simplest assumption is that the whole of both nuclei first 
combine into a single excited compound nucleus, ‘1Ca. This will be des- 


- 


a 
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cribed as the ‘fusion’ nucleus to distinguish it from other possible 
compound nuclei involving parts only of the initial pair. The excitation 
energy can then be lost by fission or by evaporation of heavy particles. 
(This excitation energy, calculated from the “N energy distribution, will 
usually be in the range 35-50 Mev.) Evaporation and fission processes 
are not of course completely distinct ; fission here will be taken to mean a 
disintegration in which both fragments are heavier than alpha particles. 

In their paper Wyly and Zucker (1953) tacitly assume fission ; Miller 
(1951) appears rather to favour evaporation. We believe that fission may 
generally be regarded as improbable on the basis of the photographic 
emulsion evidence in § 5. It does not seem very likely that such reactions 
could be markedly more frequent with aluminium than with any of 
the elements common in photographic-emulsions. 

A multiple evaporation, particularly of «-particles, would be consistent 
with the photographic plate evidence and energetically is entirely possible. 
The production of *8K, for example, could well occur by loss of p+ 2n from 
41Ca, and that of *4Cl and °°P by the same process with loss of one or two 
_extra «-particles respectively. It seems curious, however, that the 
chief product of the N reaction should be *P, which involves the 
evaporation of two «-particles and a proton without loss of a neutron at all, 
and still more difficult to believe in the production of *1Si (two «-particles 
and two protons) or 24Na (four «-particles and one proton) by this means. 
Even if such a production were possible it is difficult to see why the yields 
(see table 1) should be so nearly comparable over the whole range from 
24Na to 38K. The fair yields of 4C and “N are further awkward facts. 

The evidence is strong, then, that even if evaporation processes account 
for the heavier positron emitters, they cannot alone be involved. 

The next possible process which will be considered is a simple * stripping 5 
process in which a neutron is removed by either nucleus from the other 
without significant interaction of the remaining parts of the nuclei. This 
process might well be favoured at low energies as the charged parts of the 
nuclei would not need to penetrate the full Coulomb barrier. This process 
would give rise to the following reactions : 


4NY 127A] > N-+26A] 
14) 127A] -> N+ 28A] 


possibly followed by 
26A] > 2Na+a 


28A] > 4Na+a 


which might account for N, Na, Na, 283A], There remain 1/C, ?’Mg, 
319}, 2P and *8 not fitting either of the above supposed processes. 

An explanation which fits the facts observed so far is as follows. It 
was noticed that the energy required to overcome the Coulomb barrier of 
aluminium in these reactions (17-4 Mev for 14N +-27A] in the centre-of- mass 
system) was comparable with the energy required for the disintegration 
4 > 3a--n-+p, which is 19-77 Mev. Let us suppose that the incident 
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14N nucleus might be thought of, at least in the field of target nucleus, as 
three independent «-particles, a proton and a neutron. The target could 
then be visualized as being struck, not by a single bullet, but rather by a 
shower of buckshot. (If the nuclei are of comparable mass, of course, 
both might be regarded at the instant of collision as composed of largely 
independent particles.) 

It is not intended to imply that the undisturbed 4N must be already 
divided up in this way or that the paths of the individual particles are 
parallel. The internal velocities of the particles must be quite comparable 
with the velocity of the whole nucleus. The separation visualized may be 
regarded as a virtual one only, unless one or more of the Apt e is 
actually captured by the target nucleus. 

If we consider this picture, it is easy to see that only when the impact 
parameter is small will there be a good chance of all the separate particles 
being captured by the target nucleus. If the impact parameter is large, 
some particles may be captured but others will miss and be deflected away. 
On simple geometrical grounds it can then be seen that the chances of 
capture of one, two or three «-particles would all be comparable. - 

A difference of behaviour between the «-particles and the odd proton 
and the odd neutron is to be expected, for the neutron will find no potential 
barrier against its capture. The deflecting force in fact will be even more 
important for the proton than it will be for an alpha particle, which for 
a given velocity has four times the kinetic energy and has to surmount only 
twice the potential barrier. The groups of particles escaping capture by 
the aluminium, then, would have the following order of probability : 

(1) wx+p, 
(2) ax«+n-+p, 
(3) xa, 
(4) va-+n, 
where x may be 0, 1, 2 or 3, though z=3 may be regarded as a special case. 

It is reasonable to suppose, where considerable energies are involved, 
that radiative capture is less likely than capture with subsequent emission 
of a heavy particle. Potential-barrier considerations would favour 
neutron emission but discourage proton emission. In the region of the 
periodic table concerned, the binding energy of either nucleon is fairly 
large, often larger than that of an alpha particle. The energy per nucleon 
needed by an alpha particle is thus a good deal less. We shall thus regard 
the emission of either a neutron or an alpha particle as more favoured than 
emission of a proton. It is worth noting that this will be particularly 
important where zx-+-n have just been captured, as the energy required to 
extract a proton will then be higher than that required for either of the 
other particles. 

On this basis, then, we can readily explain the large production of 32P, as 
follows : 

27 Al+n-+ 2a —> 36C] -> 32P+ 
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The complete set of this class of reactions, with the favoured modes of 
_ breakdown, are as follows : 


a77“Al+n—7:7 Mev (The reaction- 
‘ 24Na-+a%— 10-9 Mev energies listed 


te 731P4+n—7-9 mev refer only to the 
\\28A]1¢—9-9mey Second steps of 


Pe re tae Ueun SHanee these processes.) 
‘\ 82P %—8-3 Mev 

7*®*K +n—7:8 Mev 

‘ 86C] + «%— 5-6 Mev 
There are thus three observable radioactive products, *4Na, ?8Al and 22P. 
In the parent of *?P alone « emission is energetically favoured as against 

neutron emission, which may well account for the yield being greater than 
that of 8Al. For the first case, 24Na, it may well be that the capture of a 
single neutron, from a bound state of similar energy, does not usually 
leave enough excitation energy for « emission and either gives ?’Al or 
contributes directly to the *8Al yield. 

The group of nuclei produced by adding xx particles to aluminium 

- comes next, giving us ?°K, ®Cl and 31P. ‘The previous situation is now 
reversed ; emission of an «-particle in each case leaves a stable nucleus 
while emission of a neutron gives *°K, **Cl and °P. Table 1 shows that 

these are in fact the next group experimentally, with quite similar yields. 

The next group, in which x«-++-n-+p are captured, gives us ?°Si, 39S, 397A 
and 41Ca. Emission of either a neutron or an alpha particle from any of 
‘these leads to a stable nucleus except for ?’A which is a K-capturer and 
not observable by the methods used. In the final group, xx-+-p are 
supposed to be captured. This gives rise to *8Si, *S, °°A and “Ca. In 
every case « emission gives rise to a stable isotope and n-emission to an 
isotope too short-lived for our observation. 

We have next 2’Mg, 348i and *S, each of which requires a proton 
emission from one of the readily-produced Al+-a«-+n nuclei. It will be seen 
that the yield increases more rapidly with mass than in the other series. 
This may reasonably be regarded as due to the higher excitation needed 
for proton emission, which is more likely to occur when more particles have 

been captured. Finally, we have 3N and ''C. These may arise from 
the uncaptured particles, when these are not permanently separated. 

Tn table 1 are shown all of the observable nuclei which could be produced 
by evaporation, with the corresponding values of Q. There is some 
correlation between the Q’s for evaporation and the observed yield but 

- there are some striking discrepancies such as the large yield of **Cl and the 
* s0-far unobservably small yield of Cl. The same table shows the nuclei 
which would be expected on the buckshot picture, in the expected. order of 
yield, and the observed values. No attempt is made to predict the relative 
probabilities of the inverse reactions ; i.e. those in which nucleons are 
captured by the lighter nucleus from the heavier. Most of the nuclides 


27Al+n—> A] 


7Altotn > 


27 Al+3a-+n > °K 
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produced in this way would be stable and unobservable, and the photo- 
graphic emulsion evidence suggests that reactions leaving two heavy 
nuclei are rare. — 

It will be seen, however, that an excellent qualitative fit exists between 
the nuclei predicted by the proposed ‘ buckshot’ picture and those 
observed. Although no observed nucleus requires the emission of more 
than one heavy particle from the primary compound nucleus, it is not 
suggested that multiple emission cannot occur at all. 

Thus two neutrons might well be emitted by the first group of compound 
nuclei mentioned (28Al, ?2P, 36Cl and #°K), adding to the yield of 9°P, *4Cl 
and 38K obtained from the second group. It must be remembered, how- 
ever, that in the lighter compound nuclei the level of excitation will not 
ordinarily be very large. The more-or-less exponential energy distri- 
bution above 50 Mev of the 14N ions in the cyclotron, already mentioned, 
ensures that, at any point within the target, the majority of the ions 
energetic enough to cause a reaction have energies not very far above the 
barrier. This explains why the relative yields of different products seldom 
vary very much with depth in the stack of foils. An estimate based on the 
energy spectrum of the ions indicates that almost three-quarters of the 
ions reaching the target with energies over 26 Mev are still below 45 Mev. 
A 14N ion with 45 Mev energy in the laboratory system will carry into the 
compound nucleus an energy of only about 2 Mev per nucleon in the centre- 
of-mass system. To take a specific instance, consider the production of 
36C by capture of two alpha particles and a neutron from a #N ion with 
20 Mev available energy (i.e. 30 Mev in lab. system). Of this, since 9 
nucleons are captured, on the average about 13 Mev will go towards exciting 
the compound nucleus. Since the Q of the reaction 

27Al+ 4N—p—« > 36C] 
is only about 6 Mev, the total excitation of the °°Cl will be about 19 Mev. 
This is only a little above the binding energy of two alpha particles (16 Mev) 
and is actually below that of two neutrons (22 Mev). At such an excita- 
tion two-particle emission is, therefore, likely to be uncommon even from 
this relatively heavy compound nucleus. It can probably be neglected 
altogether in the lighter ones. 

This slightly surprising result is, of course, due to the fact that the 
particles are captured from a bound state in the N nucleus of potential 
energy not very dissimilar to that in the compound nucleus. 

While not ruling out the possibility of fusion followed by evaporation 
for the heavier nuclei formed, or of neutron-stripping for the nuclei close 
to the original ones, the partial capture picture seems satisfactorily to 
cover the whole range of products ; it may be noted, of course, that the 
other possibilities in fact arise naturally as limiting special cases of this. 

Now if there is any substance in these considerations, an immediate 
experimental test suggests itself. If we use as our bombarding nucleus, 
instead of 4N, a pure «-particle nucleus such as !2C or 160, a very different 
set of products should appear. 1208+ is less convenient as it is very difficult 
to avoid heavy contamination from the cyclotron with 2D+ or 4He?+, but 
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the resonance of 160% lies well away from any possible interfering particles 
and this was chosen. 


In the next section the results obtained with this are discussed. 


$7. BOMBARDMENT OF ALUMINIUM wrTH 160 Ions 


The best oxygen beams at present obtainable are some 50 times less 
than the best nitrogen beams, but are nevertheless quite adequate to give 
fair activities. 

Several bombardments of various durations have been carried out. 
A very marked difference from the nitrogen decay curves was observable. 
The oxygen decay curves could be analysed perfectly clearly into : 


(1) a weak 14 day activity (°?P), 
(2) a weak 15 hour activity (?4Na), 
(3) a weak 112 min activity (48F), 
(4) a 33 min activity (°4Cl), 

(5) a 74 min activity (°8K), 

(6) a 24 min activity (?°P). 


No activities longer than 14-days half-life have yet been observed, although 

this would not rule out the possibility of ?*?Na being formed in yield 
comparable to that of the other nuclei mentioned, as its counting rate 
would be very low owing to its long half-life. Chemical experiments 
proved conclusively that the 112-minute decay was due to SF, as a weak 
but definite activity of this period was found to be associated with a 
precipitate of calcium fluoride obtained from a solution of a bombarded 
foil in hydrochloric acid after addition of suitable carriers. On the other 
hand, a precipitate of silica was completely inactive indicating the absence 
of 218i among the reaction products. These results are in contrast to those 
obtained with 14N bombardments, in which the *48i activity was much 
stronger than that of *F. Separations of zirconium phosphate demon- 
strated the 24 minute activity as being due to *P. The magnetic 
discriminator showed that the activities ascribed to *4Cl, °K and *°P were 
in fact due to positrons ; any real negative electron activity could not 
have been more than 5% of the positron activity observed in each case. In 
view of the clear resolution of the decay curves we have therefore 
considered further identifications unnecessary. 

The relative yields are shown in table 2, again taking 34Cl as 100. 

It will be seen that, as expected in the ‘ buckshot ’ picture, there is a 
complete change-over of the relative activities from the nitregen yields. 
As we can no longer add anything but zx, and emission of protons or 
a-particles from the resulting nuclei leads always to stable products, to a 
first approximation we should expect the only observable reactions to be : 


27Alto > 31P > %P+n 
2741-4 20 > 801 > #Cl-++n 
27A1+ 3a > 9K > 8K+n 


26A] and 42Sc being unobservable. 
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Table 2. Active Isotopes Observed to be Produced from Bombardments 
of Aluminium with Oxygen 


Active Nucleus Yield relative to 
34C]=100 
38k Hes: 
34(] 100 
30P 35 
18h" ky 
32pP 5 
24Na, 3 


No other product observed 


In fact these are not the only, but are outstandingly the most important, 
products. The very small yield of *P may involve the emission of the 
less-likely particle >He, which is known to be produced from various 
high-energy reactions. The yields of 1*F and *4Na may be made by the 
inverse reactions in which 140 abstracts the stray nucleons from the ?’Al, 
e.g. 

27 A] }-16Q > 4M o 19 
19 > 18h 4+n 
(none of the simple particles emitted from *4Mg will leave an observable 
product). 


27 A]+- 16Q) -+ *5Mo+ 18H 
25Mo > *4Na+p 
the small yield of “F produced from 18F not being observable in the 
presence of the heavy yield of ®°P. 


§ 8. CONCLUSION 


The great contrast between the N and 1%O results, which appears as 
a natural consequence of the *‘ buckshot ’ hypothesis, is very difficult to 
explain on the basis of fission and not immediately easy on the basis of 
evaporation from the fusion nuclei 4'Ca and Sc. It would be surprising 
that evaporation at moderate energies from the former should lead to such, 
an excess of 8~ emitters and from the latter of 8+ emitters. 

We believe that the data described here, while they are far from proving 
that the ‘buckshot’ picture is correct, do show that it is a helpful 
conception when trying to see a pattern in the observed facts, and further 
that it may be useful in predicting the likely outcome of this very compli- 
cated kind of nuclear reaction. The results predicted do not differ very 
greatly from those to be obtained from the evaporation model if, after the 
first particle emitted, evaporation of « particles is very strongly preferred, 
but, where they do differ, our picture seems to be in closer accord with the 
experimental results. 
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SUMMARY : 


A method of separating nuclear disintegrations caused in photographic 
emulsion by 50-125 Mev protons into those in the light and those in the 
heavy nuclei of the emulsion is described. It is used to determine the | 
characteristics of the disintegrations of the light and the heavy nuclei. 
These results show that about thirty per cent of the low energy particles 
emitted from nuclear disintegrations are the result of a nuclear cascade 
initiated by the primary particle. 


A ae 


INTRODUCTION 


In a previous paper (Hodgson 1953, referred to as I) the results of measure- 
ments on nuclear disintegrations caused by 50-125 mev protons in Ilford 
C2 emulsion were described. These measurements refer to the 
disintegrations in the mixture of elements constituting the emulsion, and 
are consequently difficult to interpret theoretically. In Part I of 
this paper, a method is described by which it is possible to separate, with 
quite a high efficiency, the disintegrations reported in I into those in the — 
light (carbon, nitrogen and oxygen) and those in the heavy (silver and — 
bromine) nuclei of the emulsion. 

Part ITsummarizes the main characteristics of the nuclear disintegrations | 
in the light and heavy nuclei determined by measurements on events 
separated according to the procedure described in Part I. In Part III 
several models of the collision and disintegration process are compared 
with the experimental results for the heavy elements. 


Part I.—TuHer SEPARATION PROCEDURE 


$1. InrTRODUCTION 


The method of separation uses some well-established characteristics of 
nuclear disintegrations in the light and heavy elements to divide the 
observed events into the two categories. Several possible characteristics 
were tried before the energy of the emitted alpha-particles was finally 
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selected as the most suitable. In addition, several subsidiary criteria 
were used to allocate those events from which no alpha-particles were 
emitted. 

The first attempt was made with the method based on the length of the 
shortest track which had previously been used with success to separate the 
disintegrations caused by the cosmic. radiation (Hodgson 1952b). This 
method depends on the results of Harding (1949 b) who showed that most 
of the stars caused by the cosmic radiation and having recoil fragments 
are due to the disintegration of the heavy elements, while most of the 
stars without recoil fragments are formed in the light elements. When 
this method was applied to the present group of stars it did not effect the 
desired separation. This is because the average energy of the disintegra- 
tions is much less than that of those studied by Harding, with the result 
that a considerable proportion of the stars in the heavy elements have no 
visible recoil fragment. 

The separate procedure used is described and justified in the following 
section. 

§2. THE SEPARATION CRITERIA 


The separation procedure used here is based essentially on the fact that 
very few, if any, alpha-particles emitted from moderately excited silver 
and bromine nuclei have energies less than 9 Mev. 

This is to be expected from the potential barriers of silver and bromine, 
which are about 14 and 12 mev respectively in the unexcited state, when 
allowance is made for slight lowering of the barrier due to the nuclear 
excitation. It has been confirmed by calculations of the emission process 
by Puppi (1953), who found a sharp low energy cut-off at about 9 Mev 
in energy spectrum of the alpha-particles from such nuclei. The first 
experimental proof of this was given by Perkins (1949), who showed by 
subtracting the energy distribution of the alpha-particles from light 
nuclei (obtained by the gelatin sandwich method) from that of a corres- 
ponding number of disintegrations in normal emulsion that “ practically 
all alpha-particles of energies below 10 Mev can be attributed to light 
nuclei’. This conclusion was confirmed by Menon, Muirhead and 
Rochat (1950), who measured the energy distribution of the alpha- 
particles from one-track stars caused by the nuclear absorption of slow 
a--mesons and found two peaks, one at about 2 Mev corresponding to the 
light elements and one at about 10-12 mev corresponding to the heavy 
elements, with a marked minimum in the range 8-10 Mev. Further 
support is provided by the measurements of Lees, Morrison, Muirhead and 
Rosser (1953) and of Hodgson (1953) who find a peak at about 12 Mev in 
the energy distribution of alpha-particles emitted from one- and two-track 
stars caused by 110-146 and 50-125 mev protons respectively, and by 
measurements of Hodgson (1952) of the energy distribution of the 
alpha-particles from disintegrations caused by the cosmic radiation. 

To establish a rigid set of criteria, it is necessary to distinguish between 
tracks due to nuclear recoils and those due to alpha-particles of very low 
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energy. These cannot be distinguished by their ionization, and so this 
must be done by their range. The range distribution of short tracks from 
all the disintegrations was therefore determined, and it was found to 
consist of a sharp peak at very small ranges, which may be ascribed to 
the nuclear recoils, superimposed on a broad distribution extending to 
greater ranges which may be ascribed to the low energy alpha-particles. 
The tracks comprising the sharp peak were all less than four microns long, 
which agrees with approximate calculations of the range of nuclear 
recoils. All particles making heavily ionizing tracks longer than four 
microns were therefore counted as alpha-particles. After the stars had 
been assigned to the light and heavy elements, tracks shorter than four 
microns in stars in the heavy elements were counted as nuclear recoils, 
while those in stars in the light elements were counted as alpha-particles. 

All events with one or more tracks of alpha-particles of energy less 
than 9 Mev were assigned to the light elements, and events with all 
alpha-particles of energy more than 9 Mev to the heavy elements. As 
alpha-particles are very frequently emitted from nuclear disintegrations, 
the majority of the events were allocated by means of this criterion. 

The energy distribution of the protons from the remaining events was 
then compared with that of the protons from the stars already assigned 
to the light or to the heavy elements. It was found to be indistinguishable 
from that of the protons from the heavy elements and significantly 
different from that of those from the light elements. All the remaining 
events were therefore assigned to the heavy elements, with the exception 
of two types of events noted below. 

Firstly, events with two or more tracks with range less than four 
microns are almost certainly due to the disintegration of a light nucleus, 
as fission of the residual nucleus from the disintegration of a heavy 
nucleus is most unlikely at such low excitation energies. Secondly, 
although the emission of low energy protons from heavy nuclei may be 
explained by the process discussed by Le Couteur (1950) by which there 
is a probability of 0-3 that the last proton evaporated from an excited 
heavy nucleus has an energy between 0-5 and 3-5 Mev, the emission of 
two such protons from a heavy nucleus is very unlikely. Such events 
were therefore assigned to the light elements. 

Tn this way each event was assigned either to the light or to the heavy 
nuclei. After this had been done, the tracks made by nuclear recoils 
from the disintegration of heavy nuclei were no longer included in the 
star size. Thus a star of two protons in a heavy element is called a 
two-track star whether it has a visible nuclear recoil or not, in contrast 
to the classification in I, in which all tracks, including nuclear recoils, 
were included in the star size. The advantage of the new classification 
is that the assignment of star size no longer depends on whether a very 
short track is observed or not. As these tracks are difficult to see, and 
their visibility depends on the type and degree of development of the 


emulsion, this should facilitate the comparison of results obtained by 
different observers. 
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The separation procedure outlined above is not infallible but it is 
sufficiently reliable for a determination of the main features of the 
disintegrations of light and heavy nuclei to be made. Its efficiency is 
clearly greater for some types of events than for others, and this must 
be taken into consideration during the evaluation of results. Thus, for 
example, it is very reliable for stars of four or more tracks but quite 
unreliable for stars consisting of only one proton. 


§3. COMPARISON WITH OTHER SEPARATION PROCEDURES 


The method described in the last section may be compared with two 
others designed to achieve the same object. Lees, Morrison, Muirhead 
and Rosser (1953) have compared the disintegrations produced by 
130 Mev protons in normal G5 emulsion with those in emulsion containing 
more than the usual proportion of gelatin. By suitable subtractions they 
find the characteristics of the light and heavy elements separately. The 
main advantage of this method over the present one is that there is less 
possibility of systematic error, but this is offset to some extent by the 
low statistical accuracy resulting from the subtraction procedure and the 
extra labour of calibrating and processing multigel plates. 

Blau, Oliver and Smith (1953) and Germain (1953) have applied the 
gelatin sandwich method of Harding (1949a) to the study of the 
disintegrations of the light elements caused by artificially accelerated 
particles. This method has the advantage that every event observed to 
originate in the gelatin layer certainly occurs in a light element. One 
of its disadvantages is that the number of events per unit volume of 
emulsion is much lower for comparable conditions of exposure than 
that obtained by other methods, which increases the scanning time 
per event. Further, some particles, especially those of short range, 
remain wholly in the gelatin layer and must be allowed for in the 
evaluation of results. 2 

These three methods are in many respects complementary, and the 
best one to use depends on the problem being studied. The main 
advantage of the method described here is that results can be obtained 
more rapidly than by either of the other methods. 


Part II.—THE DISINTEGRATIONS OF THE 
Ligut anD HpAavy NUCLEI 


$1. Tue Size DISTRIBUTIONS 


The distributions of star sizes are given in table 1 for the four energies 


of incident protons. 

These distributions are considerably more peaked than those calculated 
from Poisson’s Law. They are remarkably insensitive to the energy of 
the incident protons, and the mean star size increases slowly with 


increasing £ ,. 
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§2. Tur RELATIVE NUMBERS OF SECONDARY PROTONS AND 


ALPHA-PARTICLES 


The proportion P, of alpha-particles among the emitted particles from 
stars in the light nuclei showed no significant variation with H, for stars 
of the same size. P, is plotted in fig. 1 as a function of star size. 


Table 1 


Number of tracks 


Light Nuclei Heavy Nuclei 


3 


34 
25 
14 
20 


4 


40 
Li 
21 
21 


UPPER AND LOWER LIMIT CURVES 


— NO FAST PROTONS 
EMITTED 


---- ONE FAST PROTON 
. EMITTED PER STAR 


Star size 


Variation of P,, with star size for the disintegrations of the light nuclei. 


Upper and lower limits to P, can be calculated from elementary consider- 
ations. ‘The upper limit is set by the number of charged nucleons in the 
target nucleus and is the reason why P, decreases with increasing star size. 
If we assume that the nucleus disintegrates completely into alpha-particles 
and protons only, then the numbers of alpha-particles and protons, and 
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hence P,,, are fixed by the charge on the target nucleus and the size of the 
star. The upper limit was calculated in this way as a function of star 
size by combining the disintegrations in carbon, nitrogen and oxygen in 
the proportions in which they are present in the emulsion and assuming 
that the cross-section for star production follows an A? law. 

The lower limit to P,, is set by the fact that if a light nucleus disinte- 
grates into a few particles, at least one of these must have charge greater 
than or equal to two, and so be classed among the alpha-particles. Thus 
[Palmin=1/n for an n-track star whenn<5. It is of course possible that 
the heavily charged particle makes too short a track to be visible, but this 
happens so rarely that the lower limit given by the formula above is not 
materially affected. 

These calculations will be affected if one or more fast secondary particles 
were emitted and not counted among the tracks of the star. They were 
therefore repeated assuming that one such particle was emitted from each 
star. 

The results of these calculations are plotted in fig. 1 and are consistent 
with the experimental results. 

The proportion of alpha-particles among all the particles emitted from 
heavy nuclei was found to have no significant variation either with the 
energy of the incident proton or with the star size. Its mean value over 
all the stars is 0-22+0-02. 


§3. Tue PROPORTIONS OF STARS WITH DIFFERENT NUMBERS OF PROTONS 
AND ALPHA-PARTICLES 


The proportions of stars with different numbers of protons and alpha- 
particles were found to be independent of the energy of the incident proton. 
Results for all energies are therefore grouped together in the following 
table. 


Table 2. Numbers of Stars with V Alpha-particles 


Light Nuclei Heavy Nuclei 
4 2 4 


2 3 
oe) 


6 
3 
2 
1 


2 
1 
0 


These distributions were compared with those given by Bernouilli’s 
Law, which corresponds to the case when the emission of one particle does 
not affect the probability of emission of another similar particle (Hodgson 
1952 a), and no significant difference was found between them. 
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$4. Tur ENERGY DISTRIBUTIONS OF THE SECONDARY ALPHA-PARTICLES 
AND PROTONS 


The energy distributions of the secondary alpha-particles and protons 
were found to have no significant variation with the energy of the incident 
protons. The variation of the energy distribution of the alpha-particles 
from the light elements with star size is shown by figs. 3 (a)-3 (d) of Lif the 


(b) Backward Hemisphere 


(2) Forward Hemisphere (b) Forward Hemisphere 


fre a ve 


(c) Backward (c) Forward 
Hemisphere 10 Hemisphere 


15 20 5 auto 
(d) Backward (e) Backward 
Hemisphere Hemisphere 


Corrected number of particles 


6 8 2 4 
(¢) Forward (e) Forward 
Hemisphere Hemisphere 


Energy in Mey. 


Energy distributions of protons and alpha-particles from light and heay 
nuclei. (a) Alpha-particles from 1 to 3 track stars in light aera 
(6) Alpha-particles from 4 to 7 track stars in light nuclei. (c) Al ha- 
particles from heavy nuclei. (d) Protons from light nuclei. (e) Pride 
from heavy nuclei. -—--~- Fitted evaporation curve. 
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peak at 10-12 mev in fig. 3 (a) due to alpha-particles from the heavy 
elements is omitted. The energy distributions of the alpha-particles 
emitted in the forward and backward hemispheres with respect to the 
incident particles are shown in fig. 2 for the light and heavy nuclei. The 
energy distribution of the alpha-particles emitted in the backward 
direction from the heavy nuclei was fitted by the best curve of the form 


(= | LW Ly 
Bape al aa exp | - wer | aT’ 


TOr*! 7 


where V’ is the effective potential barrier height and +* the nuclear 
temperature, which is approximately equal to 11/12 7 (Le Couteur 1952). 
The experimental values were V’=11 Mev, 7r=2 Mev. 

These results can be used to check Le Couteur’s relation U=A 72/11 
between the excitation energy U andthe nuclear temperature r. Inserting 
the known values of 7 and A, we obtain U=34 mev which is of the correct 
order of magnitude for a group of disintegrations caused by 50-125 mev 
protons. 

The alpha-particles emitted in the forward hemisphere contain a 
considerable proportion of direct knock-ons (see Part III) and so they 
were excluded from the above comparison. 

The energy distributions of the protons did not vary significantly with 
star size and are shown in fig. 2 for all stars together, the distributions in 
the forward and backward hemispheres being shown separately. They 
were found from those protons ending in the emulsion, and so do not 
extend beyond 10 Mev. 

Tn all cases the usual corrections for particles passing out of the emulsion 
were made. 


§5. Tue ANGULAR DISTRIBUTIONS OF THE SECONDARY ALPHA-PARTICLES 
AND PROTONS 


The angular distributions of the secondary alpha-particles and protons 
were found to have no significant variation either with the energy of the 
incident protons or with the star size. The results for all stars together 
are therefore plotted in fig. 3. The particles ending in the emulsion and 
those leaving it before coming to rest, which correspond roughly to the 
lower and higher energy particles, are shown separately in these figures. 
The angle 7 plotted is the projection on the plane of the emulsion of the 
angle between the direction of emission of the secondary particle and 
that of theincident protons. Isotropic emission therefore corresponds to a 


uniform distribution of %. 


§ 6. THe Proportion oF Stars IN THE Heavy NucLEr 


The proportion of all stars that are assigned to the heavy nuclei is 
plotted in fig. 4 as a function of incident proton energy for each star size, 
The lines are drawn through the experimental points to emphasize the 


— 


TO t and heavy © 
n he einen . 
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main trends but, owing to the considerable statistical uncertainties, their 
positions are known only approximately. It can be seen that the 
proportion of stars in the heavy nuclei increases with the energy of the 
incident protons for each star size. The results of Blau, Oliver and 
Smith (1952), who investigated the disintegrations of light nuclei by 
300 Mev neutrons using gelatin sandwich emulsions show that this trend 
continues to higher energies, for they found that only 17% of all the 
events occurred in light nuclei. These results also agree with those of 
Lees, Morrison, Muirhead and Rosser (1953) who examined disintegrations 
produced by 130 Mev protons. 

This diagram shows that a considerable preliminary separation of the 
events into those in the light and those in the heavy nuclei may be 
effected by selecting the incident proton energy and the size of the stars 
examined. 

Part III.—INTERPRETATION 
§ 1. INTRODUCTION 


In this Part, the results given above are compared with those 
predicted by simple models of the collision and disintegration process. 

The simplest model of the process is to suppose that the incident proton 
enters the target nucleus and is brought to rest inside it by a series of 
nucleon—nucleon collisions, thus transferring all its energy and momentum 
to it. The target nucleus then moves forward with a momentum equal 
to that of the incident proton, and evaporates secondary particles until 
it returns to the ground state. 

If this were so, the excitation energy of the nucleus would be propor- 
tional to the energy of the incident proton, so that the mean energy of 
the emitted particles or the mean star size or both would increase markedly 
with primary energy. This is in disagreement with the results of Part IT, 
which indicate that the primary particle often emerges from the target 
nucleus with a large fraction of its initial energy. This is confirmed by 
observation of the tracks of high energy particles emitted from many 
of the disintegrations and by numerous experiments on the attenuation 
of beams of particles through various materials. 

The simplest development of this model is to suppose that only the 
primary particle goes on with high energy, while all the other particles 
are emitted at low energy by an evaporation process. The ability of 
this model to account for the observed features of the disintegration of 
the heavy elements is investigated in the next section. 


§2. Tue SzeconD Move or THE COLLISION PROCESS 


Suppose that the primary particle gives the target nucleus sufficient 
momentum to give it a velocity V, and let P(v) dv be the number of 
secondary particles with velocities in the range v to v+-dv, with 


, P(v) dv=l1. 
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Then it may easily be shown that the forward fraction F for the low- 
energy particles is given by 


P= 5. (1+ =) P(e) dot { Pe) dv. 


In the case of the heavy elements V <v, and so F is very insensitive 
to the assumed value of the momentum transfer. If rough values of 
V and v are inserted in the above relation, it may be shown that the 
mean value of F' lies between 0:50 and 0-51, and even if extreme 
assumptions are made does not much exceed 0-52. 

These figures are inconsistent with the observed values of 0-66=- 0-06 
and 0-63-0-03 for alpha-particles and protons respectively, which 
indicates that a considerable proportion of.the low energy particles are 
not emitted by an evaporation process. A modified model of the 
collision process by which this can be taken into account is discussed in 
the next section. 


§ 3. Mopirrep MopEet with NucLEAR CASCADE 


Suppose that the incident proton initiates a nuclear cascade in the 
target nucleus, as a result of which a number of neutrons, protons and 
alpha-particles are emitted. When the cascade has finished, the nucleus 
returns to the ground state by evaporation as before. 

This model of the collision process is similar to one proposed by 
Bernardini, Booth, and Lindenbaum (1952) to account for their observa- 
tions of the disintegrations of heavy nuclei caused by 350 Mev protons. 
These investigators, however, used G5 emulsion and so they were unable 
to distinguish between low energy alpha-particles and protons. They 
were therefore unable to establish the participation of alpha-particles in 
the cascade process and so in their calculations they take account only of 
the neutrons and protons. The present results, however, provide strong 
evidence that the alpha-particles take part as well. This can be 
interpreted as implying the existence of alpha-particles in the nucleus, 
orit may be due to the nucleons of the cascade picking up further nucleons 
on their way out of the nucleus. 

The detailed calculations of this cascade process are best treated by the 
method of Goldberger, but a few preliminary results can be obtained by 
making some simple assumptions. 

Assume that all the particles of the nuclear cascade are emitted in the 
forward hemisphere in the laboratory system. This is likely to be very 
nearly the case, although it is possible that a few particles are emitted 
in the backward direction as a result of multiple collisions within the 
target nucleus. As a result of this assumption, the calculated percentage 
of all the emitted particles that are emitted in the cascade is too low. 
The forward fraction for the evaporated particles can, without sensible 
error, be taken as 0:5. 

Let C be the proportion of charged cascade particles among all the 
emitted charged particles, and let p,“ and p,” be the proportions of 
alpha-particles among the cascade and evaporated particles respectively. 


Caused by 50-125 Mev Protons : II 201 


Then we have 
eee eae 0) ee 
a Cpa l= C)p,2 : 
ee OL Pye Ope) 
- CUE Ps ae tC) (1 2") : 
P,=Cp,°+(1—C)p,”. 

Substituting the experimental values of F,, F,, and P,, we obtain 
C=0-26+0-06 ; p,°=0-27-0-14 and p,”=0-20+0-05. Thus about 25% 
of the low energy particles are emitted as the result of a nuclear cascade 
in the target nucleus. This is probably a little low for reasons pointed 
out above. The results also show no significant difference between the 
proportions of alpha-particles among the cascade and evaporated 
particles. These results may be compared with those of Bernardini, 
Booth and Lindenbaum (350 Mev protons) who estimate that “at least 
25% of the black prongs from AgBr are knock-ons” and with those of 
Blau, Oliver and Smith (300 Mev neutrons) who conclude that ‘‘ at most, 
70% of black prongs in heavy elements are due to nuclear evaporation ”’. 

The existence of the nuclear cascade means that the velocity of the 
excited nucleus cannot be calculated from the forward fraction /’, although 
this has often been done. 

These calculations cannot be applied to the disintegrations of the 
light nuclei because the number of nucleons in the nuclei is so small 
that the distinction between cascade and evaporation particles can no 
longer be made. 
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SuMMARY 


The antisymmetry of the total wave function of a Fermi—Dirac gas leads 
to an apparent repulsion between the gas particles. In other words, the 
average probability that two particles will be found at a distance r from 
one another falls to zero as r tends to zero. This pair distribution function, 
Dir), is calculated (§ 2) at all temperatures by means of a variation method 
founded on a minimum property of the free energy. It is not possible to 
make a straightforward calculation using Fermi—Dirac statistics on account 
of the integrals which arise. Instead we use an approximation which has 
already been employed to meet a similar difficulty in the problem of 
finding the influence of exchange energy on the specific heat of electrons in 
metals. 


§ 1. INTRODUCTION 


Iv is well known that the distance correlations between the molecules of an 
ideal Fermi—Dirac or Bose-Einstein gas are due to the symmetry properties 
of the wave functions (see for example the review article by de Boer 1949, 
especially § 6(iii)). One may summarize the known results by saying that 
there is an apparent attraction between Bose—Hinstein particles and an 
apparent repulsion between Fermi—Dirac particles. Thus Uhlenbeck and 
Gropper (1932) showed that in the limit of high temperatures and low 
densities, the probability of finding a particle within a distance interval 
(r, r+-dr) of a given particle is 47r2D(r)dr, where 


D(r)= > (Lex (—21%/X%)}, Mee ey ee TEM, 


in which the + sign is taken with Bose—Kinstein statistics and the — sign 
with Fermi—Dirac statistics. Here N denotes the number of particles in a 
volume V and A is the de Broglie wavelength corresponding to the 
temperature 7’, 
A= h/(27mkT)12, 
At lower temperatures it is no longer possible to handle both cases together 
and in the present paper we shall consider only a Fermi—Dirac gas. The 
Bose-Einstein case has been treated by London (1943). 


SS ee 
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The expression for D(r) in the opposite limit of very low temperatures 

and high densities was obtained by Wigner and Seitz (1933). 
Din 741-3 (= cos kyr—sin “)"} 

Vi (kor)? ‘ 
where ky=7(6N/7V)1/3. The one-particle wave-functions % are plane 
waves normalized over a volume J, i.e. 4,=V-'? exp i(k. r), and’ the 
quantity k) is the wave number of the highest occupied level at absolute 
zero.* It is related to the Fermi energy «) through the equation 
€p= hk?) /82*m. 

Now although we have these expressions (1.1) and (1.2), valid at very 
high and very low temperatures respectively, there are no calculations 
which bridge the gap between these limits. Bhatnagar and Singwi 
(1949) have, however, shown that (1.1) and (1.2) are the leading terms in 
two different series and have calculated additional terms in these 
expansions, thereby approaching the region of intermediate temperatures 
simultaneously from both limits. Their starting point is a general formula 
for D(r) namely 


mo Bll enteey pel). = 0 


where f, is the average occupation number of the state 4. (As usual 
we shall replace these sums by integrals, assigning Vdk/(27)* states to 
each volume dk of k-space.) At the present time it is not possible to 
evaluate (1.3) at all temperatures, if for f, we use the exact Fermi—Dirac 
distribution function. In particular the region of intermediate temper- 
atures around kT/ey~1 is covered neither by the high nor by the low 
temperature series. It therefore seemed of interest to calculate D(r) 
using an approximate distribution function, one which would allow closed. 
expressions to be obtained at all temperatures. Previous work (Lidiard 
1951 a and 1951 b, particularly § 2 and appendix A) has already shown 
how even quite a crude form for /,, can give reliable results if the specifying 
parameters are chosen so as to minimize the free energy (Lidiard 1953). 
The details of the calculation are described in the next section and 
numerical results are exhibited graphically in fig. 2. 


(1.2) 


§ 2. CALCULATION OF D(r) 
First we discuss the approximate distribution function to be used. 
We note, to begin with, that the free energy function, F, of an ideal 
Fermi-—Dirac gas may be expressed in terms of the distribution function as 


follows : 
hy (° RTV (® 
pa EE [pay meant Sy J {fiastapin aphid. — @ 


* We are assuming that only one particle at a time can be in the state #,. 
It is a trivial matter to generalize the equations to the case where the occurrence 
of spin degeneracy allows more than one particle per state. 
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The first term in (2.1) is the total kinetic energy of the gas and the second 
term is —7' multiplied by the entropy. The total number of particles is 
fixed, so we must have 


Vi (ce) 
N= 55) fk) kak, ee 


whichever form we choose for f(k). If we make F an absolute minimum 
by setting the variation SF equal to zero, subject only to the condition 
(2.2) then, as is to be expected, one finds that f(/) must be the usual 
Fermi—Dirac function. It is however, possible to obtain approximate 
equations (which remain thermodynamically consistent) even if we place 


Fig. 1 


Approximate forms for f(k) in place of the Fermi—Dirac function: (i) a>1; 
(ii) «<1. 


additional restrictions on the form of f(k) which enable us to evaluate 
(2.1), (2.2) and (1.3) without difficulty. This has already been demon- 
strated for the function shown in fig. 1 (Lidiard 1951 b); we adopt this 
function here too, i.e. for s=1/x<1 we take, 


f(k)=1, 0<k<(a—1)/B, 
=a—Bk, (a—1)/B<k<a/Bp, 
=O) ESC Vy. Pee eee 


while for x>1, 
f(k)=a—Bk, 0<k<a/B, 
=0, eS 8, me Sy OE rhode be ore iaede ame 
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Substitution of (2.3) and (2.4) into (2.2) gives 


Vi (a\8 3x Ake 
Vi fa\3 1 
=> Ag B ms > pl, Go OG. 6" oe ea Es (2.6) 
Similarly the evaluation of (2.1) with the elimination of («/8) by (2.5) 
and (2.6) yields io Sy okT 
Vers (~)— Den C(O) eee ee ee (2.7) 


where A(x) and B(x) are the following functions of x, in which 
P=(1—32/2+-22—23/4)-1, 
A (x)= P*8 [1—(1—x)*]/6xz, x<1l, 
= 27/3 42/3/3, a>, 
C(x)=P (x—x?+ 1123/36), Belt 


Do) 


Z) 2 ‘< x 
26% Ta? | 2a8 
a EP eEy 
Finally we fix the value of x by using the condition that F/ should be 
a minimum, i.e. df/dx=0. This gives the implicit equation 
Rep Ae DCL) ee eee (248) 
where the prime denotes differentiation with respect to x. Limiting 
forms of these equations have been given previously in the paper cited 
above and will not be repeated here. This paper also contains a com- 
parison of the approximate with the accurate functions computed by 
Stoner (1939), using the tables of FermiDirac integrals previously 
drawn up by McDougall and Stoner (1938). Unfortunately there are no 
tables available which could help us to evaluate (1.3) accurately. We shall 
proceed with the calculation of D(r) using (2.3) and (2.4) for f(k) with 
(«/8) and x given by (2.5), (2.6) and (2.7). The integrations are perfectly 
straightforward and abbreviating kyr by p we find that, 


—Saaay 


3 : : 
where, J (x, p)= pale) sin {p(1—x)P1/3}—sin {pP'/3} 


2p"! P-1 cos {p(1—a) P13} 


a9 Mle 


D(r) 


= oe Pt aic08 {pw | eee lee eee oe nen 20) 
3 : 2 

= xp | —sin {ati \—E queria, Waco (4x p}) | 3 
We lon as, ogra ae gh nt hag ee (TO) 


The parameter x is related to the temperature through equation (2.8). 
VD(r)/N as calculated from (2.9) and (2.10) is plotted as a function of p 
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in fig. 2 for five different temperatures, namely kT'/ey=0-0, 0-5, 1-0, 2-0 and — 
4-0. The oscillatory behaviour which one would expect from the presence ~ 
of the sine and cosine functions is too strongly damped to show up in a — 
: ; 
drawing of the scale of fig. 2. 
Fig. 2 


ie) 10 2.0 p 3.0 40 G0) 


The pair distribution function for a perfect Fermi—Dirac gas at several different 
temperatures. The numbers on the curves are the corresponding values 
of RT/«, and the variable p is kyr, where ky is the radius of the Fermi sphere 
at absolute zero. The dashed line for RT/e, equal to 4-0 has been cal- 
culated from the high temperature expansion given by Bhatnagar and 
Singwi. 

It is not difficult to find a simple limiting form for D(r) at very low 
temperatures. In this region ~=3kRT7'/e) and when ap <1 we obtain, 


De =! [1—9{ eee p—p cos aa 1) 


Except for the occurrence of 3RT'/e) in place of 7RT/e, this expression is 
identical with the low temperature formula derived by Bhatnagar and 
Singwi (1949, eqn. (19)) by means of the Sommerfeld lemma (note that our 
k,) is equal to their 27kp). 

§ 3. Discusston 


The replacement of 7kT'/eg by 3RT'/e) in formula (2.11) appears to be a 
general consequence at low temperatures of the use of the approximate 
function (2.3); for it is known to occur also with the formulae for the 
specific heat (Wohlfarth 1950) and electron spin paramagnetism (March 
and Donovan 1953). We may be sure therefore that formula (2.9) is 
accurate to a few per cent in the region 0<RT'/e,<0-3. Unfortunatel it 
is difficult to form an estimate of the general accuracy at high tone 
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since (2.10) does not go over into any simple form which could be compared. 
directly with (1.1). The following observations are however, relevant : 

(1) The accurate curve for RT/e,=4-0, calculated from eqn. (12) of the 
paper by Bhatnagar and Singwi, lies across that calculated from (2.10)— 
being greater for small values of p. 

(2) Nevertheless at higher temperatures the approximate and the 
accurate curves come into coincidence for small p. Thus expanding 
(2.10) we get 


Nf 1 (kT 101 /RT\2 
Dm=F|5(=) sa (—) ea ye) aoa) 
whereas (1.1) becomes 
NP 1 (kT 1 /kRT\2 
Be ee (ae \ete ee fees ; 
D(r)= ris(—)e (| p =: Be ee) 


(3) A? (x, p) has roughly the same shape as the accurate expression 
exp (—R7'p?/2c)) (except for the suppressed oscillations) and we may calcu- 
late the second moment of 4? with a view to finding an equivalent 
Gaussian form exp — Yp?. In this way we obtain 42 ~exp (—R7'p?2/1-76¢,) 
compared with exp (—RT'p?/2¢,). | 

We conclude from (2) and (3) that the accurate and approximate 
expressions are in reasonable agreement at high temperatures although 
our equations predict that 4? falls to zero too rapidly. The radius of the 
Fermi hole which (2.10) gives as 


pPr= 27/41/38 c, 


4a [ &, \1/2 kT 
= (=) for (=) >1, eee wet (3:3) 


seems to agree well with (1.1). 
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ABSTRACT 


An account is given of experimental work on the reaction between 
bromine and (1) thin films of silver sulphide, (2) composite thin films of 
silver sulphide with silver and gold. It has been established that, under 
these conditions, silver sulphide reacts with bromine to form silver and a 
sulphur bromide before silver or gold is attacked. Hickman postulated 
the existence of this reaction in 1927 to account for the increase in sensi- 
tivity which follows the digestion of photographic emulsions with sulphur 
sensitizers. On the basis of the experimental observations, a new mechan- 
ism is advanced for the formation of the latent image in modern high speed 
negative emulsions. 


§1. INTRODUCTION 


THE discovery by Sheppard (1925) and his co-workers that the sensitivity 
of washed silver halide emulsions could often be considerably enhanced by 
digesting them at temperatures up to 65°c with suitable sulphur com- 
pounds represented a great advance in photographie science. They 
deduced from their experiments that silver sulphide was formed during 
digestion and favoured the hypothesis that it was present on the surfaces 
of the crystals of silver halide in the form of discrete specks which provided 
centres for the concentration of silver atoms produced by the action of 
light (see Sheppard 1928 a). Gurney and Mott (1938) later described a 
two stage mechanism for the concentration of the silver atoms, which 
could be applied to Sheppard’s hypothesis if it were assumed that silver 
sulphide trapped electrons. 

An alternative explanation of the action of silver sulphide in increasing 
sensitivity was advanced by Hickman (1927) who suggested that it acted 
as a halogen acceptor by the liberation of metallic silver and the formation 
of sulphur monobromide. He was unable to obtain any experimental 
evidence to support this assumption, but wrote: “Even so, it is possible 
that, when the reaction takes place on the small scale governing latent 
image formation, the bromine liberated atom by atom by light in contact 
with an excess of silver sulphide may indeed liberate metal.” He 
demonstrated the feasibility of his ideas by making experiments with 

eee eee 
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pieces of sensitive photographic materials which had been visibly darkened 
by immersion in a dilute solution of sodium sulphide. Strips were 
exposed to light until a print-out image was produced and then immersed 
in a solution of sodium thiosulphate and potassium ferricyanide, which 
dissolves silver bromide and metallic silver but not. silver sulphide. 
The silver sulphide was found to have disappeared from the exposed parts. 
This would scarcely be expected if the only role of the silver sulphide was 
to concentrate photolytic silver atoms, and Hickman concluded that 
bromine liberated by light attacked the silver sulphide. This conclusion 
was not generally accepted. In the discussion following the presentation 
of Hickman’s paper to the Royal Photographic Society, Renwick remarked 
that “ Dr. Hickman could not expect anyone to accept his hypotheses, 
contrary as they were to many experimental facts, without fresh evidence 
in theirfavour”. Sheppard (1928 a) also rejected this interpretation of the 
sensitizing action of silver sulphide, although he admitted the importance 
of halogen removal to which Hickman had drawn attention. Sheppard 
expressed the opinion that the molecules of gelatine around the silver 
halide crystals would combine with the extremely small amount of halogen 
liberated during an exposure of the duration necessary to form a latent 
image and this viewpoint was also taken by Gurney and Mott. 

The purpose of this paper is first to describe experiments which were 
designed to test Hickman’s hypotheses and then to describe further 
experiments which appear to lead directly to a possible new mechanism 
for the formation of the latent image in modern high speed negative 
emulsions. We have observed that a thin film of silver sulphide is 
converted into a thin film of silver by treatment with bromine at a very 
low vapour pressure. The sulphur combines with bromine, apparently 
quantitatively, to form a sulphur bromide. The silver film which is 
formed is attacked by the sulphur bromide or by excess bromine at a 
slower rate and is then converted into silver bromide. This work lends 
direct support to one of the hypotheses advanced by Hickman. 

Prior to Sheppard’s work, Liippo-Cramer (1921, 1927) had proposed that 
the digestion of washed emulsions with gelatine resulted in the formation 
of particles of colloidal silver which were adsorbed to the surfaces of the 
silver halide crystals. That this was the cause of the increase in sensi- 
tivity during digestion came to be widely accepted, because it was reason- 
able to assume that the particles could increase in size during exposure 
and thus be converted into development centres. This conception, 
however, gradually faded into the background following Sheppard’s 
discovery and interpretation of the sensitizing action of sulphur com- 
pounds. Clark (1927) presented evidence in favour of composite sensi- 
tivity specks consisting of silver and silver sulphide, but his experimental 
technique was criticized by Sheppard (1928b). Sheppard (1928 a) 
nevertheless admitted that sensitivity specks might consist of silver and 
silver sulphide; at the same time, he insisted that their principal function 
was to concentrate silver atoms liberated during exposure. 
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Interest in chemical sensitization by silver has recently been revived by 
the publication of a series of papers by Mitchell (1948), Mueller (1949), 
Lowe, Jones and Roberts (1951), Hedges and Mitchell (1953) and Wood 
(1953). Of these, the one most relevant to the present work is that by 
Lowe and his co-workers, who showed that, with a positive emulsion of 
the type which they used, the effects of sulphur sensitization and silver 
sensitization were additive. The increased sensitivity may be attributed 
to the formation of layers adsorbed at the surfaces of the crystals of silver 
halide in which atoms of silver and molecules of silver sulphide are closely 
associated. We therefore studied the action of bromine on superimposed 
thin films of silver and silver sulphide and found that bromine at very low 
vapour pressures removed the sulphur from the silver sulphide with the 
formation of a sulphur bromide and left a thin film of silver of increased 
thickness. This film was itself attacked at a slower rate by sulphur 
bromide or by excess bromine. 

Modern negative emulsions owe their high speed to the sensitizing action 
of complex salts of monovalent gold such as ammonium or alkali aurous 
thiocyanates or thiosulphates which was discovered by Koslowsky (1951) 
in 1936 in the Agfa Laboratories at Wolfen, Germany (see also Chilton 
1946, Mueller 1949). Although Hedges and Mitchell (1953) have demon- 
strated that a crystal of silver bromide can be sensitized for the formation 
of a surface latent image by the deposition of a thin film of gold on its 
surface, it is most probable that composite layers consisting of adsorbed 
atoms of gold and molecules of silver sulphide are formed on the surfaces 
of the microcrystals during the digestion of the emulsion with these added 
sensitizers. For this reason, the action of bromine at very low vapour 
pressures on composite films of silver sulphide and gold was studied. 
It was found that the silver sulphide again reacted with bromine to form a 
sulphur bromide and that a composite metallic film of silver and gold 
remained which reacted much more slowly with the sulphur bromide or 
with excess bromine to form a mixture of silver bromide and aurous 
bromide. 

The possible application of these observations to the discussion of the 
formation of the latent image in high speed negative emulsions will be 
considered in a later section. The important conclusion is that the 
reaction between silver sulphide and bromine occurs in two stages, the 
first resulting in the formation of a sulphur bromide and metallic silver and 
the second in the formation of silver bromide. Fresh evidence has thus 
been obtained in favour of Hickman’s hypothesis that silver sulphide 
reacts with bromine during the exposure of sensitized silver halide 
crystals to liberate silver and form a sulphur bromide. 


§2. ExperimentaL Mrtuops 
The thin films of silver and gold which were used in all the experiments 
were deposited on the inner surfaces of highly evacuated thin walled 
spherical bulbs of Hysil glass, 10 em in diameter, by heating beads of the 
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metals supported on tungsten filaments near the centres. The beads were 
prepared and processed by the methods described by Mitchell and 
Mitchell (1951). About one third of the experiments were done with bulbs 
which had a single bead mounted on a pinch with two tungsten leads. 
Two beads were required for all the experiments with composite films. 
They were screened from each other with a disc of molybdenum and the 
two tungsten filaments supporting them were spot welded to four 
tungsten leads passing through a pinch. This pinch was sealed in so that 
the dise was at the centre of the bulb. The two beads were symmetrically — 
spaced on opposite sides of it so that the shadow of the disc cast by each 
_ bead covered approximately one third of the internal surface; this 
meant that each component film covered one third of the area of the bulb 
while the composite film covered the remaining third where the two films 
overlapped. With this arrangement, it was possible to make reliable 
estimates of the relative rates of many reactions. A side arm opposite the 
pinch was drawn down to a long fine capillary through which liquids could 
conveniently be admitted to the evacuated bulb. The bulbs were 
provided with two other side arms and one of these was sealed to the high 
vacuum system through a constriction. This carried several glass bubble 
seals, through which the bulb could be sealed for re-pumping after the 
admission of gases, and also a long plain tube and a U tube either or both 
of which could be immersed in liquid air to condense volatile reaction 
products. Tubes with breakable tips containing bromine, sulphur mono- 
bromide, hydrogen sulphide, air and nitrogen were sealed to the opposite 
side arm as required. Short lengths of iron rod or ball bearings, sheathed 
in glass, were used for breaking the tips and the glass bubbles. 

Hydrogen sulphide was prepared in the usual way and purified by 
fractional distillation. Nitrogen was obtained from a cylinder of the 
oxygen-free gas. Analar bromine supplied by B.D.H. was used through- 
out the experimental work. The tubes with breakable tips were filled 
with small quantities of bromine, using an all glass apparatus which was 
evacuated while the side arm containing the supply of bromine was 
immersed in liquid air and then sealed from the pumping system. Sulphur 
bromide was prepared by reacting an excess of vacuum sublimed Analar 
sulphur with Analar bromine in a previously evacuated all glass apparatus. 
There was no significant evolution of heat when the reactants were mixed 
and, after a period of refluxing, the product was fractionated by immersing 
the side arms in turn in liquid air while the bulb containing the product was 
at room temperature. When a fraction had been collected, this bulb was 
also immersed in liquid air to reduce the vapour pressure to a negligible 
value before the side arm was sealed off. The sulphur bromide which was 
used appeared to be a product of definite composition. Bromine could not 
be separated from it by slow vacuum fractional distillation at low temper- 
atures and no residue of sulphur remained after the complete distillation of 
a sample between room temperature and the temperature of liquid air. 
Every possible precaution was taken to ensure that it did not contain free 


bromine. 
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The completed bulbs were pumped for several hours and repeatedly 
torched during this period. The beads were meanwhile heated to a 
temperature just below that at which the metal would have begun to 
evaporate at a significant rate. When the pressure had fallen to 
10-7 mm Hg, a bead was heated until the metal began to evaporate slowly 
from it. The thickness of a film was judged visually from its colour by 
transmitted light and deposition was interrupted when the required 
thickness was reached. If composite films were being prepared, the 
second bead was heated and the corresponding film deposited before the 
bulb was sealed from the pumping system. 

After sealing off, the films could be exposed to hydrogen sulphide by 
breaking the tip of the tube containing it. The bulb was heated and the 
excess hydrogen sulphide was then condensed as far as possible by 
immersing the tube in liquid air before it was sealed off. The bulb itself 
was then sealed to the pumping system through one of the glass bubbles 
which was broken when the pressure in the system had fallen to 10~-* mm 
Hg and pumping continued until the pressure had fallen to 10-7 mm Hg. 
The second component of a composite film could then be deposited and the 
bulb again sealed from the pumping set. 

Films of silver sulphide, silver and gold, and composite films of silver 
sulphide and silver, silver sulphide and gold, and of silver and gold were 
prepared by these methods and their reactions with bromine and sulphur 
bromide both under good vacuum conditions and in the presence of dry and 
moist air and nitrogen at pressures up to one atmosphere were investigated. 
In all the experiments the tube containing the bromine or sulphur bromide 
was immersed in liquid air before the tip was broken. It was then 
allowed to warm up very slowly in the space above the surface of the 
liquid air in a deep Dewar vessel and, as soon as the appearance of the 
film began to change, was plunged back into liquid air so that the vapour 
pressure fell rapidly. The tube was often sealed from the spherical bulb 
as soon as the conversion of silver sulphide to silver appeared to be 
complete. The sulphur bromide could then be condensed in the plain side 
arm cooled in liquid air and fractionated at a low temperature. If the 
bulb had previously been filled with air or nitrogen, it could be sealed to 
the vacuum system at this stage through a glass bubble seal, pumped out 
through the U tube which was cooled in liquid air, and then sealed off again. 
The volatile reaction products could then be condensed either in the side 
arm or in the U tube. 


§3. EXPERIMENTAL OBSERVATIONS 


The colour transmitted by a thin film of silver passes through the 
sequence, pale yellow, yellowish brown, reddish violet, bluish violet, and 
blue fading into bluish grey as the mass of silver deposited per unit area of 
the substrate is increased. The bluish violet films are the first in the 
sequence to give a white metallic reflection. The bulk of the experimental 
work described in this paper was done with unannealed reddish violet or 
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bluish violet films of silver which could be deposited reproducibly and 
showed the most pronounced colour changes when they reacted with other 
substances. Bromine attacked thin films of silver rapidly at room temper- 
ature to form clear films of silver bromide which recrystallized on heating 
becoming pale yellow and opalescent. Reddish or bluish violet films of 
silver were converted instantly into yellowish brown films of silver 
sulphide by exposure to hydrogen sulphide. Films of silver sulphide 
prepared in this way were stable and could be heated in vacuum to 
300°C without any condensation of sulphur in a side arm cooled in liquid 
air. Thin films of silver were attacked by sulphur monobromide with the 
formation of silver bromide and sulphur. The sulphur could be condensed 
in a side arm cooled in liquid air by warming the bulb. 

The films of silver sulphide reacted with bromine at a very low vapour 
pressure to form films of silver and a sulphur bromide. The films of silver 
produced by this reaction had the same colour as those which had been 
converted into silver sulphide for the experiment. The sulphur bromide 
could be partially condensed in a side arm cooled in liquid air or else 
completely dissolved out by admitting dry benzene, chloroform, or carbon 
tetrachloride to the bulb through the capillary. The solution of silver 
nitrate produced by dissolving the silver film remaining in the latter case 
in hot nitric acid gave no precipitate with barium nitrate, showing that the 
sulphur had been removed within the limits of the test. A precipitate of 
barium sulphate was obtained when a control experiment was run with a 
film of silver sulphide. If the tube containing the bromine was not im- 
mersed in liquid air and sealed off, immediately the transformation of the 
film of silver sulphide to silver appeared to be complete, the film of silver 
was slowly attacked by the sulphur bromide and by bromine and converted 
to silver bromide. In many experiments dry or moist air or nitrogen was 
admitted to the bulb at pressures up to one atmosphere before the tip of 
the bromine tube was broken. The presence of air or nitrogen did not 
appear to influence the primary phenomena in any way. 

The reaction between a thin film of silver sulphide and sulphur mono- 
bromide was next studied. Sulphur monobromide appeared to react 
catalytically with silver sulphide to form silver and sulphur at a rate which 
was not very different from the rate of reaction between silver sulphide and 
bromine. The same result was obtained when the bulb was filled with 
dry air or nitrogen before the admission of the sulphur monobromide. 
An excess of sulphur monobromide reacted with the silver to form silver 
bromide and sulphur. 

Following these experiments with thin films of silver and silver sulphide, 
a composite film of silver and silver sulphide was prepared as described in 
the previous section. A reddish violet thin film of silver was first deposited 
over two thirds of the inner surface of the bulb and converted to a yellow 
brown film of silver sulphide by exposure to hydrogen sulphide. The bulb 
was then exhausted and another reddish violet film of silver was deposited 
from the second bead so that it overlapped half the area of the silver 
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sulphide film. The bulb was evacuated or else filled with air or nitrogen in 
the different experiments of this series. The introduction of bromine at a 
very low vapour pressure into the system converted the yellow brown film 
of silver sulphide back to a reddish violet film of silver before it had 
affected the reddish violet film of silver on the other side of the bulb in 
any way. The composite film of silver and silver sulphide was changed 
from a dull greyish violet to a clear blue violet by exposure to bromine 
corresponding to the formation of a thicker film of silver. When more 
bromine was admitted, it was observed that the silver film derived from 
the silver sulphide was attacked at the same rate as the other silver film 
with the formation of silver bromide. The rate at which the thicker silver 
film derived from the composite film was attacked was not sensibly 
different. These experiments were repeated many times. 

We turned now to the investigation of the properties of thin films of 
gold. With increasing mass of gold per unit area, the colour transmitted 
by the films passed through the sequence, reddish violet, bluish violet, blue, 
bluish green with a white reflection, and bluish green with a golden yellow 
reflection. When such films were heated to temperatures above 150°c, 
they were transformed into pink films which showed the same colour by 
transmitted as by reflected light. The latter are known from other work 
to have an infinite electrical resistance and to consist of discrete particles 
of gold. These thin films, whether previously converted to pink films by 
heating or not, reacted with bromine to form clear pale green films of 
aurous bromide and then yellowish brown opalescent films of auric 
bromide. If the tube containing the bromine was immersed in liquid air 
and the bulb was heated to a temperature in the neighbourhood of 120°c, 
the film of auric bromide broke down forming first a green film of aurous 
bromide and then a pink film of gold. Coloured thin films of gold of all 
thicknesses were changed irreversibly into pink films by reaction with 
bromine followed by thermal decomposition of the gold bromides. No 
conditions could be found under which thin films of gold would react with 
hydrogen sulphide, which is in striking contrast to the immediate reaction 
which occurs with silver films. The reaction between bromine and a 
reddish violet film of silver, a bluish violet film of gold and a deep violet 
composite film of silver and gold was next studied to determine whether, 
under identical conditions, there was any significant difference in the rates 
of attack. It was found that all three films were attacked at the same rate 
by bromine with the formation of silver bromide and aurous bromide. 
Experiments with sulphur monobromide instead of bromine led to a 
similar conclusion. The same observations were made when the bulb was 
filled with air or nitrogen before the admission of the bromine. 

The reaction of composite films of silver and gold with hydrogen sulphide 
was studied in the next series of experiments. In half of them, the gold 
was superimposed on the silver and in the other half, the silver was 
superimposed on the gold. The silver films were bluish violet by trans- 
mitted light and the gold films were bluish green and showed a white 
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reflection. The composite films were a deep bluish violet. Treatment of 
either type of film with hydrogen sulphide produced exactly the same 
result : the silver was converted to silver sulphide and the gold remained 
unaffected. The whole of the area covered by the gold film became blue 
green by transmitted light and the boundary between the gold film and the 
composite film almost vanished. The bulb was then evacuated and 
sealed off or else filled with air or nitrogen. Exposure toa very low vapour 
pressure of bromine caused the film of silver sulphide and the film in the 
composite area to return to the colours which they had before treatment 
with hydrogen sulphide and sulphur bromide could be condensed in the 
side arm cooled in liquid air. The gold film was unaffected at this stage. 
Further reaction with bromine transformed the silver into silver bromide 
and the gold into aurous bromide and finally, with excess bromine, into 
auric bromide. The same changes were observed whether the silver or the 
gold of the composite film was deposited first and whether the bulb was 
evacuated or filled with air or nitrogen before the admission of the bromine. 
In later experiments, sulphur bromide was admitted to the bulb instead of 
bromine. The same changes were observed but auric bromide was. not 
formed. For the last group of experiments, a composite film of silver 
sulphide and gold was prepared in a different way. A thin film of silver 
was deposited over two thirds of the surface of the bulb and converted to 
silver sulphide. A thin film of gold was then deposited to cover one half 
the film of silver sulphide. The appearance of the bulbs was the same as 
in the previous series and when bromine and sulphur bromide were 
admitted the same changes were observed. 


$4, DiscUSSION OF THE EXPERIMENTAL RESULTS 


In recent years, the properties of thin films of silver and other metals 
have been extensively investigated at Bristol (see, for example, Holloway 
1951, Dorling 1953). The accumulated evidence has led to the working 
hypothesis that the unannealed conducting films are built up from 
crystalline granules consisting of small groups of atoms. The pink 
non-conducting films mentioned above have larger discrete groups of 
atoms. We have therefore good grounds for assuming that, in the thin 
films of silver and gold which we have used for our experiments, surface 
atoms of metal are adsorbed to the glass substrate in the same way as 
atoms of silver and gold are adsorbed to the surfaces of the crystals of 
silver halide in a sensitized photographic emulsion. For this reason, 
results obtained with thin films of silver and gold may be directly 
applicable to the discussion of the formation of the latent image in the 
photographic system. 

A granular film of silver is most probably converted into a granular 
film of silver sulphide by treatment with hydrogen sulphide. In striking 
contrast to the immediate reaction which occurs with a thin film of silver, 
a thin film of gold does not react with hydrogen sulphide. Particles of 
gold may also be stable in a photographic emulsion under conditions 


216 J. H. Burrow and J. W. Mitchell on Experiments with 


where particles of silver would be converted into silver sulphide, as, for 
example, when sodium thiosulphate is added to oa emulsion at a Pr 
greater than 7-5 and at a temperature of about 60°C. The experiments 
in which composite films of silver and gold were treated with hydrogen 
sulphide and found to react in the same way whether the silver or the 
gold was deposited first provide further evidence for the existence of 
a granular structure. If the films were continuous, it would be anticipated 
that an overlying gold film would protect a silver film from attack by 
hydrogen sulphide. 

The observations on the reactions of films of silver sulphide, and of 
composite films of silver sulphide and silver, and silver sulphide and gold 
with bromine, are all consistent with the assumption that they have a 
granular structure. Silver sulphide whether present by itself or in 
association with silver or gold is converted into silver by reaction with 
bromine at very low vapour pressure. The sulphur of the silver sulphide 
reacts to form a sulphur bromide. This reaction occurs under conditions 
where particles of silver or of gold are not attacked. The experimental 
conditions approximate very closely to those which are likely to hold 
at the interface between a crystal of silver bromide and the gelatine 
during a normal photographic exposure. We may therefore conclude 
that the increase in sensitivity of a photographic emulsion which results 
from sulphur sensitization is at least partly due to the formation of silver 
sulphide which can liberate silver when it reacts with bromine as was 
suggested by Hickman in 1927. The details of his mechanism for the 
formation of the latent image require revision, however, in the light of 
modern knowledge and this will be attempted in the next section. 

There is no experimental evidence that sulphur monobromide is the 
primary product of the reaction between silver sulphide and bromine. 
The primary product could equally well be sulphur dibromide which 
could subsequently condense with itself to form sulphur monobromide 
and bromine. The reaction between sulphur monobromide and silver 
sulphide is of considerable interest. The sulphur monobromide was 
most carefully freed from excess bromine and could be distilled unchanged 
in vacuum at room temperature. It is, nevertheless, possible that it 
suffered heterogeneous catalytic thermal decomposition on the surfaces 
of the films in the spherical bulb into sulphur and bromine. Bromine 
would then be made available for reaction with the films and sulphur 
would be deposited on the surface. The results seem to suggest that a 
thin film of silver sulphide might be converted into silver and sulphur 
by catalytic reaction involving less than the equivalent amount of bromine 
or sulphur monobromide. This matter requires further detailed investi- 
gation. Sulphur monobromide reacts more rapidly with silver sulphide 
than with particles of silver or gold of approximately the same size under 
the same conditions. 

All the experiments involving bromine were carried out using both 
bulbs which had been evacuated before the admission of bromine and 
bulbs which had been filled with air or nitrogen at pressures up to one 
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atmosphere. No difference in the primary reactions due to the presence 
of these gases was ever observed. The primary reactions were also not 
affected by the presence of traces of water vapour. The experimental 
results may therefore be applied with some confidence to the photographic 
system. When finely divided sulphur is produced in a reaction, it is 
likely that it will react with oxygen in the presence of moist air to form 
sulphur dioxide and, ultimately, sulphuric acid but this was not estab- 
lished in the present series of experiments. 


§5. THE APPLICATION OF THE RESULTS TO THE PHOTOGRAPHIC PROCESS 


The material of this section is inevitably speculative and has therefore 
been separated from the account of the experimental results. Any 
satisfactory theory of the formation of the latent image must explain 
how the most sensitive crystals of a modern high speed negative emulsion 
become developable after the absorption of very few quanta. With 
the results of the present series of experiments, it is possible to account 
for this. It is also possible to explain the function of the two main 
components of the sensitizing layer, silver sulphide and gold, and the 
effectiveness of gold sensitization in reducing high intensity reciprocity 
failure (see, Mueller 1949). 

For purposes of discussion, we shall assume that the sensitizing 
substance forms a layer, perhaps a monolayer, of molecules of silver 
sulphide with atoms of gold in addition, principally at the surface of 
the sensitivity centre (see Mitchell 1953) during the digestion of the 
washed emulsion with appropriate amounts of sodium thiosulphate and 
aurous thiocyanate. In the first stage of the formation of the latent 
image in a microcrystal, the quanta may be absorbed (1) by the lattice 
of the silver bromide, (2) by bromide ions at the surface of the crystal 
or (3) by atoms or molecules of the products of chemical sensitization. 
(2) and (3) are responsible for the inherent long wave or red sensitivity 
of the emulsion. Positive holes and electrons, both of which can diffuse 
through the crystal are liberated in the first case (1). The absorption of 
a quantum by a surface bromide ion liberates an electron and leaves a 
bromine atom and an excess silver ion at the site of absorption. An 
electron may be ejected from an atom of gold or of silver by the absorption 
of a quantum leaving an aurous or a silver ion behind. By the absorption 
of a quantum of the appropriate energy an electron may also be ejected from 
a molecule of silver sulphide. The positively charged molecule may then 
dissociate into a silver ion, a silver atom and a sulphur atom. 

We shall first consider the case where the absorption of a quantum of 
relatively short wavelength liberates a positive hole and an. electron. 
Mitchell (1953) has recently pointed out that the lifetime of the hole is 
probably less than that of the electron. The hole is most likely to be 
trapped either by a bromide ion at the surface of the senetyy Gy centre 
or by an atom or molecule of a product of chemical sensitization. The 
trapping of a hole by a bromide ion produces a bromine atom and creates 
a local positive charge in the form of an excess silver ion. The electron 
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will be held in the neighbourhood by electrostatic attraction but cannot 
immediately combine with the excess silver ion (see Mitchell 1953) ; 
neither can it combine with the bromine atom once the vibrational 
excitation energy of that atom has been dissipated. The bromine atom 
may, on the other hand, combine with a molecule of silver sulphide to 
release a silver atom and form a molecule of silver sulpho-bromide, 
AgSBr, which may subsequently break down into a molecule of silver 
bromide and an atom of sulphur. The silver atom may now combine 
with a gold atom to form a pair of atoms which may in turn associate 
themselves with the excess silver ion to form a positively charged triple 
ageregate. The positive charge of this aggregate can then be neutralized 
by the electron so that a neutral aggregate consisting of two silver atoms 
and one gold atom can result from the absorption of a single quantum. 
This aggregate may then combine with a silver ion from the surface of 
the silver halide crystal to form a positively charged aggregate of four 
atoms, probably in a tetrahedral arrangement, which could trap an 
electron directly and by a succession of similar processes continue the 
building up of the latent image. The positively charged aggregate of 
four atoms may also satisfy the minimum requirements for a development 
centre. The mechanism outlined here does not involve any long range 
motion of silver ions for the neutralization of space charges and may 
therefore be of the greatest significance under high intensity conditions 
of exposure. 

The same sequence of events would follow the absorption of a quantum 
of lower energy by a bromide ion which was adsorbed or occupied a 
kink site at the surface. Ifa quantum of suitable wavelength is absorbed 
by a molecule of silver sulphide, an electron may be ejected leaving a 
positively charged molecule which may break down into a silver ion, 
a silver atom and a sulphur atom. Here, the observations of Rigollet 
(1895) are of interest; he found that electrons could be released from a 
thin film of silver sulphide in contact with silver by the absorption of 
radiations with wavelengths between 1:32 4 and 0:45 x. Should electrons 
be ejected from molecules of silver sulphide, the mechanisms already 
discussed would apply : the silver ion and the silver atom could associate 
with the electron and a gold atom to form a neutral aggregate of three 
atoms which could then combine with a further silver ion to form a 
positively charged aggregate of four atoms. This mechanism would 
also operate if the positive hole resulting from the absorption of a 
quantum by the silver bromide lattice were trapped by a molecule of 
silver sulphide. The ejection of an electron from an atom of silver 
or of gold by the absorption of a quantum could only result in the 
formation of a pair of metallic atoms. This would also apply if the 
positive hole were trapped by an atom of silver or gold, either in an 
emulsion which had been sensitized with silver sulphide and silver or 
gold, or in an emulsion which had been sensitized only with silver or gold. 
In dye sensitized emulsions, the quanta are absorbed by the dyestuff, 
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which is adsorbed in the form of a monolayer over the whole surface 
of the crystal and the excitation energy is transferred, with greatest 
probability, either to bromide ions occupying surface sites of low positive 
electrostatic potential (for example, kink sites) at the sensitivity centre, 
or else to atoms or molecules of the products of chemical sensitization 
also at the sensitivity centre; electrons are then ejected. The exciton 
must move about in the adsorbed layer of dyestuff until this happens, 
unless the excitation energy is degraded to thermal energy. The subsequent 
processes are the same as if the electron had been ejected from the bromide 
ion or from the sensitizing atom or molecule by the direct absorption of 
a quantum and have already been discussed. 

We shall next consider the formation of the latent image in crystals 
of silver bromide, which have been subjected to both sulphur and 
reduction sensitization but have not been gold sensitized. Lowe, Jones 
and Roberts (1951) have shown that the sensitivity of certain relatively 
fine grained emulsions which have been sensitized with sulphur compounds 
can be considerably enhanced by subsequent sensitization with suitable 
reducing agents. We shall assume that, in this case, the sensitizing layer 
consists of a close association of molecules of silver sulphide and atoms 
of silver. The sequence of events already described for a gold sensitized 
emulsion would lead to the formation of a pair of silver atoms, followed 
by a positively charged aggregate of three silver atoms which could then 
combine with the electron to form a neutral triple aggregate. 

We are naturally led to inquire why the treatment of such an emulsion 
with a gold salt which will only replace the silver atoms by gold atoms 
produces a very considerable increase in sensitivity. It is important 
to remember here that the same increase in sensitivity is obtained 
whether the emulsion is treated before or after exposure with a salt 
of monovalent gold (see Mueller 1949, James 1948, James, Vanselow 
and Quirk 1949), so that at least part of the explanation must be concerned 
with the mechanism of development. As a result of their work with large 
single crystals of silver bromide, Keith and Mitchell (1953) suggested 
that chemical development results from a combination of two processes : 
first, the transference of silver ions from the surface of the silver bromide 
erystal to the metallic nucleus produced during exposure and, secondly, 
the transference of electrons from molecules of the developer to the 
positively charged nucleus. The same mechanism applies to physical 
development, but here the silver ions are derived from the solution. 

There seems to be no reason why a group of three silver atoms should 
not be formed during the exposure of an emulsion which has been both 
sulphur and reduction sensitized following the absorption of a single 
quantum by exactly the same mechanisms as in an emulsion which has 
been sulphur and gold sensitized. In seeking an explanation of the 
fact that a triple nucleus which includes a gold atom is more effective 
in initiating development than a nucleus consisting entirely of silver 
atoms, we have therefore to consider the next step in the mechanism 
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in which the nucleus combines with a silver ion from the surface of the 
crystal to form a positively charged aggregate of four atoms. The 
compensating negative charge will reside on the surface probably in the 
form of two adjacent kink sites both occupied by bromide ions. The 
question is whether more energy is liberated when a silver ion combines 
with a triple nucleus containing gold atoms than when the nucleus 
consists solely of silver atoms. When the aggregates are positively 
charged, the energy of van der Waals interaction must make an important 
contribution to the cohesive energy. This component will be greater 
when the aggregate contains gold atoms than when it consists only of 
silver atoms because the polarizability of an aurous ion is greater than 
that of a silver ion. We may therefore conclude that an aggregate 
including gold atoms will have a greater tendency to become positively 
charged at a smaller size than an aggregate of silver atoms. This could 
account for the enhanced developability of the latent image which is 
formed under conditions of minimum high intensity exposure in gold 
sensitized emulsions compared with reduction sensitized emulsions. It 
could also account for the post-exposure latensifying action of gold salts 
on the latent image formed in emulsions which have not been gold 
sensitized. 

We have finally to consider the formation of the latent image in 
emulsions which have been sensitized with a sulphur compound under 
conditions which give a minimum of reduction sensitization. We shall 
assume that the sensitizing material consists only of molecules of silver 
sulphide adsorbed at the surfaces of the crystals. The sequence of 
changes which we have already discussed will now lead to the formation 
of a pair of silver atoms which will be much less likely to combine with 
a silver ion from the surface than a triple aggregate. This seems to 
account for the lower sensitivity of emulsions which have been sensitized 
only with sulphur compounds. 

We have assumed that isolated atoms of silver and gold are produced 
mainly at the surface of the sensitivity centre during reduction and gold 
sensitization. There is no experimental justification for this and pairs 
or larger aggregates might equally well be formed. Although a pair of silver 
atoms might dissociate at room temperature, a pair of gold atoms should 
form a stable unit (see Mitchell 1953). The absorption of a single quantum 
by a microcrystal of silver halide which had a distribution of molecules 
of silver sulphide and pairs of atoms of gold at the sensitivity centre 
could lead to the formation, under the most favourable circumstances, 
of a neutral group of four metallic atoms. We have so far restricted 
the discussion to the formation of the latent image under conditions of 
minimum exposure and it will be convenient at this point to summarize 
the conclusions. After the absorption of a single quantum by a micro- 
crystal of a non-sensitized silver halide emulsion, not more than one 
silver atom can separate at the surface. When molecules of silver 
sulphide are present, the absorption of a quantum can give rise to a 
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pair of silver atoms. If single atoms of silver or gold are associated 
with molecules of silver sulphide in the sensitizing layer, groups of 
three metallic atoms can result from the absorption of a single quantum. 
Such triple groups can then combine with silver ions from the surface 
of the crystal during or after exposure or when the crystal is immersed 
in a developer to form positively charged groups of four metallic atoms 
which can function as development centres. Positively charged groups 
containing gold atoms should be more stable than groups consisting 
only of silver atoms. These considerations provide a reasonable 
explanation of the progressive increases in sensitivity which follow the 
application of the different types of chemical sensitization. If the 
mechanisms proposed are indeed applicable to the photographic system, 
it should not be possible to obtain the same increases in sensitivity by 
digesting a pure silver iodide emulsion with a sulphur compound because 
iodine does not combine with sulphur to form a stable sulphur iodide. 
The adsorbed silver sulphide molecules would therefore have either to 
absorb quanta which would cause the ejection of electrons or else to trap 
positive holes before their silver atoms could contribute to the formation 
of the latent image. It might also be predicted that it would not be 
possible to optically sensitize a sulphur sensitized silver iodide emulsion 
as efficiently as a pure silver bromide emulsion. The excitation energy 
of the dyestuff would have to be transferred directly to molecules of 
silver sulphide before the sensitization could be effective at the longer 
wavelength. For the same reasons, the efficiency with which silver 
-bromo-iodide emulsions can be sensitized for longer wavelengths by a 
combination of sulphur sensitization and dye sensitization should be 
decreased as the concentration of iodide increased because an ever 
increasing proportion of kink sites would be occupied by iodide ions. 

It must be emphasized that the speculations of this section apply to 
only one aspect of a very complex problem. They are concerned with the 
detailed behaviour of the products of chemical sensitization during 
exposure and are supplementary to one of the author’s recent speculations 
on photographic sensitivity (Mitchell 1953). These products are most 
likely to be concentrated at the surfaces of the sensitivity centre which 
is, according to Mitchell, the region of maximum localized imperfection 
in the crystal. In a tabular crystal, the sensitivity centre will be the 
part of the crystal in the immediate neighbourhood of the inevitably 
imperfect nucleus from which the crystal grew and will consequently 
usually be located near the centre of the crystal. Now, if a molecule 
of bromine were formed there as a result of a rapid succession of events 
during a normal exposure, it would be able to react with a molecule of 
silver sulphide to form a pair of silver atoms and a molecule of sulphur 
dibromide. After dissociation, this molecule could then react with a 
further molecule of silver sulphide and the process could be continued 
so that under favourable circumstances the absorption of at least two 
quanta might cause the liberation of silver atoms from more than two 
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molecules of silver sulphide. This type of reaction would probably only 
occur at the surface of the sensitivity centre where bromine atoms are 
most likely to be formed by the trapping of holes and where there is 
the greatest surface concentration of products of chemical sensitization. 
The silver atoms together with gold atoms, silver ions and electrons 
would then condense to form development nuclei. It is perfectly clear 
that the formation of the latent image under average conditions of 
exposure in a fully sensitized modern high speed silver bromo-iodide 
negative emulsion is likely to be an extremely complicated process. 
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In the course of experiments on fatigue, we have made some observations 
which seem pertinent to the question of the origin of the fatigue crack. 

The experiments were done on pure copper fully annealed, and tested 
at constant stress in push-pull at about 1000 cycles per second. The 
diameters of the cylindrical specimens were 3-8 mm, and the grain size 
was about 0-!mm. The specimens were carefully prepared so as to be 
free from a cold-worked surface layer, and were finally electro-polished. 
When tested at stresses which gave a fatigue life of the order of 
10%10’ cycles, the surface showed, at a very early stage of the test, a 
number of slip bands. These became more numerous, and some of them 
became more intense, as the test proceeded. 

It was noticed that the eventual fatigue crack frequently appeared to 
have originated at such an intense slip band. The course of the crack 
in both directions from this centre point was marked by an abnormally 
high density of slip bands, while the centre itself was comparatively 
clean. Similar features were observed on other, smaller cracks, when 
these were present. 

The detection of a small crack under the microscope is made more 
difficult by the high general density of slip markings. Accordingly the 
expedient was adopted of electro-polishing the specimen some time before 
failure was expected. This removed the slip bands and revealed the 
presence of markings which we believe to be incipient cracks. Figure 1 
(Pl. 2, figs. 1-7), for example, shows a photograph taken after 
2-8 10% cycles. This crack spread and eventually caused fracture after 

3-65 x 108 cycles. 
: ee to Sasi: the point at which such cracks could first 
be detected, a specimen was electro-polished at intervals throughout a 
test, and photographs were taken each time of a number of fields, 
Figure 2 shows the appearance of one area after 2:7 10° cycles ; this 
was on the same specimen which failed, elsewhere, after 3-65 x 108 cycles. 
Figure 3 is the same field after electro-polishing to remove a layer about 
2 microns thick from the surface. It will be seen that most of the surface 
markings have been removed, but a few of the more intense bands remain. 
A second electro-polish, removing about the same amount as before, left 


Q2 


224 Correspondence 


the appearance substantially unchanged. Further periods of stressing 
in the fatigue machine, followed by electro-polishing, showed the same 
marks persisting, and new ones of similar character appearing. Some of 
these persistent slip bands began to spread into neighbouring grains 
(see fig. 1) and eventually gave rise to unmistakable fatigue cracks. One 
of these would grow in length more rapidly than the others, perhaps 
joining up with another in the process, and eventually giving rise to the 
fatigue fracture. 

One fatigue test was stopped when the above technique revealed about 
a dozen such incipient cracks, ranging in length from 3 to 10 grain 
diameters, together with a considerable number of the persistent slip 
band markings only one grain long. This specimen was then given a 
5%, slow tensile extension. All the ‘ cracks ’ opened out in the manner 
shown in figs. 4 and 5, but none of them showed any tendency to increase 
in length. The marks confined to a single grain were difficult to locate 
on the deformed surface, but it can safely be said that their tendency to 
open out in the above manner was very much less, and perhaps not 
present at all. 

At what stage a persistent slip band can be called a crack is perhaps 
rather a question of nomenclature. The above observations show the 
presence of persistent slip band markings after less than 10% of the 
fatigue life, and of the indubitable cracks into which they develop at 
about 60°, of the life. The markings cannot simply be large steps on 
the surface, since such steps, which are produced in abundance by 
moderate unidirectional strains (say 10°), were removed by the electro- 
polishing treatment used. Again, it was found that, if a specimen 
were subjected to alternating stress for, say, 10°% of its estimated 
fatigue life and then polished to show the persistent markings, and finally 
vacuum annealed at 600°¢ for 1 hour and re-polished, the same markings 
were still present. These results suggest a possible explanation of the 
observation of Sinclair and Dolan (1952) that a (brass) specimen can be 
annealed repeatedly during a fatigue test without any significant effect 
on the fatigue life. This. conclusion is supported qualitatively by a few 
results obtained by us on the copper specimens used in the above work. 

Preliminary observations on single crystals suggest a similar pattern 
of behaviour. A broad slip band, visible at, say, 10°% of the life, will 
be removed by electro-polishing except for a number of localized regions 
(figs. 6 and 7). One or more of these, on further stressing, will eventually 
appear as an obvious crack. The crack, although thus starting in the 
slip band, will later deviate from it, to pursue its more usual course 
approximately in the plane normal to the tensile axis. 

The observations are being continued, and a fuller account will be 
published in due course. 
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The Quantum Yield of the F+Z 1 Conversion in KCl with Sr Impurity 


By P. Camaenr,* G. CuraroTti, F. G. Fumr* and L. Grunorro 
Istituto di Fisica, Universita di Pavia, Italy 


[Received December 15, 1953] 


To study the interaction as electron traps of the Schottky defects and of 
the impurities in the alkali halides, we have investigated in detail the 
conversion of the F band to the simplest colour-centre band connected 
with the presence of positive divalent impurities, the Z, band (Pick 1939, 
Seitz 1951). KCl erystals with a small quantity of Sr++ added as SrCl, 
(5-5 x 10~° mole fraction) were coloured either additively in Na vapours 
or by irradiation with x-rays at room temperature and then exposed to 
the action of the 54614 line of Hg, which falls in the spectral region 
common to the F and Z, bands. 

In additively coloured crystals we found that practically only Z, band 
is formed while the F band decreases. Therefore we assumed the 
equality between the numbers per cm® of Z, centres formed and of 
F centres destroyed at the end of the conversion in a given crystal, and 
determined in this way the oscillator strength for the Z, band using 
Smakula’s formula (Seitz 1940). The value of 0-84} that we obtain is 
very close to the oscillator strength for the F band (Kleinschrod 1936), 
as it is reasonable it should be (Seitz 1951). 

Figure 1 gives the numbers of F and Z, centres per cm? as functions 
of the number of quanta absorbed in an additively coloured crystal. The 
numbers of centres are determined directly from the height and width 
of the F and Z, bands using the oscillator strengths. The quantum yield 
for the destruction of F centres (or the formation of Z, centres) at the 
beginning of the conversion is about 0-15 in our crystals, and goes 
gradually to zero during the conversion as the Z, band reaches a 
saturation value. 

The initial absence of the Z, band can be justified through the thermal 
instability of the Z, centre (Pick 1939). To explain the observed 
saturation of the Z, band, which is practically stable under Z, light, 
we suppose that at saturation all the Sr++ which are not associated to a 
positive-ion vacancy (free Sr++) have formed a Z, centre by capturing 
an electron and that the associated Sr++ do not act as effective electron 
traps at room temperature (Harten 1950){ : the number of unassociated 
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ions should not change rapidly at room temperature owing to the values 
of the association energy and of the jump frequency for free positive-ion 
vacancies. To investigate the validity of our interpretation we have 
exploited some of its consequences. 

If we assume that the only efficient electron traps are F and Z, traps 
and that each photon frees an electron from an F centre (Mott and 
Gurney 1938), the quantum yield for the conversion can be written 


Fig. 1 
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The numbers of F and Z, centres/em* as functions of the number of absorbed 
quanta/cm? (Hg 54614) in a erystal of KCl containing 9 x 1017 Sr++/cem3 
coloured additively. 


(Ntt_N, )o 
n= Z,/~Z, i f ; . (1) 
(N,—N y)opt+ (Nt —_N,z )oz, 


where N++ is the number of free Sr++/em? and thus the number of 
Z, centres/cm? at saturation: N,, N,, Nz, are respectively the numbers 
of negative-ion vacancies, F centres and Z, centres per unit volume ; 
o, and oy are the capture cross sections for electrons of a negative-ion 
vacancy and of a free divalent ion. The behaviour of 7 during the 
conversion as a function of Nz is found to be nearly linear and this 
indicates that o,/oz, is of the order of unity. If one takes op/oz,=1, 


one obtains N,~ 1-4 10'/em*. These results appear plausible but 
must be considered approximate. 


centres /cm? 
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The knowledge of the concentration of strontium and of free Sr++ in 
the crystal allows us to determine algo the association energy of the 
complex +-+© which is related to the equilibrium constant K of the 


association reaction through the equation 
ieeeet as COXA 14) Meme ES aac tanger (2) 


Using 7'=300°K, we obtain W=0-3 ev: this value is somewhat smaller 
than the theoretical value (0-39ev) (Bassani and Fumi, see following 
paper), as it should be (Lidiard). A cause of inaccuracy of our method for 
determining W is the uncertainty in the equilibrium temperature 7’. 


Fig. 2 
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The numbers of F and Z, centres/cm? as functions of the number of absorbed 


quanta/em® (Hg 54614) in a crystal of KCl containing 9 x 101? Srt++/em? 
coloured by X-rays. 


Figure 2 refers to the F> Z, conversion in a crystal coloured by X-rays. 
The general behaviour of the curves can be understood in relation with 
fig. 1 if one allows for the presence of holes. The observed very slow 
decrease of the Z, band is possibly due to a small ionization probability 
of the Z, centre under Z, light : this probability, however, appears 80 
small that the equilibrium to which it would lead in the F> Z, ev eclace. 
in additively coloured crystals would practically correspond to t a 
saturation of the available Z, traps. The apparent initial absence o 
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the Z, band may be rationalized if one considers that o, and oz, are 
approximately equal and that the irradiation with x-rays introduces 
vacancies (Esterman, Leivo and Stern 1949, Sakaguchi and Suita 1952). 

A description of the apparatus used and a detailed discussion of the 
results obtained will appear soon in J] Nuovo Cimento. 


The writers are indebted to Professor F. Seitz for useful correspondence. 
This research was supported by the Consiglio Nazionale delle Ricerche. 
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The Association Energy between Positive Divalent Impurities and 
Positive-Ion Vacancies in NaCl and KCl Crystals 


By F. Bassani* and F. G. Fumi1* 


Istituto di Scienze Fisiche, Universita di Milano, Italy 


[Received December 15, 1953] 


THE experimental data available at present on the properties of the 
alkali halides containing positive divalent impurities give a rather 
incoherent picture of the dependence of association between divalent 
impurities and positive-ion vacancies on the nature of the impurity and 
of the crystal. This situation suggests a theoretical calculation of the 
association energies, which represent the primary cause of the dependence : 
thus far only the association energy between Cd++ and a positive-ion 
vacancy in NaCl had been computed (Reitz and Gammel 1951). Lidiard 
(in press) has pointed out that the theoretical interaction energies are 
also needed to perform a correct statistical analysis of the degree of 
association including the long-range Coulomb interactions. 

The association energy between a positive divalent impurity ion and a 
positive-ion vacancy in a crystal of NaCl or KCl is given by the difference 
of the work necessary to create a positive-ion vacancy in the perfect 
lattice and the work necessary to create it in the position of next-nearest 


* Gruppo di Fisica dei Solidi, Sezione di Milano dell’Istituto Nazionale di 
Fisica Nucleare. 


Correspondence 229 


neighbour of the divalent ion. The work necessary to remove an ion 
from a lattice is conventionally taken to be equal to the negative of the 
average of the potential energies in the position of the ion before and 
after it has been removed. Methods for performing calculations of this 
kind have been developed by Mott and Littleton (1938) and improved 
by Reitz and Gammel (1951). We have followed the latter authors 
in taking for the repulsive potential energy the form of Born and Mayer 
(1932), which considers explicitly the radii of the interacting ions, taken 
as equal to the Goldschmidt values. 

Our results are given in the table. Their accuracy is probably of the 
order of 0-1 ev in the absolute values and probably better in the relative 
values. The association energy that we compute for Cd++ in NaCl 


Association Energies of Divalent Impurities and Positive-Ion 
Vacancies (ev) 


differs somewhat from the value reported by Reitz and Gammel (0-44 ev) 
mainly owing to our more accurate calculation of the energy needed to 
create a vacancy in the perfect lattice. 

The theoretical association energies of table 1 should not be compared 
directly with the values obtained from analyses of experimental data, 
which neglect the effect of long-range interactions, since such analyses 
necessarily lead to low values for the association energies (Lidiard). So 
our value for the association energy for Cd++ in NaCl is definitely 
larger than the value obtained by Etzel and Maurer (0-24 ev) (1950) from 
a simplified analysis of ionic conductivity isotherms, but it compares 
favourably with the value obtained by Lidiard (0-35 ev), who analysed the 
same data taking into account long-range interactions. 

The relative values of the computed association energies seem signifi- 
cant. They indicate an increase of the association energy with the ionic 
radius of the impurity, and actually a larger radius implies a stronger 
opposition to the inward displacement of the surrounding negative ions, 
which tends to compensate the excess positive charge of the impurity 
ion. They are also slightly smaller in KCl than in NaCl, and this agrees 
with the larger inward distortion in the immediate neighbourhood of the 
impurity caused by the contact between positive and negative ions 
occurring in KCl but not in NaCl, and also with the fact that the impurity 
ions are smaller than Kt but larger than Na’. 


Re See 
+ The value AE —0-24 ev is obtained from the data of Etzel and Maurer if 
one uses the appropriate equation L/F?=12 exp (4H/kT). 
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A detailed discussion of our calculation and of our results in relation 
to available experimental data will appear soon in J1 Nuovo Cimento. 


The writers are indebted to Dr. Lidiard for stimulating correspondence 


on the problem. 
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Coulomb Excitation of Levels in 1°F by Alpha-Particles 


By G. A. Jonzs and D. H. WILKINSON 
Cavendish Laboratory, Cambridge 


[Received December 29, 1953] 


CouLOMB excitation has been discussed theoretically (Mullin and Guth 
1951, Huby and Newns 1951, Ter-Martirosyan 1952, Alder and Winther 
1953) and demonstrated for protons bombarding heavy nuclei (McClelland 
and Goodman 1953, Huus and Zupanéié 1953). We have shown that 
the process occurs in the bombardment of a light nucleus with 
alpha-particles. 

A thick AIF, target was bombarded with alpha-particles and the 
gamma-rays detected in a crystal of Nal(Tl) (1 emx1 cm x1} mm) ; 
a typical spectrum is shown in fig. 1. The two large peaks represent 
gamma-rays of 112-++1 kev and 195+1-5 kev. Levels have been located 
at 113-+8 and 192+12 kev by Mileikowsky and Whaling (1952) using 
2INe(da)!°F and at 108 and 196 kev by Temmer and Heydenburg (1953) 
using 1°F(ax’)!9F, (They report that this reaction proceeds at high alpha- 
particle energies through compound nucleus formation but that there is 
evidence of Coulomb excitation at lower energies ; our own work is at 
far too low an energy to consider compound nucleus formation and 
proceeds entirely via Coulomb excitation.) The small bulge below the 
lower peak of fig. 1 is probably due chiefly to radiation back-scattered 
from material in the vicinity of the detector since its relative size 
depended on the proximity of the shields placed to absorb unwanted 
stray x-rays from the H.T. set. The relative intensity of the two strong 
lines of the spectrum is strongly dependent on the bombarding energy 
while the bulge remains in constant proportion to the total radiation ; 
it cannot therefore be ascribed to a cascade from the 195 kev level via 
that at 112 kev but it could well mask a slight amount of such a cascade. 

Figure 2 shows the excitation function for the higher energy radiation 
together with some theoretical curves for El and E2 excitation. The 
measurements were made at 90° to the alpha-particle beam ; the angular 
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distribution at 1-10 Mev is isotropic to within 10%. Either El or E2 
transitions give a good account of this excitation function (and also of 
that for the lower energy radiation). 

Rough estimates have been made of the absolute yields at 1-20 Mev : 
they are 3-5<10-1! and 1:4 10-19 gamma-rays per alpha-particle for 
the low and high-energy gamma-rays respectively, corresponding to 
(upward) E2 transition probabilities of 3x 104 and 8x 10® sec-! or to 
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Pulse spectrum from Nal(TI) due to gamma-rays following the bombardment 
of fluorine with alpha-particles of 1-10 Mev. 


El transition probabilities of 3 10° and 4x 101" sect, These represent 
| M |? values (in Weisskopf (1951) units) of about 0-2 and 2 for H2 
transitions or 0-002 and 0-04 for E1 transitions. 

It seems improbable that 1°F shows a change of parity at such low 
excitation and it also seems improbable that magnetic Coulomb 
excitation should be sufficiently effective here. If, then, we suppose 
the excitation to be E2 the two excited states must be (3+) or ($+). 
If the second excited state were (3+-) it should decay preferentially by 
M1 radiation to the 112 kev state and so it appears likely to be (3++). 
(If it were ($+) we should have | M |? <10-° for the M1 transition and 
this appears rather unlikely.) The (presumed) M1 decay of the second 
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excited state to the ground state competes favourably with the (presumed) 
M1 transition to the 112 kev state and so it appears probable that if a 
change of orbit is involved between the ground state and one excited 
state the same change does not occur for the other excited state since 
this situation would favour the cascade relative to the direct ground 
state transition from the 195 kev state, the former not requiring a change 
of orbit for the M1 transition. This makes it difficult to adopt the 
otherwise attractive assignment (2s,/.)* for the ground state of 1°F. 


Fig. 2 
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Excitation function for the 195 kev radiation. The three lines give theoretical 
functions for E2 excitation (Alder and Winther). The theory contains 
the parameter 4H/H where AE is the energy lost by the particle of 
energy #. In the dotted curve Z is the energy of the ingoing particle ; 
in the dashed curve # is its outgoing energy ; in the full curve Z is (as 
seems most realistic) its mean energy. The curves are normalized at 
1:25 Mev. The theoretical curve (Mullin and Guth) for El excitation 
differs inappreciably from the full curve for E2 excitation. 


If the ground state were (1d;/.)?2s,. as is frequently supposed it 
would be possible for the first excited state to be the (2+) lowest state 
of (Id5/.)°. The EK2 excitation should then have a strength (relative to 


the ‘ single particle ’ strength) of 
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following Lane and Radicati (1953) ; in this expression we have regarded 
the (0+) state of (1d;),)? as the principal parent of the ground state 
of °F. This compares with the experimental value of 0-2 Weisskopf 
units. The second excited state could then be perhaps a (2+) coupling 
of (28,,2)"1d5,.; alternatively the (3+) state of this mixed configuration 
may be this second excited state, its M1 transition to the 112 kev state 
being very powerfully forbidden by the need for two particles to change 
their configuration and by the double orbital change. In either event 
we should expect the E2 transition to be a strong one as observed. 
With such assignments these states would suffer considerable inter- 
configurational mixing and the magnetic moment of the ground state 
may be difficult to explain (B. H. Flowers, private communication). 
We may note in passing that the excitation energy of these states is 
much too low to associate them with the rotational levels of Bohr and 
Mottelson (1953). 

We intend to measure the decay lifetimes directly and to make a more 
careful search for cascade radiation. 
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XXVI. Notices of New Books and Periodicals received 


Electron Optics. By O. KtemperEeR. 2nd ed. (Cambridge: University Press.) 

Price 50s. net. 

Tux second edition of Dr. Klemperer’s monograph will be welcomed not only 
by those seeking a concise introduction to the various branches of Electron 
Optics, but also by those already familiar with the subject who nevertheless 
need an up-to-date guide to its voluminous publications. The extent to which 
the subject has developed since the first (1939) edition may be gauged from 
the fourfold increase in the numbers of pages and of references. 

The first seven chapters, dealing with electron lenses and their aberrations, 
may be seen to have their origins in the first edition but are greatly expanded 
to give account of the more important developments over the intervening 
years. The author draws attention, for instance, to suggestions which have 
been advanced for eliminating spherical aberration and to current trends in 
the design of B-ray spectrometers. 

The remaining chapters deal with topics hardly touched on in the first 
edition. Chapter 8 presents a useful account of the influence upon electron 
lenses of their intrinsic space-charge which brings out the importance of 
neutralization and of electrode structure. Chapter 9 sets out categories and 
properties of electron guns, Chapter 10 summarizes the properties of strip 
beams and their appropriate lenses, and Chapter 11 deals with the focusing 
properties of deflecting fields. The last chapter lists and describes very 
briefly the many applications of Electron Optics. 

It is, of course, impossible to survey so wide a field in so short a space 
without restricting the enquiry in some respect. In this case the author has 
chosen to describe in some detail experimental investigations but to report 
theoretical developments briefly, and where possible in non-mathematical 
terms, if at all. This sometimes leads to confusion. For instance, it is only 
in the absence of theoretical guidance that the divergence of gauze lenses 
seems ‘ remarkable ’ or doubt arises about the existence of spherical aberration 
in deflecting fields. It may also be regretted that the chapter on ray tracing 
did not describe a modern computational procedure and that the derivation 
of the erroneous formula for the focal length of thin electrostatic lenses should 
have been reproduced. 

It is probable that this book, which is well illustrated and thoroughly 
indexed, will prove most useful to research workers who come into contact 
with electron-optical equipment for it provides brief and lucid descriptions 
of the various components of electron-optical apparatus together with adequate 
references to relevant publications. BRASS: 


Introduction to Tensors, Spinors and Relativistic Wave-Equations. By 
E. M. Corson, Pa.D. (Blackie & Son Ltd.) [Pp. xiv-+-221.] Price 55s. 


Tas book dealing with a very abstract subject in a mercilessly abstract style 
will be wisely shunned by anyone primarily interested in the physical aspect of 
the problem of elementary particles. But even those of us who have drunk from 
the cup of abstract algebra more than perhaps is healthy for a physicist will, I 
am afraid, be disappointed by it. 

Professor Corson has made a gallant attempt to clear the jungle of literature 
on the subject to which he alludes in his Preface, and he has certainly disen- 
tangled to a large extent the intricate logical relationships between the various 
independent developments of the theory; but the result of his endeavour 
somehow falls short of the synthesis which might perhaps have been achieved by 
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a more decided handling of the issues involved. As it is the book is too difficult 
to be recommended as an introduction to the subject to research students, 
for example, and on the other hand the physicist likely to consult it will find 
only the better known parts of the subject discussed at great length, whereas the 
more subtle and fundamental investigations are only briefly summarized with 
an invitation to refer to the original papers. This, however, is a rather serious 
shortcoming, since the papers in question make use of little known aspects of 
modern algebra, for which they again refer to mathematical textbooks. The 
raison Wétre of a work like the one under review would just be to supply the 
physicist with the required mathematical background and to show more 
directly than can be done in original papers how it is applied to physical 
problems. 

On the credit side, it is only fair to praise the author for the neatness of hig 
handling of special points and for the great care that he has obviously devoted to 
rigour and accuracy of treatment. The publishers must be congratulated for a 
very high-class production. L. R. 


Investigations in Physics, No. 1: Ferroelectricity. By E. T. Jaynus. 
(Princeton: University Press; London: Geoffrey Cumberlege, Oxford 
University Press.) [Pp. 142.] Price 2$.00; 12s. 6d. 


Tas book is an expansion of a doctoral thesis and is concerned mainly with 
the theory of ferroelectricity. It opens with a brief account of the properties 
of Rochelle salt and KH,PO, and the theories put forward to account for their 
behaviour. The author then discusses at greater length the properties of 
BaTiO, and the theories relating it, in particular to his own electronic theory. 
There is then a brief account of the general thermodynamic properties of 
ferroelectrics, and the remaining one-third of the book is devoted to methods 
of evaluating internal fields in crystals and finding the conditions for 
spontaneous polarization. ; - 
The book suffers somewhat from being a mixture of a review and an original 
paper. Thus the space devoted to internal fields is excessive in a general 
account of ferroelectricity though as it is new material it is interesting to the 
specialist. Nevertheless the book gives a good and readable account of the 
present state of the theory of ferroelectrics. ee ed) 


Principles of Transistor Circuits. Edited by Ricarp F. Sama. (New York : 
John Wiley and Sons; London: Chapman and Hall.) Price 88s. net. 


Tas volume, which has been written by members of the Electronics Laboratory 
of the General Electric Company, Syracuse, N.Y., provides the first systematic 
account of transistor circuits. It is clearly written and contains the funda- 
mental theory needed for a full appreciation of the potentialities of the different 
types of transistor as circuit elements. Each chapter ends with a group of 
problems and a bibliography. The inclusion of problems in a monograph of 
this type is worthy of note at a time when interest in technical education is 
increasing ; the transistor and related devices will inevitably be discussed in 
the electronics courses of technical colleges and universities in the immediate 
future. The extensive bibliography provides a valuable guide to the literature 
of the subject. The book will need no recommendation to physicists and 
electronics engineers who are concerned with transistor developments. It will 
be widely read by directors of research in industry and service laboratories 
who are concerned with the assessment of the importance of the transistor in 
their particular fields ; this group will appreciate the balanced eg the 
subject which is presented. J.W.M. 
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The First Fifty Years, 1903-1953. THe Farapay Society. (Aberdeen 
University Press.) [Pp. 86.] Price 10s. 6d. 


In the words of the Editor: “‘ This little publication is not a set piece; it is 
more in the nature of a birthday-book in which people have been good enough 
to write interesting things about the Society now on the threshold of its fiftieth 
anniversary.”’ But he is to be congratulated upon it as a collection of short 
articles from fourteen of our leading physical chemists including one of its 
original members, Professor F.G. Donnan. The work also contains a number of 
interesting portraits of Past-Presidents. These, collected in one volume, 
help to emphasize the prestige and influence that the Society has acquired in 
the first half-century of its life. A.M. T. 
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